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The California Institute of Technology 
Electric Analog Computer 


Introduction. The California Institute of Technology Analog Computer 
has been developed to provide a general-purpose device having as wide a 
field of application as is considered practicable using the electric analog 
principle. The vast majority of problems in the general field of engineering 
analysis usually require no more than a one per cent accuracy in their solu- 
tions. With this lenient accuracy requirement, many simplifications can be 
made in the computer and its operation. Electrical resistance, inductance, 
capacitance, and transformer circuits can be used to simulate the linear 
terms of both algebraic and differential equations.'~7 This choice of circuits 
is in contrast to other electric analog computers, such as REAC (Reeves 
Electronic Analog Computer), which employ electronic amplifiers for all 
integration and differentiation operations. The electronic analog computers 
can provide higher accuracy and can be used as simulators on a one-to-one 
time base, but they are more restricted in their range of application than 
computers of the type discussed here. 

With electric analog computers only one independent variable can be 
continuously represented. However, by the use of finite-difference methods, 
such as those employed with digital computers, more independent variables 
can be handled. The computer at the California Institute of Technology is 
designed to handle linear partial differential equations with up to three 
independent variables and non-linear partial differential equations with 
two independent variables. Two basic circuit techniques are employed. In 
the first, a single independent variable is represented continuously as time 
on the computer. '5 In the other, discrete values of several independent 
variables are represented in a form suitable for the application of finite- 
difference methods; in this case, the analogous circuit becomes a steady- 
state a-c or d-c mesh.” Both techniques are applicable to both ordinary and 
partial differential equations. More complex electronic equipment is required 
in the application of the first principle. It involves devices to generate 
suitable functions of the independent and dependent variables, multipliers 
for multiplying any two variables together, devices for simulating arbitrary 
functions of a dependent variable, and amplifiers to provide the proper 
power output loads from these devices and to produce negative impedance 
terms and the unilateral or unsymmetrical terms of a matrix or servo- 
system. 

Description of the Computer. In July 1947 the Institute’s computer was 
placed in active service although it was only partially completed.'~* Since 
then it has been in full-time use for the solution of commercial engineering 
problems and for Institute research. The computer is best described by 
referring to Fig. 1 and Table*! which show its general layout and list the 
component parts. Its design embodies 110 separate sets of precision induc- 
tors, resistors, and capacitors, and 25 special transformers. These circuit 
elements are used to set up circuits for solving ordinary or partial differential 
equations. The inductors have a ““Q”’ of better than one hundred under all 
conditions of operation, and every element can be set to within one per cent. 
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General View of the California Institute of Technology Computer. 
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Photograph of Amplifiers and Arbitrary Function Generators. 


. Negative gain servo-amplifiers. 

. Positive gain (d-c) amplifiers. 

. Voltage and current limiters. 

. Dependent variable functions generator—Type I. 
. Dependent variable functions generator—Type II. 


There are, at present, 32 separate forcing functions for introducing (as 
known voltages, currents, or charge) the known functions of the independent 
variable. These may be steady-state, variable-frequency, sinusoidal func- 
tions; square wave transient functions; or perfectly arbitrary functions of 
time. As listed in the Table, 20 amplifiers are now available to represent nega- 
tive impedance terms and unilateral terms such as the unsymmetrical terms 
of a matrix. For simulating the non-linear terms of an equation, ten multi- 
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pliers are provided to multiply any two variables, and 11 arbitrary function 
elements are provided to form arbitrary functions of a dependent variable 
(see Fig. 2). Five current- and voltage-limiting devices are available for 
special non-linearities such as missile rudder stops or similar abrupt changes 
on a dependent variable. The dependent variables that constitute the 
desired solutions are in all cases simulated by circuit voltages, currents, or 
charges. For transient solutions, they are recorded or viewed on specially 
adapted cathode-ray oscilloscopes, while for steady-state or algebraic solu- 
tions, they are measured with vacuum-tube or dynamometer-type meters. 
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Fic. 3. Symbolization of Basic Computer Elements. 


The synchronous switches are important elements of the computer 
used to switch in the forcing functions, to remove energy from the circuits 
in order to restore them to their proper initial condition, and to make other 
sudden changes in circuit conditions. Thirty-six of these switches are now 
available. They are required only for transient solutions, which may be 
obtained in either of two ways. One way is to apply the transient just once 
and record it photographically or with some other memory device such as 
magnetic tape. The other way is to apply the transient cyclically (usually 
ten times a second) so that the solution appears as a steady image on a 
cathode-ray oscilloscope. In both techniques, each switch can be very 
accurately timed for opening and closing operations. Other important 
features of this analog computer are the plug-board arrangement and the 
metering circuits which permit quick connection and rapid metering of all 
solutions.® Proper voltages or currents can be readily inserted at the bound- 
aries or mesh points of circuits which simulate difference-equation approxi- 
mations to partial differential equations. 
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A special effort was made to design all elements of the computer so that 
additions can be readily made if sufficiently important problems arise that 
require greater facilities. 

Examples of Computer Techniques. The required specifications for the 
analog elements are best considered by first illustrating their application 
to specific simple problems and describing the analog techniques employed. 
The basic analog elements are listed in Fig. 3 in terms of the symbols em- 
ployed for circuit layouts. Typical analogous circuits for both ordinary and 
partial differential equations are given in Figs. 4 to 7. Figs. 4 to 6 describe 
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Fic. 4. Illustrative Analogy for Ordinary Differential Equation. 


techniques in which time continuously represents one independent variable 
on the computer. Fig. 7 illustrates the solution of partial differential equa- 
tions with a steady-state network, a very useful iterative method of solving 
many non-linear partial differential equations.” 

Fig. 4 illustrates the basic technique of forming an equation by one or 
both of KiRCHHOFF’s circuit laws, the use of inductors and resistors for 
linear differential terms, and the use of both the multiplier and non-linear 
function devices for non-linear terms. Fig. 5 shows the use of transformers 
for mutual bilateral terms and the isolator and amplifier for unilateral or 
non-symmetrical terms of a matrix of simultaneous equations. 








506 ELECTRIC ANALOG COMPUTER 

















anes ends aks victim Sent vans oes) ag) eas ea nates eat 
q 
L33 
pP—$ UHH — —— 
| ome | 
a3 
+A,! 
fa(t) Lee Ree 




















é 4 ] 

f, (ths ku Ge + ky SO + hyp SE @ = «ae 
dy2 dy! 

tlt) = Kag@2S + bgp Me. + ne + — --- 


a*y3 
tlt) © Ky OP + heave a 
falt)= -—--— 


Fic. 5. Illustrative Analogy for Simultaneous Equations. 


In Fig. 6, a partial differential equation is considered which requires 
multipliers and represents one of the independent variables continuously 
as time on the computer. This technique has limited application to non- 
linear partial differential equations, since too many special electronic ele- 
ments would be required for general application. 
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COMPLETE TIME ANALOGY FOR PARTIAL DIFFERENTIAL EQUATION OF HEAT TRANSFER. 
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Fic. 6. Complete Time Analogy for Partial Differential Equation of Heat Transfer. 
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The general technique illustrated in Fig. 7 for non-linear partial differ- 
ential equations has so far been found to have a surprisingly wide field of 
application to equations containing terms up to the fourth order of differ- 
entiation and to pairs of simultaneous equations. Since circuit element 
settings and inserted mesh currents can be changed rapidly, iterative 
methods of solution can be set up in a very short time. 
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As illustrated in the examples, RLC circuit elements are usually em- 
ployed for all linear self-impedance terms. They are also employed for 
bilateral mutual terms, with junction points between the meshes simulating 
the simultaneous equations, and for transfer functions. Special transformers® 
are used, wherever possible, for bilateral mutuals requiring circuit isolation 
and, in conjunction with amplifiers, to provide a four-terminal unilateral 
element. The RLC circuit elements set the speed of computation. In order 
to have the elements behave as near-ideal resistors, inductors, and capaci- 
tors, the natural frequencies of the analog circuits must lie within certain 
limits. The normal operating frequency range was chosen to be from ten 
cycles to one or two kilocycles per second." * Over this wide frequency range, 
transformers have certain limitations that require the occasional use of 
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‘“Ssolators” to make them operate properly from zero frequency (d-c) up 
to the limit of one or two kilocycles per second (see Fig. 3). 

The circuit impedance base must be kept low and the power level 
reasonably high to minimize stray coupling and stray pickup problems 
and enable the quick set-up of a wide variety of circuits. Unintentional 
impedances, such as the output impedance of a voltage amplifier, must be 
kept below a few ohms in series elements and above several hundred thou- 
sand ohms (or in some cases 10 or 20 megohms) in the case of shunt im- 
pedances. It is difficult to keep the ratio of the series-to-magnetizing 
impedance of.a transformer sufficiently low over such a wide operating 
frequency range. This, together with the tendency of transformer coupled 
feedback circuits to ‘‘motor-boat,’’ constitute the major difficulties en- 
countered in the use of transformers. The ‘‘motor-boating”’ is only important 
when transient solutions are being obtained by repeating cyclically the 
forcing functions at the basic 10-cycle-per-second repetition rate in order 
to obtain a standing-wave solution on an oscilloscope. It can be overcome 
by applying the transient just once and recording the solution photograph- 
ically or with a magnetic-tape recorder. 

All of the electronic equipment is required to have less than one per cent 
amplitude distortion throughout the operating range. Furthermore, the 
frequency response of the equipment must not vary more than .1 decibel 
throughout the frequency range of the computer. The response is considered 
to be satisfactory if the phase shift is less than one degree throughout the 
entire frequency range. However, since these devices are frequently used in 
cascade and feedback circuits, the phase and amplitude characteristics 
must be carefully controlled to eliminate oscillations at both high and low 
frequencies. To prevent low frequency ‘‘motor-boating,” all electronic 
equipment except the initial forcing-function generators was designed for 
operation from zero frequency (d-c) to the upper frequency limit of the 
computer. The high-frequency limit of all equipment was extended as far 
as practicable to make it easier to eliminate high-frequency oscillations. 

The electronic elements discussed above have been described in detail.” 1” 
It might be of interest, however, to discuss briefly the circuit elements for 
non-linear functions. 

Generator of an Arbitrary Function of the Dependent Variable. One of 
the generators of non-linear functions of the dependent variable (Type 1) 
uses crystal diodes and batteries in the manner shown in Fig. 8. The ad- 
justable resistors are relatively large compared to the forward impedance 
of the crystals, so that the variations of this impedance with current do 
not affect the over-all impedance appreciably. If the switches are all thrown 
to the left in Fig. 8, each parallel branch will not conduct current until the 
terminal voltage rises above the battery voltage of that branch. The result- 
ing current-voltage characteristic is composed of a series of straight lines, 
and the slope of successive segments decreases as the voltage rises. Such a 
characteristic is shown in curves A of Fig. 8. The two curves labeled A 
are both of the same characteristic, but they indicate that both the positive 
and negative parts of a function may have entirely different characteristics 
if desired. 

If the switches are thrown to the right, a different characteristic results. 
Each branch produces a circulating current which passes through Ro. 
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The battery voltage E» gives zero terminal voltage when no external current 
enters the terminals. For a given applied voltage, the current is determined 
by the effective resistance of all the branches in parallel. As soon as the 
terminal voltage exceeds the lowest branch battery voltage, the current in 
that branch ceases to flow, and the effective impedance of the device is 
raised. When the terminal voltage exceeds the highest battery voltage, the 
impedance becomes equal to Ro, the maximum value. The resulting current- 
voltage characteristic is one of constantly increasing slope, as shown in 
curve B. By combining branches of each type, a characteristic may be 
achieved whose slope is alternately increasing and decreasing, provided it 
always remains positive. Curve C, Fig. 8, is typical. 
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Fic. 8. Schematic Diagram—Dependent Variable Device Type 1. 


It has been assumed that the voltage applied to the device is propor- 
tional to the dependent variable and that the current is proportional to an 
arbitrary function. If, however, the generator is driven by a current source, 
then the voltage across the device represents the arbitrary function. In this 
case, the type of characteristic will be reversed because of the interchange 
of axes (see Fig. 8). 

The generator which is used with the computer has 22 parallel branches, 
11 for positive voltages and 11 for negative voltages. A current-voltage 
characteristic composed of 22 straight-line segments can then be obtained. 
Such a characteristic can be used to approximate many typical functions 
with surprisingly good accuracy. Simpler and more specialized versions of 
the generator have also been developed in the form of abrupt current and 
voltage limiters. 

The other arbitrary function generator (Type 2) is based upon a prin- 
ciple developed elsewhere.'*®° It employs a cathode-ray tube and photo- 
electric cell. With no signal, the electron beam is directed to one side of the 
tube. The photo-cell output is amplified and applied to the plates with the 
proper polarity to drive the beam to the other side. If an opaque template 
is placed on the screen of the tube, as shown in Fig. 9, the beam spot will 
not appear until the light reaching the photo-tube produces a large enough 
signal to deflect the beam above the edge of the template. 
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The generator has several notable limitations. There is some distortion 
for which compensation is difficult, since the spot is not infinitely small 
and because the amplifier gain, the spot intensity, the height of the pattern, 
and the quiescent position all may vary from time to time. In addition, the 
halo on the screen contributes an appreciable signal, which may cause 
clipping of sharp, high peaks of the template, since the halo signal can be 
picked up on both sides of the peak with the spot completely hidden. Since 
the actual spot is on the rear of the cathode-ray screen while the observed 
pattern is on the front, there will be a third source of error in the resulting 
parallax. Careful proportioning of over-all height and width of the pattern, 
location of quiescent position, and adjustment of gain minimize these 
difficulties, but a definite accuracy limitation exists. The inability of the 
beam to rise or fall abruptly introduces another limitation. Thus, if the 
dependent variable passes rapidly through a region where the arbitrary 
function changes abruptly, the output will lag behind the ideal curve. This 
effect can be minimized by reducing the vertical scale of the template. Since 
the time lag is very small compared to the period of the highest frequency 
used on this computer, the effect on solutions is negligible. Despite these 
inherent limitations, it has been possible to develop a device that will form 
an arbitrary function of a dependent variable to within one per cent of the 
function’s crest magnitude. 
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Fic. 9. Schematic Diagram of Arbitrary Forcing Function of Dependent Variable—Type 2. 


Multiplier. The multiplication of two arbitrary voltages can be accom- 
plished in many ways. The method which has been adopted for the CIT 
Computer uses double modulation followed by detection.? A balanced 
modulator produces sidebands by modulating a carrier with one of the 
product terms. The resulting voltage is of the form E, cos (wt). This voltage 
is used as a carrier for a second modulator which gives a sideband output, 
E,E, cos (wt). The carrier voltage is, of course, suppressed through use of 
a balanced modulator in both cases. Addition of a carrier voltage permits 
detection, with the output voltage proportional to the product, E,E. The 
d-c voltage resulting from detection of the added carrier is balanced by 
detection of another carrier signal and subtraction of the outputs. For con- 
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venience, a varistor bridge is used for the first modulator and vacuum tubes 
for the second modulator. Detection is accomplished with germanium diodes, 
thus minimizing the power supply requirements of the unit. Isolation be- 
tween input and output is automatically supplied by the high-frequency 
transformers. One of the input voltages is not isolated from the system 
ground. If such isolation is required for both inputs, an additional isolator 
must be used. Care must be taken to minimize phase shift in the sidebands 
since this results in a phase shift of the detected signal. The first modulator, 
with untuned transformers, performs very well in this respect. Since tuned 
circuits are required in the second modulator for several reasons, the phase 
shift is minimized by over-coupling of the transformer to the plate circuit 
of the modulator. In the case of the multiplier, the accuracy is much better 
than the over-all one per cent requirement. The groups of multipliers and 
d-c isolators in the computer are supplied by a common high-frequency 
generator, thus reducing the space required for each unit. 

Present Fields of Application. Although the computer has been in active 
service for almost two years and has been applied to a wide range of prob- 
lems in many fields of engineering and science, it is not yet possible to predict 
its complete range of usefulness. It may be more appropriate to merely list 
some of the more important types of problems for which it has been used. 

In the field of servomechanisms and automatic control,*!*%1!4!7 the 
computer makes possible the rapid solution of control equations describing 
even the more comprehensive cases of missiles or autopilots with several 
degrees of freedom, and non-linearities in both the control system and the 
aerodynamic equations. 

General analyses have been made of many simple and complex linear 
and non-linear ordinary differential equations pertaining to mechanical 
vibration analysis.* "6 Equations have been solved occurring in such fields 
as packaging of delicate apparatus, the investigation of the design and per- 
formance of turbine generators, short-circuit torques, airplane landing shock, 
acceleration through resonance, and earthquake shock effects on buildings. 

A practica! method has now been developed for solving accurately and 
rapidly the more general partial differential equations of elastic beams and 
plates.'* This facilitates the solution of the general types of transient beam 
and plate problems, as well as static loading and steady-state vibration 
problems. 

The computer has been used for the solution of several linear thermal 
conductivity problems involving time and two space variables, and of non- 
linear conductivity problems involving time and one space variable. Com- 
plex boundary conditions and heat generation within the media can be 
readily handled. 

Other fields for which the California Institute of Technology Electric 
Analog Computer has so far been used include compressible fluid flow 
studies, complex electric circuit analysis, aircraft landing stability anal- 
ysis, basic mechanical system backlash and non-linear friction analysis, 
rideability analysis of automobiles and railroad cars, and many others. 

G. D. McCann 


California Institute of Technology 
Pasadena, California 
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acceleration through critical speeds,” Jn. Appl. Mech., v. 16, 1949 (in press). 

17G, D. McCann, C. H. Witts, & B. N. Locantut, “Application of the Cal. Tech. 
Electric Analog Computer to non-linear mechanisms and servomechanisms,” AIEE, 
Trans., 1949 (in press). Techn. Paper, no. 49-165. 

3G. D. McCann & R. H. MacNEAt, “Beam vibration analysis with the electric 
analog computer,”’ Jn. Appl. Mech., v. 16, 1949, p. A (in press). 

19D. M. MacKay, “A high-speed electronic function generator,’’ Nature, v. 159, 
Mar. 22, 1947, p. 406. 

20D. J. MYNALL, “Electrical analogue computing,’’ Electronic Engin., London, v. 19, 
1947, p. 178-180, 214-217, 254-262, 283-285. 
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TABLE 


Basic Elements of the California Institute of Technology 
Electric Analogue Computer 


A. Ten Passive Circuit Elements Cabinets 


1. 100 high “Q” inductors 0 to 1.0 henrys in 0.001 henry steps 
2. 20 high “Q” inductors 0 to 1.50 henrys in 0.01 henry steps 
3. 100 resistors 0 to 7,000 ohms in 2 ohm steps 

4. 20 resistors 0 to 500,000 ohms in 200 ohm steps 
5. 60 capacitors 0 to 1.0 mfd in 0.01 mfd steps 

6. 40 capacitors 0 to 4.0 mfd in 0.01 mfd steps 

7. 10 capacitors 0 to 80 mfd in 10 mfd steps 

8. 25 special transformers* 


B. Five Plugboards 


1. Plugboards for all elements of two circuit elements cabinets 
2. Receptacle patch-board for 500 element connections and 25 main busses 


C. Two Forcing Functions Cabinets 


1. 10 sinusoidal variable-frequency and phase-position-voltage sources (0 to 50 volts, 
one ampere, one ohm) 
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2. 12 separate square wave voltage sources (0 to 12 volts, or one 72 volt source in 2 volt 
steps 

3. 5 arbitrary functions of independent variable 

4. 5 arbitrary current sources (controllable from any voltage source) 


D. Two Metering Desks 


1. 3 dynamometer-type meters and 2 vacuum tube meters for metering sinusoidal 
current, voltage, power, etc. 
3 cathode-ray oscillographs for transient solution (current, voltage, or charge) 
2 three-wire metering circuits (and selector system) for metering all passive circuit 
elements and source forcing functions 
2 manual metering selector circuits for 16 of 25 main busses 
Circuits for connecting 36 synchronous switches to any of 26 main busses 


ae WN 


E. Amplifier Cabinets (portable) 


1. 5 negative gain (0 to 50) d-c amplifiers (15 ohm impedance), each having an RC 
time-delay circuit for servo-problems 


2. 15 positive gain (0 to 100) d-c amplifiers with input and one ohm output impedance 
3. 10 d-c isolators and preamplifiers 


F. Multiplier Cabinets (portable) 
1. 10 multipliers with isolated inputs* 


G. Arbitrary Functions of Dependent Variable (portable) 


1. 3 voltage limiters 

2. 2 current limiters 

3. 1 arbitrary function of dependent variable—Type I 
10 arbitrary functions of dependent variable—Type II 


* Additional elements are contemplated. The Table lists only elements of the com- 
puter in its present operating form. 
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658[A-D, Q, R].—ItTa.y, IstiTuUTO IDROGRAFICO DELLA R. MARINA, Tavole 
Logaritmiche a Cinque Cifre Decimali. Sixth ed. Genoa, 1941, xxii, 512. 
19.7 X 27.8 cm. Full cloth, and gilt titles. 


Contents: T. I (p. 1-57).—Mantissae of logarithms and cologarithms of numbers from 
1 to 10 080; T. I (p. 60-179)—Logarithms of the absolute values of the trigonometric func- 
tions 0(1’’)2°; T. III (p. 182-361)—Logarithms of the absolute values of the trigonometric 
functions 0(15’’)90°; T. IV (p. 364~-382)—Addition and subtraction logarithms; T. V (p. 
384—428)—-V alues of the natural trigonometric functions 0(1’)90° or 0(4*)6*; T. VI-XIII 
(430-464)—Conversion of common logarithms tq natural, and vice versa; n?, n*, nt, ni, In n, 
1000 n-', xn, 42n®, n = 0(1)1020; table of arcs, altitude of sector, chord for 0(1°)180°, also 
arcs corresponding to 0(1’)60’, 0(1’)60”; conversion of centesimal grades, minutes and 
seconds into sexagesimals; conversion of arcs into time and vice versa; conversion of sex- 
agesimal minutes and seconds into decimals of degrees, and vice versa; geodetic elements 
relative to the international ellipsoid ; values of international normal gravity, g, for  (geogr. 
latitude) = 0(15’)90°. Constants and units of measure (p. 467-480). Formulae (p. 483-512) 
—algebra, plane analytic geometry, differential and integral calculus, plane and spherical 
trigonometry, series, interpolation, astronomy, geodesy, least squares. 


EpitroriaL NotE—The first edition of this work was published in 1913; the second, 
1916; third, 1922; fourth, 1930; fifth, 1937. There was a “nuova edizione, ridotta,” 128 p., in 
1939, and a copy may be seen in Library of Congress. 
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659[A-E, J—L].—Joun Bortuwick DALE, Five-figure Tables of Mathematical 
Functions comprising Tables of Logarithms, Powers of Numbers, Trigo- 
nometric, Elliptic, and other Transcendental Functions. Second edition, 
London, Edward Arnold & Co., 1949, viii, 121 p. 13 X 21.6 cm. 6 shill- 
ings. New York, Longmans, Green & Co., $1.50. 


This popular work, first published in 1903 (xv, 92 p.) was reprinted in 1905, 1908, 1909, 
1911, 1913, 1917, 1919, 1922, 1927, 1932, 1937, 1942, 1943, 1945. It was prepared with the 
idea of meeting the requirements of workers in Physical Science and Applied Mathematics. 
For the new edition the type has been reset, two new tables have been added (IA (p. 14-18), 
logarithms of numbers from 1000 to 2999, and XVIIA (p. 88-90) logarithms of hyperbolic 
functions), and other tables have been rearranged in place and in forms more suitable for 
ready use. It would seem as if many errors had been made in resetting the tables; here are 
five random illustrations: In 4 = 1.38629 (instead of 1.38269); 3.98 = 1.5848 (instead of 
1.5484); 4.55! = 1.6571 (instead of 1.6751); 5.49! = 1.7641 (instead of 1.7461); p. 112, 
last col. x = .87, .83424 (instead of .83244). 

The Tables are still numbered I-XX XIX: I-VII (p. 10-53), logs of numbers, antilogs, 
squares, square roots, cubes, cube roots, reciprocals; VIII-XIV (p. 54-81), the natural 
trigonometric functions and their logs, circular or radian measure; XV-—XVIII (p. 82-91), 
natural logs—to base e, nat. log 10*", exponential and hyperbolic functions, log e*; XIX 
(p. 92), gd0; page 93 is blank; XX-XXII (p. 94-102), F(0, ¢), E(0, ¢), complete elliptic 
functions (the rounding error of F,(88°) in the final digit, 1 instead of 2, MTAC, v. 3, p. 262, 
still persists); XXIII-XXV (p. 103-105), gamma function, logs of factorials, zonal surface 
harmonics, P,(x), m = 1(1)7, x = 0(.01)1; XXVI-XXVII (p. 106-107), Jo(x), Ji(x), 
x = 0(.1)15, Ja(x), m = 0(1)12, x = 0(1)20 (the error in Jo(5.9) in the first ed. has been 
corrected); XXVIII (p. 108-109), I,(x), m = 0(1)11, x = 0(.2)5, Lodge, 1889; XXIX 
(p. 110), J, = 0, first 9 roots m = 0(1)5 (here are 27 end-figure errors of 1 to 25 units) ; first 
6 roots of Jo(x) and Ji(x), jo,./m and j:,./" (Stokes 1851 tables); XXX (p. 110), Jn(x), 
n=v+4, +2n = 1(2)11; XXXI (p. 111), Erf(x); XXXII (p. 112-113), Ei(+x), Ci(x), 
Si(x); XXXIII-XXXIV (p. 114-116), binomial coefficients for interpolation by differences, 
and proportional parts; XXXV-XXXVII (p. 117), reciprocals of factorials, Bernoulli 
numbers (B:-B20); conversion of time and angular measure; XX XVIII (p. 118), decimal 
equivalents, fractions of Julian year of 365} days; XX XIX (p. 119-121), numbers used in 
calculations, Euler’s numbers, Euler’s constant; 10D values, for » = 1(1)20, of S, = 1-* 
+2-"*+3-"---,s, = 1" —2"*+4+3"%—.:--, T, = 1-*°+3-*°+5-"*+---; English 
and metric equivalents. 

Dale does not give the source of any of his tables. 

In conclusion it may be noted that in the first ed. [1913, “sixth impression’’], p. 27, 
arg. 993, for 3.9766, read 9.9766; p. 91, 1. 20, for 3.0205, read 3.0255. 

_ = om © 


660[B-E, M, N].—Anpré SaInTE-LaGtE, Mathématiques, p. A30-A110, in 
Techniques de l'Ingénieur, directeur général de la rédaction C. MONTEIL. 
I: Généralités. Paris, Techniques de |’ Ingénieur, 1947, 24.3 X 29.4 cm. 


There are various tables of powers, quarter squares, interest, logs, antilogs, natural and 
hyperbolic trigonometric functions, etc., p. A30-A60, 43 p.; table of indefinite integrals, 
p. A80-5 to A80-14. Nomography, p. A110-12 to A110-18. 


661[C, D, K].—Anpr& Danyjon, Tables des Fonctions Trigonométriques. 
Valeurs naturelles a 6 décimales de centidme en centiéme du degré nona- 
gésimal. Paris, Librairie Hachette, 1948, ii, 126 p. 18.4 K 27.2 cm. Stiff 
cover. 850 francs. 


This is the third publication, within a year, in which are tables of natural trigonometric 
functions in the nonagesimal system with the arguments taken at interval 0°.01. Comrie’s 
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tables (RMT 631), and NBSCL (RMT 662) are the others. The further publication of tables 
using minutes and seconds, instead of decimals of a degree, is likely to become increasingly 
rare. Thus we are now harking back to tables such as nearly 350 years ago HENRY Briccs 
thought were the most useful he could prepare. 

Here are 6D tables (p. 6-95), with A, of sin, csc, tan, cot, sec, and cos. In connection 
with sines and cosines linear interpolation can be used throughout, and in the case of tan- 
gents and secants from 0 to 75°.8; the method for dealing with later angles is indicated. 
If one is content with 6S, linear interpolation is legitimate up to 87°. On p. 5 are interpola- 
tion tables for trigonometric functions, including S and T. There is also a loose card for 
11(1)200(100)900 in a proportional parts table. 

The present tables were prepared from a 10D manuscript. The author states that the 
proofs were carefully compared with the tables in Trigonometria Britannica, 1633, of Briggs 
and for each 0°.05 with ANpover, Nouvelles Tables Fondamentales (1915-1918). 

On p. 96-99 are tables for the conversion of (a) sexagesimal fractions of a degree into 
decimal fractions; (b) decimal fractions of a degree into sexagesimal fractions; (c) hours into 
degrees and decimal fractions of a degree; (d) degrees into hours. A 6D table of logs 1000(1)- 
9999 occupies p. 100-117, and various mathematical formulae and constants, p. 119-125. 
On p. 126 are reprinted two statistical tables, Normal distribution of Laplace-Gauss, and 
Distribution of x* from R. A. Fisher & F. Yates, Statistical Tables, third ed., 1948. 

The work presents a neat appearance. The one-column “Directions for the use of the 
tables” is in four languages: French, Spanish, English and German. 


R. C. A. 


662[D].—NBSCL, Table of Sines and Cosines to Fifteen Decimal Places at 
Hundredths of a Degree. (NBS Applied Math. Series, no. 5.) Washington, 
Govt. Printing Office, 1949. viii, 95 p. For sale by the Superintendent of 
Documents, Washington, D. C. 40 cents. 19.8 K 25.9 cm. 


We have here an original valuable table of sines and cosines to 15D, &, throughout the 
quadrant, at interval 0°.01. A preliminary manuscript to 18D was checked by differencing 
and then rounded to 15D and typewritten. The final manuscript was again checked by 
differencing. 

The introduction by H. E. Sauzer includes a discussion of direct and inverse interpola- 
tion, with examples, and the usual interpolation tables, of E(p) and F(p), are given at the 
end of the volume, p. 94-95. 

The appearance of the volume would have been decidedly improved by the omission 
of the tautological headings, ‘“‘Table of Sines and Cosines,”” on 90 pages. 

It is more than 300 years ago since tables such as these were first computed and pub- 
lished. About 1601 Henry Briccs (1561-1631) finished the calculation of a remarkable series 
of tables at interval 0°.01, throughout the quadrant, including those of sines and cosines to 
15D, A. These tables were posthumously published by HENRY GELLIBRAND (1633). We have 
already indicated (MTAC, v. 1, p. 129-130) that Briggs was led to the decimal choice of 
such parameter by perusal of a publication of Vieta. In the first part of this 1633 volume 
(p. 43-44) there is also a table of sines for each 0°.625 in the quadrant to 19D, A. 

Many years earlier than Briggs RuEticus (died 1576) computed an extraordinary 
table of sines for each 10” of the quadrant, 15D, A* or A*; compare N109. Hence we may say 
that this table includes the values of sines and cosines to each 0°.125 of the quadrant. 
This table was first published by Pit1scus* in 1613 along with a Rheticus table of sines at 
interval 1”, from 0 to 1° and from 89° to 90°, 15D, A’, and hence giving values at interval 
0°.0625 in this range. 

In the work under review there is no reference to these Rheticus tables or to the AN- 
DOYER? tables (1915) of sines and cosines covering the same range, 15D, A, and hence includ- 
ing values to each 0°.125 of the quadrant. 

There are three further tables of sines and cosines to at least 15D, as follows: by HERR- 
MANN‘ (1848), [1°(1°)89°; 30D], and reprinted, p. 92-93 in the work under review; by J. T. 





516 RECENT MATHEMATICAL TABLES 


Peters® (1911), [0(10’)45°; 21D], and therefore including a table at interval 0°.5; and 
SPENCELEY & SPENCELEY® (1947), [1°(1°)90°; 25D]. 

On checking the first 50 values each of sines and cosines of the NBSCL tabie and that 
of Briggs, it was found that among the 50 sine values of the latter 42 differed from the former 
only in the fifteenth place, and by 1 to 4 units (20 unit errors, 5 2-unit errors, 5 3-unit errors, 
and 12 4-unit errors). Among the values of the 50 cosines there were only 11 unit errors in 
the fifteenth decimal place. 


R. C. A. 


1H. Briccs & H. GELLIBRAND, 7 rigonometria Britannica, Gouda, 1633. 

2G. J. Rueticus, ed. by B. Pitiscus, Thesaurus Mathematicus, Frankfort, 1613. 

3H. ANDOYER, Nouvelles Tables T: rigonométriques Fondamentales, v. 1, Paris, 1915. 

4 HERRMANN, “Bestimmung der trigonometrischen Functionen aus den Winkeln und 
der Winkel aus den Functionen, bis zu einer beliebigen Grenze der Genauigkeit,’’ Akad. d. 
Wissen., Vienna, Math. naturw. Kl., Sitzb., v. 1, “second unchanged ed.,” 1848, p. 477-478; 
a similar table of tangents and cotangents is given on p. 479-480. In this v. the pages are 
numbered consecutively from the beginning to the end. These page numbers agree with 
those given in the Royal Soc. Cat. Sci. Papers. Yet the references given by Mr. SALZER, p. v, 
are entirely different, namely: pt. IV, p. 176-177, and 178-179. It would therefore seem as 
if the five parts of the volume were originally issued with separate paging but when assembled 
in a volume were paged consecutively. 

5 J. T. Peters, “Einundzwanzigstellige Werte der Funktionen Sinus und Cosinus zur 
genauen Berechnung von zwanzigstelligen Werten samtlicher trigonometrischen Funk- 
tionen eines beliebigen Arguments sowie ihrer Logarithmen,’’ Akad. d. Wissen., Berlin, 
Phys.-math. Cl., Abh., 1911, p. 12-18. 

6G. W. & R. M. SPENCELEY, Smithsonian Elliptic Functions Tables. Washington, 1947. 


663[D, H, M].—Harry D. Huskey, “On the precision of a certain procedure 
of numerical integration’; Appendix: D. R. HARTREE, “Note on sys- 
tematic rounding-off errors in numerical integration,” NBS, Jn. of 
Research, v. 42, 1949, p. 57-62. 19.7 X 26 cm. 


Summary: An example of numerical integration (illustrated by three trigonometric 
tables) is given and shows very systematic effects in the less significant digits. This lack of 
randomness gives rounding-off errors that exceed the predicted standard deviation by a 
factor of three. 

The example considered shows that systematic rounding-off errors can occur in numeri- 
cal integration, irrespective of the number of digits kept in the contributions to the integral. 
In the Appendix this phenomenon is examined, and criteria are set up to detect the cases 
in which it may arise to a serious extent. 


664[D, I].—J @RGEN RyBNER & K. STEENBERG SGRENSEN, Tabel til Brug 
ved Addition af Komplekse Tal. Table for use in the Addition of Complex 
Numbers. Copenhagen, Jul. Gjellerups Forlag, 1948, xiv, 95 p. 23.6 31.5 
cm. 20 Danish crowns. American agent: Scandinavian Book Service, 
P.O. Box 99, Audobon Sta., New York 32. $5.50. 


The addition of two complex numbers can be carried out by means of a table by divid- 
ing by the number having the larger numerical value (in KENNELLY's notation /¢ = e’? 


=cos ¢ + isin ¢): 
23=2a+22= ri/ or + ro/ oo = ri/difl + (r2/ri)/b2 — $1), 
to S11 = 1i/gi(1 + 7/$) = ri/o1°R/a = nR/¢1 + a. 


Subtraction is reduced to addition by changing the negative sign of z2 into an alteration 
of its phase angle by +180°: 


23> ri/or - r2/2, = ri/or + r2/d2 + 180°. 
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The table represents the function 
R/a =-1+ r/o 
giving R(S5D) and a(3D) as functions of r and ¢, for r = 0(.01)1, @ = 0(1°)180°. 
=1+7r+2rcos¢, sina=rsing/R, tan($¢—a) = [(1 —7)/(1 + 1) ]tan }¢. 
The computations were started by calculating R and a directly from these formulae at 
intervals .05 for r, and 5° for @. These calculations were carried out by means of the seven- 
place trigonometrical tables of PETERS,’ R being determined to 7D and a to SD, that is, 
in both cases two more places than given in the final table. The other values of the table were 
found by interpolation. 
Pages 2-93 are devoted to tabulation of R and a for r = 0(.01)1, ¢ = 0(1°)180°. Differ- 


ences d,, dg are given throughout. At the bottoms of the R and a pages are the interpolation 
formulae: 


R(ro + xh, oo + yk) = Roo + xd, + ydg + xydrg? + B”(x)d,? + B’(y)dgg? 
a(ro+ xh, oo + yk) = coo + xd + ydg + xydrg? + BY’ (x)d,? + B’ (y)dgg? 
and on each page are the ranges of values of d,,*, d,;*, dgg?. On p. 95 is a table of interpolation 
coefficients B’’(u) = — 4u(1 — u), u = x or y. 

The introductory pages include a number of worked out examples illustrating the table’s 
use; for example: the values of R and a are found when r = .5082, @ = 136°.59. 

It is obvious that great care was taken to attain accuracy in preparation of the volume. 

In an earlier work by Professor Rybner,? which we reviewed, MTAC, v. 3, p. 174-175, 
there was a nomogram for R/a =1i+ r/o but the present table was necessary because it 
was not found possible to construct the corresponding nomograms, with an accuracy suffi- 
cient for practical calculations. 

A small table of R/a =1+ r/¢ was given by F. Empe* but the angles are measured in 
fractions of a right angle, 0(.05)2, which is inconvenient, since this unit is not commonly 
used. 

i 
1J. T. Peters, Siebenstellige Werte der trigonometrischen Funktionen von Tausendstel zu 
Tausendstel des Grades. Berlin, 1938. 
27. RyBNER, Nomograms of Complex Hyperbolic Functions. Copenhagen, 1947. 
*F, EmMDE, Tables of Elementary Functions. Tafeln elementarer Funktionen. Leipzig, 1940. 


American reprint 1945, p. 32-33. See MTAC, v. 1, p. 384-385. The table is also given in 
JAHNKE & Empe, Tables of Functions, e.g. 1945, p. is-19 of the second part. 


665[E].—HARVARD UNIVERSITY, COMPUTATION LABORATORY, Annals, v. 20, 
Tables of Inverse Hyperbolic Functions, By the staff of the Laboratory, 
H. H. Aiken, director. Cambridge, Mass., Harvard University Press, 
1949, xx, 290 p. 19.5 X 26.7 cm. $10.00. Offset print. 


These tables give 9D values with first and second forward differences of the inverse 
hyperbolic sine, cosine, and tangent. More precisely, T. I (p. 3-47), from which T. I and 
T. III were derived, gives tanh x for 0 < x < 1. For x < .9 the values of the argument are 


x = 0(.001).5(.0005).75(.0002).9. 
Only 39 percent of the table is devoted to 90 percent of the unit interval. The singularity of 


tanh x at x = 1 forced the tabulation for .9 < x < 1 to proceed by finer and finer intervals 


in order to preserve the adequacy of second order interpolation. The final entry is for 
x = .99999. 


T. II (p. 51-66) is straightforward and gives sinh x for 
x = 0(.002)3(.005)3.495. 


T. III (p. 69-103) gives cosh x for 1 < x < 3.5. Here trouble is encountered at x = 1 
so that one-half of the table has to be devoted to 6 percent of the entire range of x. Beyond 
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x = 1.15 the argument values are 
x = 1.15(.0005)1.4(.001)1.8(.002)3.498. 


Finally T. IV (p. 107-290) gives sinh x and cosh x for 3.5 < x < 22980 at various 
intervals ranging from .005 to 20 keeping the second differences under about 2000 units 
in the 9th decimal place. 

The accuracy of the whole volume is to within .9 units in the 9th decimal. The computa- 
tion time was only 12 days, of 24 hours each. This represents an increase in speed of the 
Automatic Sequence Controlled Calculator due to an improved printing mechanism. 

The fact that the table was produced by automatic methods accounts for a number of its 
features which would have puzzled table makers of a decade ago. For example T. IV has 
the luxurious feature of tabulating two functions so nearly alike that their values differ by 
only a few hundred units in the last place of decimals and their second differences are prac- 
tically indistinguishable. This goes on for more than 150 pages and would be an unthinkable 
waste of time and paper, not to mention an added source of typographical error, had the 
table been hand-made. In fact, by past standards an extensive table of these three functions 
has no doubt been considered unworthy of much effort in view of the fact that their values 
can be obtained from a good table of logarithms.' At any rate the bibliography of 17 titles? 
(p. xvi-xx) included mostly trivial tables. The well known table of Hayashi* is perhaps the 
only one approaching the present one in extent. The former table, however, is poorly ar- 
ranged for these functions and is unreliable. The present volume, therefore, represents a 
significant advance in the tabulation of hyperbolic functions. 


- =. 
1Sinh™ x = In[x + (x? + 1)4] = csch™ (1/x), —-o <x< ow 
cosh! x = + In[x + (x? — 1)!] = sech™ (1/x), 1i<x<o 
tanh x = $In[{(i + x)/(1 — x)] = coth™ (1/x), —-1<x<+1. 


2 After listing the tables of inverse hyperbolic functions in M. Bott, Tables Numériques 
Universelles. Paris, 1947, p. 475, there is the remark: ‘“This single page of tables is replete 
with errors.” See MTAC, v. 2, p. 336-338; v. 3, p. 466-467. 

3K. Hayasul, Sieben- und mehrstellige Tafeln der Kreis- und Hyperbelfunktionen, Berlin, 
1926, p. 9-201. 


666[F].—D. H. LEumer, ‘‘On the converse of Fermat’s theorem,’’ Amer. 
Math. Mo., v. 43, 1936, p. 347-354; v. 56, 1949, p. 300-309. 


Fermat’s theorem states that if a is any integer, and p a prime, then a? — a is divisible 
by p. For a = 2, the strict converse, attributed to the ancient Chinese, would assert that if 
n divides 2" — 2, then n is a prime. SARRUs in 1819 disproved this by the counter-example 
n = 341 = 11 X 31; 2"! — 2 is divisible by 341. A true converse of Fermat's theorem, 
due to Lucas, states that if m divides a* — 1 for x = m — 1, but not for x a proper divisor 
of m — 1, then n is prime. In both theoretical and practical work it is frequently required to 
settle whether a fairly large number is prime or not. No general usable test is available 
(WiLson’s theorem is unusable). ‘‘One may ask for a criterion which is almost character- 
istic of primes and then tabulate those relatively few composite numbers which also satisfy 
this condition.” The false Chinese converse is such a criterion. But to find the composite 
exceptions, like 341, factorizations of large numbers are frequently demanded. The author's 
long experience in this enabled him in the first paper to publish a list of composite exceptions! 
between 107 and 10® whose least factor exceeds 313. The smallest prime factor of each is 
listed. In the second paper the list is extended? from 108 to 2-108. The Army Ordnance’s 
ENIAC all-electronic digital computer made the work possible in a remarkably short time. 
Incidentally it is proved that there is an unlimited number of n’s dividing 2" — 2 that are 
the products of three primes. 

E. T. BELL 
California Institute of Technology 


1 This list contained 526 entries believed to contain all exceptional numbers; no number 
having more than two prime factors, each in excess of 313. In the second paper it was dis- 
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covered after recomputation that in this list there were 2 entries to delete and 7 to insert 
(MTAC v. 2, p. 279) so that there are exactly 531 exceptional numbers in question less than 
108.—EpiTor. 

2 In the second list are 327 entries—320 two prime-factor numbers and 7 three prime- 
factor numbers. Thus there are in all 858 exceptional numbers <2-10*.—Eprtor. 


667(F].—D. H. Leumer, “On the partition of numbers into squares,” 
Amer. Math. Mo., v. 55, 1948, p. 476-481. 


P;(n) denotes the number of partitions of m into k integral squares = 0. Partitions differ- 
ing only by permutations of the squares are not counted as distinct. The author is concerned 
with the solutions of Pi(m) = 1. It has long been known that P:(n) > 0 for k= 4. The 
author solves completely P,(m) = 1, and finds all numbers that can be partitioned into 4 
squares in precisely 2 ways. For P;(n) = 1 he finds the solutions 4°m, a = 0, for m = 427. 
Any further values of m exceed 10,000. 

E. T. BELL 


668[F].—Nus Pipp1nG, “Diagonalkettenbriiche,” Aabo, Finland, Akademi, 
Acta, Mathem. et Phys., v. 16, no. 5, 1949, 23 p. 


The conclusion of this paper (p. 16-23) presents a table of the expansion of D# in a 
diagonal continued fraction for non-square integers D from 501 to 1000. This is an extension 
of a previous table (MTAC, v. 3, p. 96). In 97 of the 491 expansions the diagonal continued 
fraction is regular. 


D. H. L. 


669[F].—D. YarpDEN & A. Katz, ‘‘Luah sedarot nesiga linariyot binariyot 
misseder 3’’ [Table of binary linear recurring sequences of order 3]. 
Riveon Lematemattka, v. 2, p. 54-55, 1948. 


Four recurring sequences defined by 
Us, = Un_-2 + Un-s Uo = 0, U; aon 0, U2 =1 
Vin Vn-2 + Va-s Vo = a Vi = 0, V2 = 2 
U,= — Ua-2+ Un-3 Uo = 0, U; = 0, U:=1 
YT, =a ot + _ a V, = 3, V; = 0, v; =l ae 2 
are tabulated for +n = 0(1)64. The factorizations of the numbers are also given with the 
exception of Vez = 2-9803919989. Prime values of these four functions are underlined. 


D. H. L. 


670(I|.—H. E. Sauzer, Table of Coefficients for Inverse Interpolation with 
Central Differences. NBS, Mathematical Table, no. MT 27. Second ed., 
June 1, 1949. Washington, D. C., ii, 15 p. 14.9 K 23.2 cm. 


This table, previously reviewed in MTAC, v. 1, p. 315-316, is here reprinted with per- 
mission of the editors of Jn. Math. Phys. No changes were made in the text or tabular ma- 
terial. 


671[K].—Haro_p JEFFREYS, Theory of Probability. Second ed., Oxford, 
Clarendon Press, 1948, viii, 412 p. 15.4 X 23.5 cm. 30 shillings; New 
York branch office price $9.00. Appendix, p. 396-404, ‘Tables of K.” 


The author uses K as a generic term for the ratio of probabilities P(q|@H)/P(q'|@H), 
where q is the initial hypothesis, g’ the alternative, H the previous information, and @ the 
observational evidence. The actual seven tables given are double-entry tables, entered with 
values of n, or v, and of K, to read values of x’, # or z, depending upon the table. The par- 
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ticular function K is separately defined in each table. Each resulting table gives, of course, 
selected values of a specified mathematical function of two arguments, and does not depend 
for its calculation upon any particular theory of probability. The values of K are 1, 10-#, 
10-1, 10-4, 10-, in each table. In T. I-IV, n or »(= nm — 1) ranges from 5 to 100,000 (save in 
T. II, where n starts at 7). The entries proceed by one unit up to 20, by ten units to 100, 
thence the numbers go by successive factors of 2 and 5, up to 100,000. The functional values 
to be read are to one place of decimals, and hence to two or three significant figures. In 
T. IIIA» = 1(1)9, T. IVA, » = 1(1)8, and in T. V, » = 1(1)10(2)20, 50, while z is read to 
2D. In T. IIIA, IVA, and V, the values of K for ¢ = 0, or z = 0, are also listed for the given 
values of ». The argument, K, is related to m or vy and the function to be read, as follows: 


T. 1: K = (2n/x)' exp (—4,°); reading x’. 


T. Il: K = 4anty exp (— 4,2); reading x’. 
T. II: K = (x»/s)§{1 + (2/v) PO; reading #. 








~ ae ad ns??? dv 
ie - ™ Se} { - 
T. WIA: K+ = Q/x) J. sR {40 — 9), 4, - mh ew - a mh oe 
and F(a, 7, x), is the confluent hypergeometric function, 1 + = fe ale + 1) Sie 
y 2by(y7+1) 
t = vit/s’; reading ?. 
T.1V: K = 4eat[1 + (2/v) T-"; reading #. 
»Ki=-2,1f° { = n(a® + B)v? exp { -; nso 
T. IVA: K2 = 287 f-" R41 — 4, 1, - Br eeH ee eT 


d: 
x or where 1F,(a, y, x) is as given above, and a” + 0 replaces 2”; reading ?. 
T. V: K7 = 24-7 9 Llu + 1)/(u* + 1) Ju" exp{int?(1 — u*)}du where b = e*; reading z. 


In general these tables serve to indicate with what degree of accuracy it is worth while 
to carry out the reduction of a given set of observations. Though differing in principle from 
the use of probability integrals, they lead to substantially the same numerical estimates. 


ALBERT A. BENNETT 
Brown University 


672[K].—N. V. Smirnov, “Table for estimating the goodness of fit of em- 
pirical distributions,” Annals Math. Statistics, v. 19, 1948, p. 279-281. 
17.4 X 25.2 cm. Reprinted from N. V. Smirnov, “On the estimation of 
the discrepancy between empirical curves of distribution for two inde- 
pendent samples,’’ Moscow, Univ., Bull. Math., série internationale, v. 2, 
fasc. 2, 1939, p. 15-16. 14.7 X 21.8 cm. 


This is a table of L(z) =1—2 3B (-1)e?"" = z(- 1)’e2"** = (29)tg 
v=] 


x Be 2-1)? 7/82? mostly 6D, for z = .28(.01)2.5(.05)3. The values vary from .000001 for 


for zs = .28, to .999 99997 for z = 3. 

It may be noted that in the Annals, p. 178 and 279 (W. FELLER) the expression for L(z) 
in the first of the three forms given above is incorrect, since it has e~”***. This error is copied 
in Math. Reviews, v. 9, p. 599. The second and third forms, given by Smirnov, are correct. 
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673[K, L].—J. I. Marcum, Tables of Hermite Polynomials and the Derivatives 
of the Error Function. U.S. Air Force. Project RAND no. P-90. Santa 
Monica, Cal., 29 December 1948, 241 leaves, text on one side only. 
Hektographed. 20.5 X 28 cm. Not available for distribution at present. 


Introductory text: The error function may be defined by (x) = E’(x) = (22)~te-*’. 
The n* derivative of this function is then 
o%(x) = D*E'(x) = (20)-td*e-** /dx*. 
The following tables give values of the derivatives of the error function so defined for 
nm = 1(1)10, x = 0(.01)12; 6S. 
In the process of computing these tables, the Hermite polynomials, defined by 
hn(x) = (—1)"e*dne-** /dx*, 
were generated as a by-product. The tables also include the values of these polynomials. 

The specific purpose for computing these tables was for use in GRAM-CHARLIER series 
approximations to certain distribution functions for a sine wave plus random noise [see 
J. I. Marcum, A Statistical Theory of Target Detection by Pulsed Radar—mathematical ap- 
pendix. Project RAND, Report no. R-113, July 1, 1948 (restricted) ]. For use in the Gram- 
Charlier series, m should go at least to 9 and x to at least 10 in interval .01. No combination 
of available tables could fulfill these requirements. Hence these tables were computed on the 
IBM machinery at Project RAND. 

A table of exact values of h,(x) for m = 2(1)6, x = 0(.01)4 was given by N. R. J@rGeEn- 
SEN Unders¢gelser over Frekvensflader og Korrelation. Diss. Copenhagen, 1916, p. 196-205. 
There is also a table of ¢*(x) in T. C. Fry, Probability and iis Engineering Uses. New York, 
1928, p. 456-457, m = 1(1)6, x = [0(.1)4; 5D]; and in New York, W.P.A., ms., m = 1(1)14, 
x = [0(.1)8.4; 20D]. 

Epitortat Note: In the Marcum tables 5D uniformity in display for both h,(x) and 
¢"(x) is attained by supplying extra columns for each function indicating the number of 
places the decimal point is to be moved to the right or to the left. Thus for x = 12 and 


n = 10 the first 6 figures in a 11 figure value for 4,(x) are given, while for ¢*(x) we have a 
value to 27D. 


674[L].— FREDERICK W. Grover, Inductance Calculations. Working Formu- 
las and Tables, New York, Van Nostrand, 1946, xiv, 286 p. $5.75. 
15 X 23 cm. 


Among 55 tables in this volume Dr. FLETCHER has already referred to T. 13, 36, 37, 40, 
49-51 in his “‘Guide to tables of elliptic functions,” MTAC, v, 3, p. 235, 238, 239, 253. 

Auxiliary tables 3-4, p. 238-247 are of zonal harmonic functions and differential coeffi- 
cients of zonal harmonics. T. 3: P,(x), m = 2(1)8, and xP, (x), m = 3(2)7; T. 4: P,’(x), 
nm = 3(2)7, A, A*, or A*; each table giving exact values for x = 0(.01)1. Such tables are also 
to be found in A. H. H. TaLLovist, Grunderna af Teorin for sferiska Funktioner. Helsingfors, 
1905. 


R. C. A. 


675[L].—D. R. Hartree & S. Jonnston, “On a function associated with 
the logarithmic derivative of the gamma function,” Quart. Jn. Mech. 
Appl. Math., v. 1, 1948, p. 29-34. 15.3 X 23.4 cm. 


Here is a table of Re[y(k) — Ink + $k] = Rey(10 + i|k|) — log |k| — || 
1 2 9 
+ --- +——_], where y¥(z) = I'’(z)/T(z), for y = k?* 
Pll +1 * By +1 ly] +1 ood tdcerkasaild 
=[—1(.01) +1; 8D], A*. The values were evaluated to 10D. The maximum error in a 
tabulated value should not be more than .6 in the eighth decimal. Linear interpolation is 
adequate to give 6D accuracy. 





Extracts from text 
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676[L].— VERA HuckKEL, Tables of Hypergeometric Functions for Use in Com- 
pressible-flow Theory, U. S. National Advisory Committee for Aero- 
nautics, Technical Note, no. 1716, Oct. 1948, 13 p. 19.9 K 26.3 cm. 


In the hodograph method of treating plane potential compressible flows the differential 
equation, originally obtained by CHAPLYGIN in his study on gas jets, plays a significant role. 
This paper tabulates various hypergeometic functions which arise as particular solutions of 
Chaplygin’s differential equation. The tables should prove useful in the tabulation of other 
auxiliary functions which may arise in various compressible-flow problems. The adiabatic 
index for air has been taken as 1.4. 

Chaplygin’s differential equation is 


(1) or(1 — r)d*¥z/de* + ((R +1) — (8 +1 — B)r]dYi/dr + 36k(R +1)V% = 0 


where 8 = 2.5, M? = 5r(1 — 1), M being the Mach number, and Y;(r) ~ 1, as r > 0. 

The tables of numerical values have been prepared for a selected number of the com- 
plete set of solutions of (1). These solutions extend the results of Chaplygin into the super- 
sonic range and to negative values of the index k. 


Tables 1-2 are for the functions ¥;, k = — 15(.5) + 15, 7 = [0(.01).5; 5D], with 
corresponding Mach numbers. 
Tables 3—4 are for the functions d Y¥;/dr, k = — 15(.5) + 15, r = [0(.01).5; 4D], with 


corresponding Mach numbers. 
REFERENCES 
I. E. Garrick & Cart Kaplan, On the Flow of a Compressible Fluid by the Hodograph 
Method II—Fundamental Set of Particular Flow Solutions of the Chaplygin Differen- 
tial Equation. NACA, Report, no. 790, 1944, 21 p. 
S. A. CHapitycin, Gas Jets, NACA, Technical Memorandum, no. 1063, 1944, 112 p. +3 
leaves. Translated by S. Reiss from the Russian in Moscow Univ., Uchenyia Zapiski, 
Otd. Fiz.-Mat., v. 21, 1904, p. 1-121. 
Extracts from text 
EpitoriAL Note: Miss HUCKEL has acknowledged to us that her statement regarding 
what is tabulated in Tables 3-4 is incorrect, that what she really has tabulated is 


—(2/B8k)d Y;/dr, as suggested by Dr. J. C. P. MILLER in Math. Reviews, v. 10, May 1949, 
p. 329. 

In the Report of Garrick & KapLan there are a number of Tables including a 5D 
table of Y; for k from —5 to +5 and with M as parameter varying from .1 to 2, and with 
corresponding values of r. 


677(L].—A. N. Lowan & WILLIAM HORENSTEIN, On the Function H(m, a, x) 
= exp (—ix)F(m + 1 — ia, 2m + 2; ix). NBS, Mathematical Table, no. 
MT 19. Second ed., June 15, 1949. Washington, D. C., ii, 20 p. 14.9 
X 23.2 cm. 


This table, previously reviewed in MTAC, v. 1, p. 156, is here reprinted, with permis- 
sion of the editors of Jn. Math. Phys. No changes were made in the text or tabular material. 


678[L].—WILHELM MacGnus & FRITZ OBERHETTINGER, Formulas and 
Theorems for the Special Functions of Mathematical Physics. Translated 
from the German by JOHN WERMER. Published and distributed in the Public 
Interest by authority of the Attorney General under License No. A-1341. 
New York, Chelsea Publ. Co., 231 West 29th St., 1949, viii, 172 p. 
$3.50; 15% discount to teachers and members of the Amer. Math. Soc. 
on direct purchase from the publishers. Offset print. 
The first German edition of this valuable work (1943), from which this English transla- 

tion was made, was reviewed in MTAC, v. 3, p. 103-105. Since we there gave the detailed 
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table of contents it will not be necessary to repeat them here. Most of the errata which we 
then listed have been corrected in this translation, but those noted below still remain. The 
translation and original correspond page for page throughout. 

But we reviewed also the second enlarged and greatly improved German edition (1948), 
in MTAC, v. 3, p. 368-369, which may be procured comparatively easily. Hence many 
people who might have purchased the English edition may now prefer the better German 
work. 

The errata referred to above are the following: 


P. 24, 1. 13, for u3, read w*. 

P. 46, |. 13, before =0, insert y. 

P. 73, exponents yu, of lines 3 and 5, are badly placed. 

P. 83,1. 11, for ay are realanda +1 > y > 1, read wand yare real, and y > 0,a > y — 1; 
1. 14, delete. 

P. 88, 1. —1, for 0 + 1, read 0, +1. 

P. 108, 1. 3, for x[3 times], read ft. 

P. 119, heading, for §2, read §1; col. 2, i. 5, for y|, read |y|. 

P. 131, 1. —1, for the image function of the original function f’J,(at) in the Laplace trans- 
form table, the second edition has 


Tutr+i 

pers oe MF te ti), Mete+2); wt+1; —a*/p*), 
1 

not (—1n ATED p niyo + at)Re(v + m) > — 1. 


(P? + a?)hot) 
me. A, 


679([L].—NBSCL, Tables of ithe Confluent Hypergeometric Function F(}n, 4; x) 
and Related Functions. (NBS, Applied Mathematics Series, no. 3.) 
Washington, D. C., 1949, xxii, 73 p. 19.8 X 26.1 cm. For sale by the 
Superintendent of Documents, Washington, D. C., 35 cents. See MTAC, 
v. 3, p. 414, 483. 


The confluent hypergeometric function F(a, 7; x) is defined by 
(1) F(a, ¥;*) = 2 P(y)(@ + j)xi/(T(@)l(y + A)7!1. 
i- 


It is easily verified that it is the solution of the differential equation 

(2) xy” + (y —x)y’ — ay =0 

satisfying the initial conditions y(0) = 1 and y’(0) = a/y. KumMer’s formula 

(3) F(a, y;x) = eF(y — a, y; —x) 

is a direct consequence of this property and may be verified by substitution into (2). In this 
review we use the abbreviation F, = F(4n, $;x). Using the series expansion (1) for the 


right side in (3) it is readily seen that when m is a positive odd integer, then F, equals e 
multiplied by a polynomial of degree }(m — 1). In particular 


(4) F, = é, F; = e*(1 + 2x), F; = e*(1 + 4x + $x°). 
The ratios g, = F,/Fn-2 satisfy the recurrence relation 
(5) En42 = (2x + 2m — 1)/n — (m — 1)/ngn. 


Formulae (4) and (5) permit to calculate successively all F,(x) for odd n, and. this method 
was used in the preparation of the Tables. Various other relations were used to check the 
accuracy of the tables, and it is believed that the error is less than one unit of the last decimal 
place. 

In the present tables the subscript » of F, is always an odd integer. The original objec- 
tive was the tabulation of In F,, but it turned out that the function (2mx)-+ In F, is more 
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convenient for interpolation and has the added advantage of being of the order of magnitude 
of unity, over a considerable range for m and x. The main part of the tables (p. 14-73) is 
devoted to the tabulation of this function for m = 3(2)201 and x = .1(.01).6(.1)2(.2)7(1)- 
45(5)100. For m = 3(2)21, the range .6 to 2 is covered in steps of .05. Rows correspond to 
values of x, columns to ”. On p. 2-6 the function F, itself is tabulated for x = 0(.01).1 and 
n = 3(2)201, and on p. 8-11 we find F,/cosh{(2m — 1)x}# for the same range of x and 
n = 43(2)201. The latter table is given for interpolation purposes since in the range covered 
F, varies very rapidly. All tables are to 6D, except the table of F, on p. 3-6 which is to 7S. 

A pleasing feature of the volume is a comparatively large collection of formulae and 
asymptotic expansions for F, and certain related functions. (Most of these formulae are not 
directly connected with the tables.) An interpolation chart shows how many terms of the 
bivariate Lagrange interpolation formula are required for various and x in order to obtain 
an accuracy to 6D. 

Confluent hypergeometric functions appear in many applications. The present tables 
were designed especially for statisticians. Their primary purpose is in connection with the 
sequential é-test designed by A. WALD. It is also of great use in the analysis-of-variance and 
the ordinary /test. 

W. FELLER 
Cornell University 


680[L, M, S].—C. J. BouwKkamp, “‘On the effective length of a linear trans- 
mitting antenna,” Philips Research Reports, Eindhoven, Holland, v. 4, 
June 1949, p. 179-188. 15.4 X 23.4 cm. 

F,(H) =So¥ Ji(x) cos xdx/x, Fo(H) =So#Je(x) sin xdx/x, C(4H)=So#(1—cos x)dx/x, 
Hes = 2F\(H) sin H, Hee = 2F:(H) cosH, R = 30[C(4H) — sin* 2H], R, ~ 80|H,,.|*, 
Re = 64|Heo|?. 

In Table I are 6D values of Fi(H), F2(H), C(4H) for H = 0(.1)3. 

In Table II are given 4D values of He,., He,a, Rs, Ra, Re + Ro, R, for H = .1(.1)3. 


Extracts from text 


681[L, M].—Haro_p LEvINE & JULIAN SCHWINGER, “On the theory of 
diffraction by an aperture in an infinite plane screen. I,’’ Phys. Rev., 
s. 2, v. 74, 1948, p. 958-974. 19.8 X 26.6 cm. 
On p. 972-974, So'u-"-™ (1 — v*)4Tn4(av)Jm4;(av)dv is evaluated for m = 1, 2, and 
a = 1, 2. 


682[L, S].—J. MEIXNER, Lamé’s Wave Functions of the Ellipsoid of Revolu- 
tion. Translated by Mary L. MAHLER from ‘‘Die Laméschen Wellen- 
funktionen des Drehellipsoids,” Zentrale fiir wissenschaftliches Berichts- 
wesen der Luftfahrtforschung des Generalluftzeugmeisters (ZWB), 
Berlin-Adlershof, Forschungsbericht no. 1952, June, 1944, for the Na- 
tional Advisory Committee for Aeronautics. Technical Memorandum, 
no. 1224, Washington, D. C., April 1949, iv, 102 p. 19.8 26.3 cm. 
The basic differential equation treated is 

(t) el a-we] +( 8 - Herta) sr 

dn dn 1-7 ; 

Corresponding to given parameters » and y = kc, there exists a countably infinite set of 

eigenvalues \,“(y) corresponding to which (1) has a solution, called the eigenfunction, which 

can be represented by 
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(2) XPM(E; y) = = arab. (vy) Peer (€), 
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where P%,,(£) is the associated LEGENDRE function of the first kind and r is always even. 
It should be noted that if in (1) we set f = (1 — 7*)*w, we obtain 


(3) (1 — 9*)w” — 2(u + 1)qw’ + C— wu + 1) — Bc*g? + Aw = 0. 


If » is an integer m, (3) yields the prolate spheroidal wave functions (oblate if y = iC, 
C real). Fairly extensive tabulations of both the eigenvalues and coefficients proportional 
to an, for integral m, are available.* In fact, the author refers to this text, and his mathe- 
matical development closely follows it—at least so far as the determination of the eigen- 
values, and the derivation of the relations existing between the various solutions of (2). 
The author notes that if u in (1) is equal to —}, we obtain MATHIEU’s equation (this was 
already noted.*) Although the mathematical treatment applies to general values of yu, the 
tabular material is restricted to integral values of wu, which are designated by m, as in 3. 
Before describing the tables, it will be necessary to define a second solution of (1): 


XR OME: y) = Daal o(y)Paer(®) + Zsa (veer (7). 
In the above Onyr (é) is the associated Legendre function of the second kind, and r’ is some 
finite integer. 
It may be shown that 


Xn (EE; y) = Zier) Patel), 
where 
bn. r(y) = Ce + m +1) (in — m)!/(n +7 — m) Mn + m)! Jan r(y) = ant (y)- 
The tables comprise the following: 


Table 1. Coefficients of +’, y‘, y°, vy’, and y"® in the power-series expansion of A,™(), m=0(1)9, 
n = m(1)9; 10D. 

Tables 2, 3, 4, and 5. Coefficients of y*, y* and 7° in the power-series expansion of ay: + /an,o 
of (2), m = 0(1)9, nm = m(1)9, +r=[2, 4, 6; 10D]. (The coefficient of y** is zero if 
k < |r|). 

Table 6. The coefficients of y*, y* and y® in the expansion of a; ,/ano; m = 0(1)2, 
n = m(1)(5 — m);r = [— 6(2) — 2; 10D]. 

Tables 7, 8, 9, 10. The coefficients of y*, vy‘, and * in the power-series expansion of bf ,/Dh 0; 
m = 1(1)9, m = m(1)9, +r = [2(2)6; 10D]. (For m = 0, bn, + = an, +). 


The author points out that these truncated power series are useful only for relatively 
smal! values of , when 2 is of low order, but that more accuracy is obtainable for high values 
of n. 

In addition to the above tabular material, the author inserted some of BouwKAMP’s 
tables of eigenvalues! for m = 0, and values of the corresponding eigenfunction or its deriva- 
tive at § = 0 and 1. There are seven tables of this character, comprising seven pages. Cor- 
rections which Bouwkamp made after publication were apparently not known to Meixner; 
as a result the last two places of these tables are totally unreliable, and more serious errors 
also remain uncorrected. It should be noted that nearly all values of 1 which Meixner used® 
were corrected and amplified? (see review‘). Since 2 is relatively accessible, Meixner’s re- 
production of Bouwkamp’s erroneous values should be discarded; it is unfortunate that the 
translator did not replace these seven pages with the more useful tables.” 

Meixner’s new results are contained in Tables 1-10; they duplicate some results of 1 
only for m = 0. Since Meixner gives no tabular values for the eigenvalues or coefficients, his 
results do not duplicate the tables in 3 but rather add to them. Some asymptotic formulae 
for the eigenvalues and functions are also given, and there is a graph of A,»(y) for m = 0, 
n = 0(1)9. The text contains a valuable compilation of theoretical results, with some dis- 
cussion of their application to various problems. 
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The NACA performed a genuinely useful service in making this work available in 
English. 


GERTRUDE BLANCH 
NBSINA 


1C. J. Bouwxamp, Theoretische en numerieke Behandeling van de Buiging door een ronde 

Coonip- Diss. Groningen, 1941 
Bouwkamp, “On spheroidal wave functions of order zero,” Jn. Math. Physics, 

v. 26, Gol: p. 79-92. 

J. A. "STRATTON, P. M. Morse, L. J. Cau, & R. A. Hutner. Elliptic Cylinder and 
Spheroidal Wave Functions. New York, 1941. See MTAC, v. 1, p. 157-160. 

4 MTAC, v. 3,.1948, p. 99-101; review of items 1 and 2. 
‘ be oe corresponding to y? = 15, 20, 25, 50, and 100 are the only ones not repro- 

uced in 2. 


683[M].—Erik HALLéEN, Iterated Sine and Cosine Integrals, Roy. Inst. 
Techn., Stockholm, Trans., no. 12, 1947, 6 p. 17.6 X 24.7 cm. 


In antenna theory there are some integral functions which might also be of interest 
elsewhere. They derive from sine and cosine integrals in the same manner as these derive 
from sine and cosine. The following are tabulated: 


Liu(x) = Cu(x) + iSu(x) = Set7L@dt, 
Lor(x) C(x) ob 4So1(x) = Sete L (28) - L(t) jdt, 


x = [0(.2)14;5D], where L(x) = C(x) + iS(x) = Set-(1 — e~**)dt, C(x) = y+Inx 
— Ci(x) and S(x) = Si(x). There is also a table of Li:(x), for x = [14.2(.2)28; 5D]. The 
tables were computed by mechanical integration. The fifth decimal is uncertain within 
two units. 


Extracts from text 


EpitoriAL Note: For a less extended 6D table of Lui, by C. J. Bouwkamp, see MTAC, 
v. 3, p. 303-304 


684[M].—J. Le Caine, A Table of Integrals involving the functions E,(x) 
= fi*e-™“u-"du. National Res. Council, Div. of Atomic Energy, Report 
no. MT 131. 1948, i, 45 leaves. 20.3 X 27.4 cm. 


The definite and indefinite integrals involve products of two or more factors E,(kx), 
E,(ax + 6), En(k|x — y|), e*, e*!-ulxm, In the single integrations x is the integration 
variable, and x and y in the double integrations. For numerical tables of such integrals see 
MTAC, v. 2, p. 260, 272 and v. 3, p. 303-304, and RMT 683. 


 <. & 


685[M, Q].—H. C. Van DE Hutst, “Scattering in a planetary atrosphere,”’ 
Astrophys. Jn., v. 107, 1948, p. 220-246. 16.9 K 24 cm. 


In this paper the problem of multiple scattering in a planetary atmosphere, both with 
and without a diffusely reflecting bottom surface, is discussed. We assume that the atmo- 
spheric scattering is isotropic with an albedo a, and that the ground surface reflects the 
radiation according to Lambert’s law, with an albedo 6. The mathematical functions occur- 
ring in the calculations are of three types. T. 1: for E,(x) = /i*t-*e~**dt, E,® (x) = St Edt 
the functions tabulated are E;(x), E(x), Ey (x) = e*, Es(x), and also In x + vy (vy being 
Euler’s constant) are each tabulated to 6D; and 2E;(x) to 5D. T. 2-3: For F,(d, x) 
=JSe"E,(t)dt, there are 4D tables for F,(—s,x), s = 10, 5, 2, 1; and of e~*F,(s, x), 
s = 0, 1, 2, 5, 10, in each case m = 1,2. T. 4: For Gmn(x) = Se*Em(t)En(t)dt, G’mn(x) 
= S¢Enm(t)En(x — t)dt, there are 5D values for m, n = 1, 2; also there are 4D values for 
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g(x) = 2 — 2E2(x) — Gio(x) — Gis’(x) and c(x) = [x — 4 + E,(x)]/x. In each of the 
tables x = .01, .02, .05, .1(.1).4(.2)1(.5)3, «. 
Extracts from text 


EprtoriaL Note: For considerably more extensive tables of E,(x) see MTAC, v. 2, 
p. 272. 


686[N].—FINANCIAL PUBLISHING Company, A. Eighth Coupon Rate Monthly 
Bond Values showing accurate monthly values on bonds and other re- 
deemable securities paying interest semi-annually. For eighth coupon rates 
only from 13% to 34%. Boston, Financial Publishing Co.; London, 
Routledge & Kegan Paul Ltd., 1949 [the printing and binding were not 
completed until Feb. 1949], x, 537 p. 25 XK 36.8 cm. Bound in red imita- 
tion leather, thumb index. $125.00. In the back of the volume is a loose, 
slightly smaller (24.5 X 36 cm.) paper covered reprint, 48 p., dated 
1948, from the second edition (p. 1109-1154, 0% coupon) of Monthly 
Bond Values, 1941, and entitled: Tables of Present Worth of 100, from 1 
month to 110 years. 
B. Executive Bond Values Tables. Desk Edition, Showing net returns on 
Bonds and other redeemable Securities paying Interest semi-annually. 
Boston, Financial Publishing Co.; London, George Routledge & Sons, 
Ltd., 1947, 1984 p. 11 X 18.6 cm. Bound in leather, full gilt, thumb 
index. $25.00. 


A. In the 8 years since the second edition of Monthly Bond Values was published (see 
MTAC, v. 1, p. 114-115), continuing low interest rates and the development of sale at 
competitive bidding of corporate bonds have produced a number of issues with }% coupons 
which were not covered by complete tables in the second edition. Hence the need for this 
new volume which is essentially supplemental to the second edition. 

The 12 coupon rates here are 14%(4%)34%. The yields are from 0 to 5.70%, at interval 
.05%, as before. Maturities, as in the earlier volume, progress by monthly intervals up to 
40 years 3 months, from 403 years to 50 years semi-annually, 50 years to 65 years annually, 
and 70 years to 110 years by 5 years. All tabulation is to 6D, adjusted from at least 8D 
computation. 

As formerly, on every page are 12 columns (one for each month) and 115 lines corre- 
sponding to 0 (.05%) 5.70%, the various yields (in black face type on each side of the page), 
making 1380 entries all arranged in groups of five, admirably spaced for ready use. Also 
in black face type are given under the appropriate columns the amounts of accrued interest 
on $100. for 1, 2, 3, 4, 5 months. Pages 536-537 contain a 6D table of accrued interest on 
$100. for any number of days 1(1)180, in 6 months, for each of the coupon rates. An 8D 
“Table of factors” is provided (p. 532-534) for computing values for 1 to 30 days, between 
the even months, for each of the yields. 

In the foreword to the volume are explained various methods of using and applying the 
tables which are largely a reprint of the second edition. Part of what was said in review of 
the earlier volume is applicable to the present book, and the comments made at that time 


. apply also to this volume. 


The manner of preparing the tables was almost wholly dictated by trade practice and 
by what would be most intelligible to bond dealers. The method of computation involved 
the addition of the constant difference between values at successive coupon rates. This 
constant difference was computed to a minimum of two decimal places beyond the six 
decimal places contained in the table, and in most cases to more than two extra decimals. 
The mathematician notes with interest that trade practice calls for the addition to every 
invoice of simple interest at the coupon rate. Hence it became necessary to compute the 
theoretically correct ‘‘flat” value and from this to subtract simple accrued interest. Thus 
when an invoice is made for bonds sold at a yield basis the sum which changes hands will 
be the theoretically correct amount, but it will be itemized in two distinct elements neither 
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of which may be considered mathematically exact. This feature is most readily observed in 
the case of any bond selling to yield its coupon rate. The value is shown as 100 on coupon 
dates only, and for intermediate dates at something less than 100; but if the accrued interest 
at the bottom of the column be added a correct result will be obtained. 

Another curious feature is that in spite of solicitude for correct semi-annual compound 
interest in the body of the book, values for less than six months’ maturity are shown at true 
discount, in accord with trade custom of dealers and banks. Again, the editors believed that 
labelling a present-worth table as the value of a 0% bond evokes a clearer concept to the 
average bond dealer, than a proper labelling would. 

Without a volume of this character, dealing in municipal securities would be a some- 
what difficult process. For example, a dealer offers a 3% coupon municipal having 20 years 
and 4 months to run on a 1.95% yield basis. What dollar price must the purchaser pay? 
Turning to the 33% coupon table, find the page “20 years, column 4 months, and on the 
margin the yield 1.95%,” and the price is given as 125.913848 or $1259.14 per $1000 bond. 
Also at the bottom of the column is given the accrued interest in monthly intervals, $.583333 
per $100 or $5.83 per $1000 bond. The methods of adjusting for fractional parts of a month, 
and for many other special cases of practice, are explained clearly in the introductory pages. 
These pages contain also quite explicit instructions involving the use of adding or account- 
ing machines where these are available. 

These 1941 and 1949 tables are the most extensive of their kind, and are of value par- 
ticularly to dealers and those handling state and municipal securities where price is usually 
quoted on a yield basis, also for insurance companies where exact computations of yield are 
essential. The tables are obviously of importance also to the actuary or accountant in mak- 
ing up amortization schedules. 

B. Since the 1941 and 1949 volumes are too large and bulky for use outside an office, 
two small volumes have earlier found considerable favor. One of these, Comprehensive Bond 
Values, 1936, had the advantage of containing tables for odd eighth coupons, now being 
increasingly used. The other, Pocket Edition of Monthly Bond Values, 1939, had the advan- 
tage of monthly maturity intervals but with no eighth coupons. The volume under review 
contains a combination of the material in these volumes. 

It covers coupon rates $%(4%)4%(4%)5%. The maturities for }% coupons are monthly 
to 20 years, quarterly to 30 years, semi-annually to 50 years, annually to 56 years, bi- 
annually 56 years to 60 years, and at 5 year intervals from 60 years to 100 years. On the 4% 
coupons maturities are shown monthly to 10 years, semi-annually to 42 years, annually to 
56 years, bi-annually to 60 years, and at 5 year intervals to 100 years. Coupon rates for 4% 
and 2% stop at 28 years. For §%, 4%, 14% and 13% they stop at 30 years. The tables are 
in condensed forms and show values to 2D only. 

In the foreword to the book are descriptions of methods of use of the tables and at the 
end are tables of factors to be used on callable bonds and tables of current yields on bonds or 
stocks, given the price and coupon or dividend rate. There are also included tables for ac- 
crued interest on a 365 day basis from 4% (4%) 4% (4%) 6% for use with government se- 
curities. The volume is very valuable for the convenience of salesmen and others in de- 
termining quick approximate figures. 

F. J. Hunt 
R. I. Hospital Trust Co. 
Providence, R. I. 


687[R].—ItaLy, Istrruto GEoGRAFICO MILITARE, Florence, A. [GIOVANNI 
Boaca ], Sulla Rappresentazione Conforme di Gauss, 1941, ii, 65 p. 
B. [Giovanni BoaGa], Sulla Rappresentazione Conforme di Gauss, 
Allegati, (i) Parte I, Tavole Numerici, 1942, vii, 25 p.; (ii) Parte II, 
Tavole Ausiliarie . . ., 1942, xiv, 49 + Errata Corrige 4 p.; (iii) Parte III, 
A pplicazioni, 1942, 38 p. 15.5 X 21.6 cm. 
This work treats the conformal representation of the earth considered as an ellipsoid of 

revolution onto the Euclidean plane. The tabular material is for that part of the representa- 
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tion which lies between latitudes 36° and 48°, the range in which Italy lies. The main treatise 
gives a theoretical development of the representation with special attention to the prepara- 
tion of formulae well adapted to calculation. 

Let \ be the longitude in seconds of a given point P on the ellipsoid, measured from a 
convenient meridian, and let B be the length of arc on the meridian from the equator to a 
reference point. Next let x, y be the rectangular coordinates of the point P’ in the plane 
corresponding to the point P. Then one has approximately 

log (x — B) = 2 log A + To + Crr* 
log y = logy + T1 + Cy" 
where the quantities To, 71, C2, C; depend on the latitude of the point P. This dependence is 
tabulated in Part I. More precisely the quantities B, T) and 7; and their differences per 
second together with C2, C;, and their logarithms are given for every minute of latitude 
between 36° and 48°. This table is based on the Hayrorp ellipsoid. B is given to a milli- 
meter, the 7’s to 7D and the C’s to 4S or 5S. 

Part II is based on the BEssEL ellipsoid and gives slightly different values for the above 
quantities. Both Parts I, II give three other coefficients used to compute the “convergence 
of the meridian,” the “modulus of linear deformation” and “angular reduction.” These are 
also tabulated for each minute of latitude. The use of the tables is illustrated in Part III. 

L. H. Swinrorp 
Univ. of California 
Berkeley, Calif. 


688[U].—GerorGeE G. HoEHNE, Practical Celestial Air Navigation Tables. 
Volume II, Latitudes 20° to 39° North inclusive. Miami, The Navigation 
Publishing Co., 1943, 256 p. 15 X 24.3 cm. $3.65. 


This volume was reviewed earlier in MTAC v. 2, p. 129-131; the copyright difficulties 
which formerly prevented its sale to the public have been taken care of and it is now gen- 
erally available at the price indicated above. In this connection, it may be pointed out that 
many of the altitudes in the volume were taken directly from the British ANT (or our H. O. 
218 which was a direct copy of it); during the war, both of these sets of tables were re- 
stricted. With the closing of the war, these restrictions were lifted and the way was open for 
arrangements to meet British Crown copyright requirements. 

It is hoped that another volume, similar to this but covering latitudes 40°-59° North, 
will soon be available. 


CHARLES H. SMILEY 
Brown University 


689[U].—O. A. DE AZEREDO RopriIGuEs, Grdficos Usados na Navegagio. 
D. N. 7-1. Rio de Janeiro, Diretoria de Hidrografia e Navegacao, 1948, 
48 p. 16.5 X 24.5 cm. 35 cruzeiros. 


This volume contains twenty nomograms useful in the ordinary computations of naviga- 
tion. The author is a Capitao-de-Corveta in the Brazilian Navy; he already has one other 
navigation table to his credit, Tébuas para Retas de Altura (RMT 516, MTAC, v. 3, p. 114- 
115). The nomograms are well designed for practical use and are clearly printed on good 
heavy white paper. The volume is bound in light green cardboard and presents a neat ap- 
pearance. 

Following a brief preface and a list of abbreviations, there is an index. The remainder 
of the volume is given over to the nomograms and to explanations of the use of a few of the 
more complicated ones. For each one except the last, one or more examples are given with 
the correct answers. 

The first nomogram provides the difference of latitude, dL, and departure, P, for a 
rhumb line, given the distance, D = 0(1')60’(2’)150’(5’)220’, and the course angle, 
C = (0-90°). The corresponding equations are dL = Dcos C, P = Dsin C. 
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The second converts departure, P = 0(1’)100’ into difference of longitude, dLo = 0 
— 220’, given the mean latitude, Ln = 0 — 90°. dLo = P sec Lm. 

The third provides the distance, D = 0 — 100 nautical miles to an object of known 
elevation, E = 0 — 2400 meters, given the vertical angle a = 0(5’)6° and the height of eye 
of the observer, 1-45 meters. A small table converts height of eye into dip of the horizon, 
B = 2’(1')12’. The basic equation is D = RE cot (a — 8) where k is the fraction of a nautical 
mile equal to one meter. The fourth yields the distance, D = 100 — 3500 meters, given the 
elevation E = 0(1)52 meters of an object and the vertical angle a = 0°30’(5’)3°30’ sub- 
tended by it. D = Ecot a. 

The fifth gives the distance, D = 3 — 1000 meters travelled in S = 1 — 60 seconds at a 
velocity, v = 4(0.5)36 knots. D = kvs where k = number of meters per second in one knot. 

The sixth yields the distance D = 12 to 880 nautical miles travelled in H = 25(10™)- 
6»(20™)24» at a speed V = 6(0.5)36 knots. D = VH. 

Number 7 provides the correction y = 0(0°.1)13°.5 to be applied to a radio bearing be- 
fore plotting it on a Mercator map, having given the difference in longitude dLo = 0(0°.5)31° 
and the mean latitude L, = 0(0°.5)65°. The corresponding equation is y = $dLo sin Lm. 

The eighth nomogram is somewhat larger than the remaining charts; it requires two 
facing pages, which are in reality a single sheet folded in the middle. Even though it has been 
placed in the middle of the book and special care has been taken in binding it, the crease will 
probably cause some trouble until the volume is limp from use. This nomogram allows one 
to find the azimuth angle A = 16° to 164° of a celestial body when the local hour angle 
t = 5° to 170°, the declination d = 75° same name —25° opposite name to latitude and the 
latitude L = 0 to 70° are given. The equation solved by the nomogram is 


cot A sec L = tandcsct — tan Ll cott. 


The PAGEL coefficient, cot A sec L, may thus be evaluated with this nomogram. 

Number 9 gives the azimuth angle A = 0 — 180° as a function of local hour angle, 
t = 0 to 180°, declination, d = 0(1°)90°, and altitude kh = 0(1°)90°. The equation here is 
csc A = csc? sec d cos a. 

Nomograms 10 and 11 provide respectively the altitude, h = 0(1°)90° and local hour 
angle ¢ = 0 to 90° of a celestial body when it crosses the prime vertical. The arguments are 
in both cases latitude and declination; in the altitude nomogram, L = 1° to 90°, d = 1° to 
90° and in the local hour angle chart, L = 0(1°)90° and d = 0(1°)60°. The equations are 
csc h = csc d sin L and sect = cot dtan L. 

Nomograms 12 and 13 give respectively the altitude, h = 0(1°)90° and the local hour 
angle ¢ = 10°(1°)90° of a body at maximum azimuth angle. The arguments in the former 
nomogram are L = 1°(1°)90° and d = 1°(1°)90°; in the latter, they are L = 0(1°)60° and 
d = 0(1°)90°. The equations are csc h = csc L sind and sect = tand cot L. 

Nomograms 14 and 15 yield the azimuth angle A = 0 — 180° and local hour angle 
t = 5° to 175° of a body which is rising (or setting). The arguments are respectively L = 0 
to 89°,d = 1°to 90° and L = 0(1°)90°, d = 0(1°)60°. The equations are sec A = csc d cos L 
and sec # = cot dcot L. 

Nomogram 16 provides the correction c = 0(1’)80’ to be applied to a circum-meridian 
altitude, given the local hour angle ¢ = 0(15’)17°30’ and the declination d = 60° same name 
to 30° opposite name h, and latitude ZL with the same range as d. 

Nomogram 17 gives the change in altitude Ah per minute of arc change in local hour 
angle At. With L = 0 to 90° and azimuth angle 0 to 90°, Ah/At = 0(0’.05)1’. The equation 
used is Ah/Ai = sin A cos L. 

Nomogram 18 yields the change in altitude Ah = 0(1’)30’ corresponding to a change in 
latitude AL = 0(1’)30’ and an azimuth angle of 0 to 180°. Ak = AL cos A. 

Nomogram 19 is a simple nomogram for multiplication and division. Chart no. 20 
is intended for plotting lines of position and obtaining fixes. On it the scales of distance and 
latitude vary with the mean latitude, while the scale of longitude is constant. This is in 
contrast to the usual charts of this type in which the scale of longitude varies and the 
other two scales are equal and constant; there is much to recommend this approach. 
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It should be pointed out that where the smallest division of a scale in the nomogram is 
uniform, the interval of the corresponding argument is indicated; otherwise the range of 
values covered by the scale is noted. It will be clear that the accuracy of a given nomogram 
may often be increased by interpolating between the marks provided. 


CHARLES H. SMILEY 


690[U].—D. H. SaDier, “Tables for astronomical polar navigation,”’ Insti- 
tute of Navigation, Jn., v. 2, 1949, p. 9-24. 15.5 K 24.9 cm. 


Summary: The special features of polar navigation are examined with a view to the 
design of tables for astronomical navigation in polar regions. The use of long intercepts and 
curved position lines is thoroughly investigated. It is shown that considerable economy of 
presentation can be achieved ; as an example, a one-page table covering all bodies and both 
polar caps is given. Making the fullest use of polar astronomy, samples are given of perma- 
nent tables for the sun and stars independent of the Air Almanac. 


MATHEMATICAL TABLES—ERRATA 


References to Errata have been made in RMT 659 (Dale), 662 (Briggs), 
665 (Boll), 666 (Lehmer), 672 (Feller), 676 (Huckel), 678 (Magnus & 
Oberhettinger), 682 (Bouwkamp, Meixner), 683 (Hallén); N 106 (Watson), 
107 (FMR, Gauss), 109 (DeMorgan, Rheticus). 


160.—D. H. LEHMER, Guide to Tables in the Theory of Numbers, 1941. 


. 14, table, col. 3, the second and third entries from the bottom should be interchanged. 
. 15, 1. 9, for Creak 1, read Creak 2. 
47,1. —9, for 195 [twice], read 295. 
74, |. 18, for 29 = p, read 17 = p. 
76, 1. 7, for Vorono* read Voronol. 
. 95, 1. 17, for Haupt Exponents, read Haupt-Exponents. 
. 102, Glaisher 16, 1. 2, for 127, read 125. 
. 109, 1. 10-11, for milia accuratis, read millia accuratius. 
. 114, in reference Ostrogradsky 1 for v. 1, 1838, read v. 3, 1836. 

P. 129, the error for m = 1019681, occurs in Burckhardt 1, mot 2. The error for 
n = 2012603 occurs in Burckhardt 2, not 3. Delete 1. 4 (446021) and 5 (446023) since these 
errata were listed by Burckhardt himself, v. 1, p. IV. 

P. 134, |. —6, for 61-330413, read 1399-14407. ¢ 

P. 134, Cunningham 28, p. 243, col. 11, 1. 18, read 9683, instead of 4683; p. 244, in tables 
giving elements y, of primes (y* + 1)/13, omit the entry 671, and for 3930 in col. 4, 1. 6, 
read 2930. 

P. 139, Euler 23, 193 factor of 812 + 1 mot of 82? + 1; 1068* + 1 = 5*-73, mot 55-73; 
773 factor of 1090? + 1, not of 1080? + 1. 

P. 156, Kraitchik 4[f1] for p = 116537, read 116337. 

P. 174-6, the following items in the index are slightly out of alphabetical order: Forms, 
quartic; Ince; Logarithms; Pierce; Sum of powers. 

A few careful readers have pointed out the carelessness of the first definition of the 
function ¥(m) on p. 7. This function has only half the value attributed to it by this definition 
whenever 1 is divisible by 8. 


a a-ha: Ba -Ba- Ba Ba- a Be-) 


D. H. L. 


161.—D. N. Leumer, “‘A complete census of 4 X 4 magic squares,’’ Amer. 
Math. Soc., Bull., v. 39, 1933, p. 764-767. 


This paper contains a list, p. 767, of the number of normalized squares with various 
first columns and last rows. 
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For (15, 2, 1) (10, 5, 3)12, (9, 6, 3)11, (8, 7, 3)2, (8, 6, 4)s, read (15, 2, 1) (10, 5, 3)15, 
(9, 6, 3)1s, (8, 7, 3)s, (8, 6, 4)4; for (14, 3, 1) (8, 6, 4)2, read (14, 3, 1) (8, 6, 4)s; for (12, 4, 2) 
(11, 6, 1)4, read (12, 4, 2) (11, 6, 1)s. These corrections change the total number of 4 X 4 
magic squares from 539136 to 549504. 
S. Urata, A. Sakat & A. HIRAYAMA 
Yamagata-Ken 
Japan 


UNPUBLISHED MATHEMATICAL TABLES 


82[F].—L. Po.tetti, Factor Table and List of Primes for the 30000 natural 
numbers nearest 15,000,000. Manuscript table deposited in the library 
of the American Math. Soc. New York. 


This table gives new information for the range 14984970—-15000000. The second half 
from 15000000 to 15015000, is also covered by W. P. Durfee’s factor table for the 16th mil- 
lion, a table which is in the same library. 

The factor table, which the author calls ‘‘Neocribrum,” is a “‘type 3 table” arranged in 
the usual way modulo 30. On p. 1 are given data on the distribution of the primes in this 
range. Thus there are 1809 primes which are also classified modulo 30. There are 159 prime 
pairs. There are 113 consecutive composite numbers following 14996687. 

Poletti is the author of Tavole di Numeri Primi entro Limiti Diversi e Tavole A ffini, 
Milan, 1920. 

D. H. L. 


83[G, I].—H. E. Sauzer, Coefficients of the first fifteen General Laguerre 
Polynomials. Ms. in possession of the author. 


The writer announced previously (MTAC, v. 2, p. 89) a manuscript giving the coeffi- 
cients of LAGUERRE polynomials, which are a special case of general Laguerre polynomials 
L,(x), namely for a = 0. The present manuscript gives the polynomials in a which are 
the coefficients of x” in the general Laguerre polynomial 

Lal(x) me ergo 1 1 (4) crane) * 2(* + *) (—x)’ 
n!\dx rag \% — » v! 
and for m=0(1)15. 





, for vy = 0(1)n, 


H. E. SALZER 


AUTOMATIC COMPUTING MACHINERY 


Edited by the Staff of the Machine Development Laboratory of the National Bureau 
of Standards. Correspondence regarding the Section should be directed to Dr. E. W. 
CANNON, 418 South Building, National Bureau of Standards, Washington 25, D. C. 
TECHNICAL DEVELOPMENTS 
Our contribution under this heading, appearing earlier in this issue, is “The California 
Institute of Technology Electric Analogue Computer” by Prof. G. D. McCann. 
DISCUSSIONS 


Procedure for the Machine or Numerical Solution of Ordinary Linear 
Differential Equations for Two-Point Linear Boundary Values 


Introduction. Increased attention is being focused on machine and - 


numerical solutions of differential equations which cannot be solved by 
ordinary mathematical methods. There is need for more information on this 
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subject for engineers and others who deal with such equations. This paper 
describes a procedure applicable to numerical or machine solutions of the 
general non-homogeneous ordinary linear differential equation with variable 
coefficients where the form of these coefficients does not easily permit of 
solution by series. The method is based on the well-known properties of 
linear differential equations. 

Ordinary differential equations of order higher than the first commonly 
describe problems where the known boundary conditions are expressed 
at two different values of the independent variable. Such problems are known 
as two-point boundary-value problems. Although a great many linear equa- 
tions, such as the BEssEL and LEGENDRE equations, may be rigorously 
handled by the method of FroBentus, there are frequently those where the 
variable coefficients of the derivatives are so complex that a series solution 
is not feasible. For such equations, recourse to solution by numerical methods 
or by some type of computing machine or analyzer may be sought. In this 
event, a difficulty is at once encountered if divided boundary conditions are 
present. In the case of an equation of order n, the dependent variable and its 
n — 1 derivatives must possess assigned values at some point within the 
interval in order that a machine or step-by-step solution may proceed from 
that point. Consider the case of a second-order equation where the two 
boundary conditions are divided between both ends of «he range of the 
independent variable. Only one of an unlimited number o/ possible values 
of the initially unknown dependent variable or one of its derivatives, as the 
case may be, at one end of the interval of solution will satisfy the boundary 
condition at the other end. For a fourth-order equation with equally divided 
boundary conditions, a double latitude of possible initial choices would exist. 

The theory of ordinary linear differential equations appears extensively 
in the mathematical literature.’ ** Application to the two-point boundary 
problem where numerical or machine solutions are involved does not appear 
to be generally well known. The purpose of this paper is to show in relatively 
simple mathematical terms and by graphical illustration how the two-point 
boundary problem may be handled. The method applies to any ordinary 
linear differential equation of order m where the boundary conditions are 
expressible in terms of linear combinations of the dependent variable and its 
nm — 1 derivatives and where K boundary conditions are known at one point 
of the interval of solution and the remaining » — K conditions are known 
at some other point of the interval. It is shown that for K < n/2 only 
K + 1 trial solutions with arbitrarily selected values for those derivatives 
which are initially unknown are required to determine the unique solution 
satisfying all m boundary conditions. It is also shown that m + 1 trial solu- 
tions will give data for any solution of the non-homogeneous equation. In 
addition, a possible graphical procedure is suggested for converging on the 
solution of non-linear equations. 

Theory. Consider an ordinary non-homogeneous linear differential equa- 
tion of order m with variable coefficients. This is 


(1) fo(d"y/dx") + fi(d”ty/dx") + +++ + fai(dy/dx) + fry = f. 


Assume fo, fi, -- + fn, f are finite, one-valued, and continuous functions of x 
in the interval a) < x < bo, and that fy does not vanish at any point in the 
interval. Under these conditions, there is known to exist a solution y such 
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that y and (its first derivatives are continuous and have unique values at 
every point in the interval. The reduced or homogeneous equation corre- 
sponding to equation (1) is 


(2) fo(dy/dx") + fild”ty/dx"™) + +--+ + fns(dy/dx) + fay = 0 


and is known to have m and only m linearly independent solutions,‘ ¥;, ye, 
-+ yn. The known complete solution of equation (1) is then 
(3) Y = C11 + Coa + o> Henn +p = Met L Ci, 
i=1 

where y, is any particular solution of equation (1) and ¢, ¢2, --- C, are 
arbitrary constants to be determined by the ” boundary conditions. 

If the quantities B, which assume boundary values are linear,’ we may 
express them in the following manner: 


j=n—1 
(4) B= YX axjdiy/dxi), k=1,2,---n, 
j=0 


where each of the coefficients, a,, ;, represents a constant or some known 
function of x, and where |a;,;| #0 since the boundary conditions are 
linearly independent. Substituting the general solution, equation (3), into 
equation (4) and factoring out the c’s, we have for the particular point, 
x = x, at which B; has a known value, 


i=n j=n—-l1 j=n—l 
(4 = Da L a, (d*yi/dxi) + DV ax, (d*y,/dx*) ) 
i=1 j=0 i=0 pines 


k=1,2,3,---m. 


In this expression for B,, all terms are of fixed value. Since there are m such 
linear expressions for the B’s, and the product | az, ;| -|d*y;/dx*| does not 
vanish, we may solve explicitly for the ” values of the c’s. Carrying this out 
one obtains linear equations of the form, 

kon 
(5) ¢ = b+ DX bi .Bi, 4=1,2,---n. 

k=1 
The term }; is a combination of the constant values of y, and its derivatives 
for the particular x involved, and 0; is a combination of the constants 
dy, ; and the fixed values of y; and its derivatives. Substituting equation (5) 
into equation (3) one obtains 


i=n k=n i=n 
(6) ¥=Vot LV doit LD Bed dinyi. 
i=l k=1 i=] 


This is a formal statement of the fact that the general solution may be ex- 
pressed directly in terms of linear combinations of the » boundary param- 
eters. It now becomes evident that for a given equation of order , if m arbi- 
trary but linearly independent solutions y; of equation (2) and any particular 
solution y, of equation (1) are obtained, then the solution to the problem 
with any desired values of divided boundary parameters, B;, Bo, ---, Bn, 
known at any points within the interval (a, bo), may be obtained by direct 
substitution into equation (6). It should be emphasized that y,; and y, are 
any solutions to their respective equations, barring linearly dependent 
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y;'s, without restriction on initial or final values. It is seen from equation (6) 
that, if the boundary conditions are all homogeneous, then each B equals 
zero and the solution is simply 

i=n 


¥=¥t+ dX divi. 


Also, if the equation is homogeneous but not all of the boundary conditions 
are homogeneous, then 


k=n i=n 
y= DBD dix. 
k=1 i=1 


The term in 5,y; is zero since 5; is a linear combination of the fixed boundary 
values of y, and its m — 1 derivatives, and these quantities exist only in the 
non-homogeneous equation. If the equation and boundary conditions are all 
homogeneous, then it is apparent that y = 0, and there is no problem. 

The advantage of expressing. the solution directly in terms of the bound- 
ary parameters will now be illustrated for a fourth-order equation with 
boundary values B, and By, existing at x = a and B; and B, existing at 
x = b. The limits a, 5 are ordinary points of the interval (do, bo). For n = 4, 
the solution y and its first m — 1 derivatives may be written from equa- 
tion (6) as 

> = BiY; + BY, + B;Y; + BsuY¥,+ Y 
dy/dx = B,(dY;/dx) + etc., d*y/dx? = B,(d*Y,/dx*) + etc., 
d*y/dx* = B,(d* Y;/dx*) + etc., 


where ¥;, Yo, Y3, Y« are linear combinations of y1, ye, ys, ys, and Y is a 
linear combination of the same y’s and y,. Now consider the problem of a 
machine or step-by-step solution starting at x = a where the two boundary 
parameters have the desired values B, = 8, and B, = 6. From the two 
boundary relations at x = a given by equation (4) we may write 


Ai = [ay, oy + @1,1(dy/dx) + 1, 2(d*y/dx") + 1, 3(d*y/dx*) Jena 
Bo = [a2,0y + Ge,1(dy/dx) + de,2(d*y/dx*) + de, 3(d*y/dx*) ],-. 


from which two of the initial values of the derivatives may be solved in 
terms of the remaining two. Thus, for any arbitrary values assigned to 
Yeua and (dy/dx)z-a, the values of the second and third derivatives at x = a 
may be calculated.* With B; and B, assigned the values 8; and f2, respec- 
tively, the solution and its first derivative become 


y = B:Y; + BsY, + Ys, dy/dx = Bd Y;/dx + Bad Y./dx + dY;/dx, 


where Y; is the new function, Y + 6,¥: + 62¥2. At x = a, Y; and dY,/dx 
(where j = 3, 4, 5) assume fixed values so that 


(7) Yuma = d3B; + dB, + ds, (dy/dx) ema = e333 + e4B, + es, 


where the d’s and e’s are constants. Equations (7) show that the initial 
values of y and dy/dx are related linearly to the boundary parameters at 
x = b. This linearity may be represented graphically as shown in Fig. 1. 
Although the existence of these linear families of boundary parameters is 
now established, their determination is still unknown for any trial solution. 
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Assume now that the problem at hand requires that B; = 6; and By = (4 
atx = b. It is readily seen that there exists a point S which defines the proper 
initial choice of yzoa and (dy/dx)z-a. It is also evident from the linearity of 
equations (7) that, in proceeding on a straight line joining any two points, 
such as 1 and 2 in Fig. 1, the values of B; and By, will vary linearly with the 
distance measured along that line. With these facts in mind, let us make 
three separate trial solutions with any arbitrary combination of yz. and 
(dy/dx)2z-2 which define three non-collinear points 1, 2, 3. The initial values 
of the second and third derivatives for each trial are of course determined 
so as to satisfy. B, = 6; and Bz = 8, at x = a. When these trial solutions 
have reached x = 6, the values of y and its three derivatives at this point 
are used to calculate B; and By. By linear interpolation point P; which 
represents a point where B; equals the desired value 8; may be located on 
line 1-2. Likewise points P: and P; which lie on the desired 8; line are located 
by linear interpolation and extrapolation along lines 2-3 and 3-1. Similarly 
points Q1, Q2, Q3 are located on the desired §, line. These two parameter lines 
may now be constructed and the desired solution point S determined by 
their intersection. 






POSSIBLE 
VALUES 
OF B, AT X=b 








Yx=a 


Fic. 1. Control plot for solution of fourth-order linear differential 
equation with equally divided boundary values. 


Although the desired boundary conditions can now be satisfied with a 
fourth solution which begins with the correct combination of yz. and 
(dy/dx) =a, it is not necessary to effect this fourth solution. For any particular 
value of x in the interval (a, 6), the correct values of y and its derivatives 
may be obtained by linear interpolation and extrapolation to the solution 
point S from the corresponding values at points 1, 2, 3. This further linearity 
is at once evident from equations (7) if we replace B; and By, by any two of 
the derivatives of y at this particular value of x. 

Consider next a sixth-order equation with three boundary conditions 
known at x = a and the other three at x = b. Here, there is a triple latitude 
of possible choices of the three initially unknown values of y and its five 
derivatives at x = a and only one unique combination will satisfy the given 
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conditions at x = 6. To solve this problem, we may visualize, in place of 
Fig. 1, three families of parallel planes in a three-dimensional space defined 
by the three initial values of y and its derivatives which are unknown. To 
determine the solution point S representing the common intersection of the 
three boundary-value planes will require four different trial solutions, each 
satisfying the three conditions at x = a. These trials will define four non- 
coplanar points in the three-dimensional space. Linear interpolation and 
extrapolation along any three non-coplanar lines joining these points will 
determine three sets of three points, each set uniquely defining one of the 
desired boundary-value planes. Again the values of y and its derivatives can 
be obtained at point S for any value of x by a three-dimensional linear 
interpolation and extrapolation from the four points. 

In the case of a second-order equation with its two boundary conditions 
divided, only two trials, each satisfying the initial condition at x = a, are 
necessary to establish the solution. The relationship is shown graphically in 
Fig. 2 where the solution point S is determined by the linearity along the line 
joining the trial points 1 and 2. 











o| DESIRED VALUE 
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«> vs y INITIAL VALUE 
<i : ! 
OF 
2 | 
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Fic. 2. Control plot for solution of second order linear differential 
equation with equally divided boundary values. 


So far, only cases of equally divided boundary conditions have been 
discussed. If unequal division occurs, the solution should be started at the 
point where the greater number of known conditions exists. In the case of a 
fourth-order equation where three conditions are known for x = a, only 
one derivative of y is unknown, and two trials will be sufficient for solution. 

Boundary conditions need not exist at the extremities of the interval of 
solution. In such cases the solution may begin at a point within the interval 
where at least half of the conditions are known. It will be necessary to re- 
verse the direction of .the independent variable in order to cover the com- 
plete range, but otherwise the procedure will be the same. 

The following rules may be stated from the foregoing development. 
They hold for any ordinary linear differential equation of order m defined 
by equation (1) where K boundary conditions are known at one point within 
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the interval of solution and the remaining  — K conditions are known at 
some other point within the interval. If K < 2/2, only K + 1 trial solutions 
with K + 1 arbitrarily selected initial values of the initially unknown de- 
rivatives are required to determine uniquely the desired solution which 
satisfies all boundary conditions. If data for the solution for any or all 
possible combinations of the m boundary conditions are required, then any 
n solutions, barring linearly dependent ones, of the homogeneous equation 
without reference to initial conditions and any one solution of the complete 
equation, or a total of m + 1 arbitrary trial solutions, must be made. 

Discussion. Although the m + 1 trial solutions are sufficient to solve all 
boundary value problems for a given equation, in most practical cases, 
where interest is centered on a given set of boundary values, the method 
described using K + 1 trials, where K < n/2, will be the simpler procedure. 
Even in this case K + 1 families of solutions may be had with the K + 1 
trials, since any boundary value at x = 5 for each of the K boundary param- 
eters and any value of } within the interval (ao, b9) may be used. All of these 
solutions must satisfy the same conditions at x = a. 

Graphical illustration has been used to describe the linearities involved 
in order to aid in visualizing the problem. The solution point and the values 
of the various functions at this point may be determined by direct solution 
of the linear relationships involved. And indeed, if it is necessary to discuss 
an equation of order higher than the sixth, our spatial visualization which is 
limited to three dimensions would not aid in this problem. If the method 
of the ~ + 1 trials is adopted, the use of determinants will facilitate the 
necessary computation of b; and b;,, in equation (6). 

There are several practical considerations which place some limits on the 
success of these methods. It is essential that the distances between the trial 
points be of the same order of magnitude as the distance from any point to 
the solution point. Although the danger of interpolation and extrapolation 
on curved lines is absent, still the errors inherent in any machine or numerical 
solution will limit the extent of accurate extrapolation. It is usually possible 
in most physical problems by approximation, comparison, and reasonable 
guessing to predict the general region of solution and thus choose trial 
points which are not too far removed from the solution point. If a reasonably 
close estimate cannot be made, any K + 1 trials will point to the approxi- 
mate location of S whereupon K + 1 additional trials in the neighborhood 
of S will yield the solution. 

In order to insure sufficient accuracy in the results, a reasonable estimate 
of the range of magnitudes of y and its » — 1 derivatives must be made in 
the case of a machine solution. It may be necessary to repeat one or more 
trials with adjusted scale factors if the estimate is far in error. 

Although K + 1 trial solutions are mathematically sufficient for solu- 
tion, one or two more may be desirable in the case of a machine solution 
as a check on the accuracy of the work. In the case of the fourth-order equa- 
tion with equally divided boundary values, the triangular configuration of 
trial points shown in Fig. 1 might well be replaced by four points represent- 
ing the corners of a square. Any three of the six possible line segments 
joining the four points may be considered as locating these points. These 
three segments then determine independently three points on the desired 
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parameter line. The remaining three segments may be considered as de- 
pendent on the first three and hence will determine three dependent points 
on the same parameter line. 

The method described in this paper has been used successfully by the 
author in the solution of a problem in the theory of shells involving a fourth- 
order equation with equally divided boundary conditions. Triangular sets 
of trial points in the machine solution used were adequate. 

In the case of non-linear equations or linear equations where the bound- 
ary conditions are non-linear a procedure similar to that described in the 
foregoing paragraphs is suggested. In this case graphical representation 
would be indicated, and, in a control plot corresponding to Fig. 1, the families 
of lines representing boundary parameters would no longer be straight lines 
or linearly spaced. However, by the principle of uniqueness of solution, it is 
evident that any one boundary parameter line of one family will not cross 
any other boundary parameter line of the other family more than once. 
Also, if a reasonably good initial estimate is made, it should be possible to 
converge on the solution with a few successive sets of trial points. 


J. L. Mertiam 
Univ. of California 


‘ Berkeley 


1E. L. Ince, Ordinary Differential Equations. New York, 1944. 
2A. R. Forsytu, A Treatise on Differential Equations. Sixth ed. London, 1933. 


31.S. & E. S. Sokotnixorr, Higher Mathematics for Engineers and Physicists. New York, 
1941 


‘ The criterion for linear independence of these solutions is that the Wronskian of the 
y's (where i = 1, 2, 3, --- m) and their m — 1 derivatives does not vanish. 

5 Boundary values are said to be linear if they can be expressed as linear combinations 
of the dependent variable and its m — 1 derivatives. 


6 If in some actual problem yz, and (dy/dx),~s, for example, took on assigned boundary 


values, then these may not be changed. Thus (d*y/dx*),.. and (d*y/dx*),.. would be the 
unknown initial quantities. 


BIBLIOGRAPHY Z-IX 


1. ANon., ““AMA—an epoch in telephone accounting,” Bell Laboratories’ 
Record, v. 27, 1949, p. 2-4, illustrs. 25 X 18 cm. 


The automatic message accounting equipment is able to ‘‘read” perfectly 80 coded digits 
per second. Paper tape, on which records of thousands of telephone calls are stored, contains 
information on the calling and called numbers, the month, day, and exact time to tenths 
of minutes at which a telephone conversation begins and ends. The tapes are run through 
“reader”. machines, and the information is assembled, translated, sorted, summarized, and 
printed so that telephone bills may be made from it. 


2. Anon., “Automatic computer for aircraft-motion mathematics,” Jn- 
struments, v. 22, n. 1, Jan. 1949, p. 69, 70, illustr. 20 X 26.5 cm. 


3. Anon., “GE analyzer for Indian Institute,” Electrical Engineering, v. 68, 
1949, p. 84. 29 XK 21.25 cm. 
An a-c network analyzer, for installation at the Indian Institute of Science, Bangalore. 


4. Anon., “Instruments guide industry,’ Business Week, n. 995, Sept. 25, 
1948, p. 70-74, illustr. 25.5 X 28.25 cm. 
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5. Davin L. ARENBERG, “Ultrasonic solid delay lines,’’ Acoustical Soc. 
Amer., Jn., v. 20, 1948, p. 1-26, bibl., diags., tables. 19 K 27 cm. 


Abstract: The development of an ultrasonic delay line for the storage of pulses one 
microsecond long for periods up to 2000 microseconds is discussed. Theoretical considera- 
tions of the piezoelectric transducer involving the quartz crystals and acoustical constants 
of the medium developed by EMsLiE, HUNTINGTON, and BENFIELD are applied to determine 
the band width and power loss of the delay line. Problems, peculiar to the use of solids be- 
cause of the presence of more than one mode of propagation and the distribution of energy 
between the different modes on reflection from an air surface, are evaluated for a wide range 
of Poisson ratios. Measurements of velocity and attenuation in some polycrystalline metals, 
single crystals, and glasses are given. Fused quartz was found to be the most promising 
material. 

Delay lines constructed using a light source and photoelastic pick-up were not com- 
pletely developed, but several types using crystal as well as transmitter pick-up were satis- 
factory. For delay times longer than 300 microseconds, multiple reflection paths in two and 
three dimensions were used and their preparation from fused quartz discussed. 


6. A. R. BootHroyp & E. CoLtn CHERRY, “Analogue calculating machine for 
functions of a complex variable,” Nature, v. 163, April 30, 1949, p. 687. 
17.4 X 25.4 cm. 


7. J. G. BRAINERD & T. K. SHARPLEsS, ‘““The ENIAC,” Electrical Engineer- 
ing, v. 67, n. 2, Feb. 1948, p. 168-172, diags. 21.5 XK 29 cm. 


8. J. A. Bronzo & H. G. CouEn, “A note on analog computer design,” 
Rev. Sci. Instrs., v. 20, n. 2, Feb. 1949, p. 101-102. 19 K 26.5 cm. 


A method is proposed for the simplification of computers designed to solve certain 
types of partial differential equations. 


9. C. Cuu, “The recent development in the electronic digital variable com- 
puter,” Chinese Institute of Engineers, Jn., v. 6, no. 1, May 1948, 
p. 58-65. 14.5 X 23 cm. 


10. Harry M. Davis, ‘‘Mathematical machines,” Scientific American, v. 
180, April 1949, p. 28-39, illustrs. 21.25 X 29 cm. 


This article is principally concerned with the mathematical machines which count, i.e., 
the digital machines, and a detailed description of the digital idea is presented. It em- 
phasizes the fact that, although these machines are powerless to do anything creative, they 
are extremely useful in that they are not subject to human intervention and are capable of 
enormous speed. Current large-scale computer developments are mentioned including: (1) 
the IBM Selective Sequence Calculator and the electronic Type 604 Calculator, which is 
ten times faster than the previous IBM commercial calculating punches; (2) the Moore 
School EDVAC, having a far more capacious memory and greater versatility than its fore- 
runner, the ENIAC; (3) the twin computers known as BINAC and the UNIVAC, developed 
by John Mauchly and J. P. Eckert, JR.; (4) the colossal computing device which is being 
designed at the Institute for Advanced Study for purposes of weather prediction and pos- 
sibly for work on questions of atomic physics and economic statistics; (5) the Raytheon 
Manufacturing Company computer; (6) the Harvard Mark III computer; (7) the MIT 
“Project Whirlwind” for the construction of a computer for military purposes; and (8) 
advanced computers now under construction by four British Laboratories. Every current 
computer will soon become outdated in this new age of electronic thinking; however, the 
older machines will not become obsolete, because problems exist suitable for every type to 
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handle. Also the relay, as a computing device, and the much speedier electron tube is dis- 
cussed with a detailed explanation of the operation of the 20 ENIAC accumulators and of 
the binary representation of numbers. The application of the binary system to logic as well 
as numbers is stressed. In particular the author mentions the various memory devices now 
being developed and outlines future applications to the 1950 census, adjustment of maps, 
engineering computations, hydrodynamics, large-scale planning in the Air Force, etc. 


11. HERMAN H. GOLDSTINE & JOHN VON NEUMANN, Planning and Coding of 
Problems for an Electronic Computing Instrument. Part I1, v. 3. Princeton, 
N. J. Institute for Advanced Study, 1948, iii + 23 p. 21.6 K 27.9 cm. 


This volume deals with the manner of incorporating subroutines into a complete prob- 
lem. After considerable explanation, the authors conclude that this is best accomplished for 
small subroutines by recording them into the main routine as needed. For larger sub- 
routines a preparatory routine automatically modifies the addresses of the subroutine 
relating to positions within itself by adding a quantity corresponding to the difference 
between the actual location of the subroutine in the memory and the assumed location 
at the time the subroutine was initially coded. There is a detailed discussion of suitable 
preparatory routines for a single subroutine and for a group of subroutines all of which 
are considered as coding problems. The volume closes with a description of proposed 
details for handling main routines and subroutines and the preparatory routine for machine 
input. 

While the subject matter is of appreciable importance and the manner of handling it 
shows ingenuity, there is some question as to whether its treatment in this report is optimum. 
The following are some of the points one might question: 


1. In order to use the method described, the coding of subroutines must be handicapped 
somewhat by certain limitations. For example, it is not permissible to use one half of a word 
for numerical storage and the other half as an order. It is also not permissible to build up 
addresses by the addition of numerical quantities, some of which occupy address positions, 
unless these quantities satisfy certain requirements. While such limitations are simple in 
character, their existence lends some importance to use of other methods which impose 
fewer restrictions to achieve the same result. 

In addition to other methods which utilize auxiliary input equipment (for example, that 
suggested by SAMUEL LUBKIN in “Proposed programing for the EDVAC,” Aberdeen 
Proving Ground, Jan. 1947), we may mention one suggested by ALAN L. LEINER and 
RALPH J. SLUTZ of the Electronic Computer Section of NBS. It is based upon the provision 
of extra binary digits in the coding as an indication of whether the corresponding address 


is to be taken as an absolute location in the memory or as “relative’’ to the memory location 
in which it occurs. Also, there is a method based upon the provision of one or more “‘base 
numbers” for automatically modifying appropriate classes of addresses (suggested by 
Samuel Lubkin primarily for utilization in connection with iterations but incidentally of 
value in connection with the present problem). In addition, it is possible to use a different 
criterion than that chosen by the authors for determining which addresses are to be modi- 
fied. For example, NBSMDL has recently coded an analogous preparatory routine for the 
EDVAC, using the first binary digit in each address to indicate whether or not the address 
requires modification. 

2. Time requirements for machine modification of subroutines are certainly negligible 
compared to those for manual operations relating to insertion of subroutines from a “‘library” 
into the memory. This does not, however, necessarily rule out special routine-preparing 
equipment, since the latter may facilitate the manual operations themselves and may make 
coding more convenient. 

3. There seems to be little reason for not including the constants of the preparatory 
routine in the main routine instead of inserting each constant manually as indicated by the 
authors. Also, the space taken up by the preparatory routine might be used by the main 
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routine at some time after the address modification is completed. For this purpose, it might 
be more convenient to locate the preparatory routine at the beginning of the memory rather 
than at the end. After reading in the subroutines and the preparatory routine (with con- 
stants), the latter could be carried out and followed by instructions to read in the balance 
of the main routine, re-using memory locations originally containing the preparatory routine. 

4. The authors mention the possible convention of having all subroutines coded as if 
they begin at 0, but they do not require this in their preparatory routine. Such a convention 
seems simpler to use than the requirement that the actual position occupied by a sub- 
routine be further along in the memory than the position for which it is coded. Use of the 
convention would cut the length of the preparatory routine. It might also be noted that one 
of the subroutines (by choice, the longest) can always be inserted in the memory in the same 
position for which it was originally coded and, in such a case, requires no modification. 

5. The coding given by the authors for the preparatory routine can be materially 
shortened. A new code of the routine has been prepared by Otto T. STEINER & Samuel 
Lubkin requiring only 36 + 2I ‘‘words” (where J represents the number of subroutines) as 
compared to the 58 + 4J ‘‘words’’ required in the original routine. It was prepared without 
the use of flow diagrams or sub-sectioning. Indeed, most MDL coders find such aids, in the 
form used by the authors, confusing and even detrimental. 

SAMUEL LUBKIN 
NBSMDL 


12. D. R. HARTREE, ‘‘A historical survey of digital computing machines” ; 
M. H. A. Newman, ‘General principles of the design of all-purpose 
computing machines”; M. V. WILKEs, ‘“The design of a practical high- 
speed computing machine. The EDSAC”’; F. C. WILLIAMs, ‘‘A cathode- 
ray tube digit store’; J. H. WILKINSON, “The automatic computing 
engine at the National Physical Laboratory’; A. D. Bootn, ‘Recent 
computer projects,’ R. Soc. London, Proc., v. 195A, Dec. 1948, p. 
265-287 + 3 plates. Discussion held 4 Mar. 1948. 17.2 X 25.6 cm. 
Reported in Nature, v. 161, May 8, 1948, p. 712-713. 15.5 X 25.5 cm. 


The discussion was confined to automatic digital general-purpose computing machines. 
The EDSAC (Electronic Delay Storage Automatic Calculator) is at the University Mathe- 
matical Laboratory, Cambridge, England. Compare MTAC, v. 3, p. 214-215. 


13. Wm. A. McCoo., An Improved Electronic Analog Computing Circuit. 
Report no. P-3423. Washington, Naval Res. Lab., Feb. 23, 1949, 15 p. 
illustr. 26.5 XK 20 cm. 


Abstract: An improved electronic analog computing circuit, employing only two vacuum 
tubes and characterized by both simplicity and accuracy of computation, is described and 
analyzed. The theoretical error is calculated for integration, differentiation, and constant 
multiplication with step-function and sine-wave forms of excitation. The results of a cursory 
experiment, made to evaluate the performance of the improved circuit functioning as a 
differentiator, agree with the theoretical results. 


14. T. K. SHArRPLEss, ‘Design of mercury delay lines,’ Electronics, v. 20, 
Nov. 1947, p. 134-138, illustrs. 28.25 20.25 cm. 


Mercury is used to obtain millisecond delays for radar, computers, and memory devices 
because it transmits compression waves relatively slowly, introduces negligible loss, and has 
an impedance comparable to that of crystal transducers. Recirculation and temperature 
compensation techniques are presented. 
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15. NORBERT WIENER, Cybernetics. Hermann et Cie, Paris, and John Wiley 
& Sons, New York, 1948, 194 p. 22.5 K 14.5 cm. 


The term, cybernetics, has been chosen by Wiener and his collaborators to designate 
the entire field of control and communication in the animal and the machine. The justifica- 
tion for a common treatment of apparently diverse phenomena in control and communica- 
tion, ranging, for example, from the operation of purely mechanical feed-back systems to 
neurological phenomena like memory and visual perception, would appear to lie in the 
identity of the mathematical approach. However, as is explained in a thirty-nine page 
narrative which serves as an introduction to the text of the book, the underlying motif of 
research in cybernetics is rather the investigation of boundary regions of science, areas 
formerly “‘neglected as a no-man’s land between the various established fields.” 

The volume, Cybernetics, contains eight formal chapters, beginninz with a chapter on 
Newtonian and Bergsonian time and ending with a chapter on information. language, and 
society. The remaining chapters are concerned, respectively, with ‘Groups and statistical 
mechanics,”’ “Time series, information, and communications," “‘Feed-back and oscillation,” 
“Computing machines and the nervous system,” ““Gestalt and universals,”’ and “Cybernetics 
and psychopathology.” In addition there is a note, refreshingly objective by comparison 
with popular press notices on the subject, in which the capability of a computing machine 
at chess playing is discussed briefly. 

The introduction gives a somewhat detailed account of the origin of Wiener’s interest in 
the investigation of the nature of control and communication in the human being and of the 
manner in which group effort was developed in the field. It is stated that the book represents 
the outcome of a program of work undertaken jointly with Dr. AkRTURO RosENBLUETH, 
formerly of the Harvard Medical School, and now of the Instituto Nacional de Cardiologia 
of Mexico. The assistance of scientists, distinguished for their outstanding work in such 
diverse fields as electro-mechanics, anthropology, sociology, and statistical opinion sampling, 
is acknowledged. 

In chapter I, ‘“‘Newtonian and Bergsonian time,” the implications of reversibility of 
time are discussed. It is pointed out that planned physical experiments introduce asym- 
metry into reversible time systems. In fact, the very nature of the questions asked, of prob- 
ability and prediction, must of necessity be asymmetrical, distinguishing between the past 
held in the present by the fixing of certain quantities, and the future. Newtonian astronomy 
and meteorology, Newtonian dynamics and Gibbsian statistical mechanics, physics, and 
biology are discussed from the standpoint of reversibility of time. The difference between 
the irreversible time of evolution and biology and the reversible time of physics was em- 
phasized by Bergson. In Wiener’s opinion, the controversy between vitalism and mechanism 
had at its kernel the fact that Newtonian physics is not the proper frame for biology. 

Chapter II, “Groups and statistical mechanics,” is, in a sense, an introduction to the 
following chapter, ‘“‘Time series, information, and communications."’ These two chapters 
alone, while they are exacting from the standpoint of the mathematical maturity and back- 
ground required of the reader, would appear to justify the cost of publishing the book. 
For example, a lucid yet concise account is presented of the foundations of statistical me- 
chanics, including the role of LEBESGUE measure and the metrical invariants of transforma- 
tion groups in the development of the ergodic theory, and a discussion of the notion of 
entropy. Also a satisfying formal exposition of a statistical theory of time series is made 
showing how the theory applies to message-noise problems in communication, to the en- 
semble of time series associated with Brownian motion, and to prediction problems. Single 
and multiple time series, both continuous and discrete, are treated from the viewpoint of 
full knowledge of the parts of the series. The role of the theory of time series in quantum 
mechanics is touched upon and the line of development of a practicable theory of time series 
involving conclusions based on a sampling of the past is indicated. 

In the remaining chapters of Cybernetics, IV, ‘‘Feed-back and oscillation,”’ V, ‘““Com- 


puting machines and the nervous system,”’ VI, ‘Gestalt and universals,”’ VII, “Cybernetics 
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and psychopathology,” and VIII, ‘Information, language, and society,” the reader of 
moderate mathematical interest will find the going easier. It is necessary, says the author, 
to use mathematical symbolism and mathematical techniques in IV, but the mathematical 
treatment is made more palatable to the average reader by extensive verbal explanation. 
In IV, feed-backs found in human and animal reflexes and homostasis are discussed. In 
chapters V, VI, and VII, discussion of both normal and abnormal action of the nervous 
system centers about the ultra-fast, automatically-sequenced computing machine as a 
model. The reviewer cannot repress a feeling that the author is unrealistic in his discussion 
of possible future lines of development of such computing devices, for example, the incor- 
poration in them of the ability to learn; nevertheless, the utility of such devices outside the 
field of numerical computation is made clear. 

The last chapter, ‘Information, language, and society,’ ends the Cybernetics on a 
sombre note, seeming to possess a tincture of misanthropy. One can only hope that in its 
impact upon sociology, the new science of cybernetics will exert a beneficent influence and 
that the social sciences will not be too greatly retarded in development by the comparatively 
close coupling of the observer with phenomena mentioned in this chapter. 


E. W. CANNON 
NBS 


EpitoriAL Norte: Professor Wiener’s book was printed in France and is filled with 
typographical errors, sometimes making the mathematics quite unintelligible. 


NEws 

Association for Computing Machinery.—Cuar_es B. Tompkins has resigned from the 
Council of the Association, and JAN RAJCHMAN has been appointed. At present, the member- 
ship totals about 600. Prof. Evy VELANDER of Sweden has been authorized to undertake the 
formation of a Swedish section, and Mr. ALBERT CAHN of NBSINA is in charge of forming a 
California section of the Association. 

About 200 members of the Association attended the meeting in Oak Ridge, Tennessee, 
April 18-20, 1949. Those who were present found it to be an informative and satisfying 
meeting. The Association greatly appreciated the welcome and hospitality of the Oak Ridge 
National Laboratory, the Oak Ridge Institute of Nuclear Studies, and the other Oak Ridge 
groups. 

The papers listed below were delivered. Plans are being made to furnish summaries (or 
the full version) of these papers to members. In some cases the full paper will be submitted 
to the editors of MTAC. 

1. ‘“‘Address of welcome,’’ A. M. WEINBERG, Oak Ridge Nat. Lab. 2. ‘Limitations of 
electronic digital computing machines,’’ E. W. Cannon, NBS. 3. “Analog and special- 
purpose computing machines,’’ S. H. CALDWELL, MIT. 4. “Inversion of large matrices,” 
GERTRUDE BLANCcH, NBSINA. 5. “Problems in linear programing,” Emit SCHELL, 
Office of the Air Comptroller, USAF, Washington. 6. “‘Computing machinery in Australia,”’ 
T. G. Room, Univ. of Sidney, Australia. 7. “The forgotten man of computing,”” MINA REEs, 
ONR, Washington. 8. ‘Statistical theory of round-off errors,”” H. D. Huskey, NBSINA. 
9. “The relations of symbolic logic and large-scale calculating machines,’ EpMuND C. 
BERKELEY, New York. 10. ‘‘Calculation of neutron age by Monte Carlo,’”’ A. S. House- 
HOLDER, Oak Ridge Nat. Lab. 11. ““Computing statistical atomic fields on the selective se- 
quence electronic calculator,’’ L. H. THomas, IBM. 12. “Binary notation in arithmetic 
algorithms and operation code partitions,’’ B. F. CHEypLEUR, Naval Ordn. Lab., White 
Oak, Silver Spring, Md. 13. “602A calculating punch—time saver in computing thermo- 
dynamic properties of gases,"” JACK Bevzer, H. L. Jounston, & Lypta SAveporr, Dept. 
Chem., Ohio State Univ., Columbus, O. 14. “The register code for the ENIAC,” R. F. 
CLiPPINGER & B. DimspALeE, BRL, Aberdeen, Md. 15. ‘‘Recent progress on the Selectron,” 
JAN RAJCHMAN; Radio Corp. Amer., Princeton, N. J. 16. “Analog computer for FoURIER 
transforms in crystal analysis,” RAYMOND PEepinsky, Alabama Polytechnic Inst., Auburn, 
Ala. 17. “The philosophy of computing machine design,’’ GEoRGE R. Stisitz, Burlington, 
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Vt. 18. “Internal conversion calculation on Mark I,” M. E. Rose, Oak Ridge Nat. Lab., & 
BERNARD SpInRAD, Argonne Nat. Lab., Chicago. 19. “Experience with the use of the 
AIKEN Relay Calculator,”” CLARENCE Ross, Naval Proving Ground, Dahlgren, Va. 20. 
“Logical design of direct-coupled vacuum tube control circuits,"” THEoporRE Kain, Com- 
putation Lab., Harvard Univ. 21. “‘The statistical mechanics of high polymers,” GiLBert W. 
Kinc, Cambridge, Mass. 22. “The NBS Interim Computer,” Samuet Luskin, NBS. 23. 
“Uses of pulse transformer in a-c coupling of computers,”” RALPH SLutTz, Electronics Div., 
NBS. 24. “‘An intercomparison of storage devices,” S. N. ALEXANDER, NBS. 

NBSINA.—Two symposia were held on June 22 through July 1, 1949, on subjects 
pertinent to the effective utilization of automatic digital computing machinery. A representa- 
tive group of the nation’s engineers, physicists, and mathematicians participated. 

The subject of the first one, June 22-25, was the Construction and Applications of 
Conformal Maps. The applications of conformal maps in such fields as aerodynamics and 
electronics were emphasized, and special attention was devoted to the actual current needs 
of research workers in these fields. Construction methods were discussed, along with the 
problem of programing such methods on existing and proposed automatic digital computing 
systems. In this connection, particular reference was made to the electronic machine now 
being designed at NBSINA under the direction of Harry D. Huskey. 

The program opened with a one-day course designed to acquaint those attending the 
Symposium with the preparation of problems for automatic digital computing machines. 
The following talks were included: ‘‘Definition of an automatic digital computing machine” 
by H. D. Huskey; “Description of a specific automatic computer” by ROSELYN SIEGEL, 
NBSINA; “Programing the solution of m simultaneous linear equations” by H. D. 
Huskey; and a laboratory session in which participants were urged to work out detailed 
routines for causing the computers to perform such operations as division, floating addition, 
solving a simple differential equation, and possibly solving a set of simultaneous linear equa- 
tions. On subsequent days the following Sessions were held: 


I. June 23, General Session, J. H. Curtiss, NBS, chairman: 

“On network methods in conformal mapping computation”’ by R. von Mises, Harvard 
Univ. 

“Conformal mapping of domains of higher topological structure illustrated by flow 
patterns’ by R. Courant, New York Univ. 


II. June 23, Session on Physical and Industrial Applications, D. V. Wipper, Harvard 

Univ., chairman: 

“Applications of conformal mapping to torsional rigidity, principal frequency, and 
electrostatic capacity”’ by G. SzeG6, Stanford Univ. 

“Some industrial applications of conformal mapping” by H. Poritsky, General 
Electric Co. 

“Conformal maps involving multiple-connected regions and their technical applica- 
tions” by G. Stern, Westinghouse Electric Corp. 

“On the use of conformal mapping in problems of two-dimensional elasticity” by I. S. 
SOKOLNIKOFF, UCLA. 


III. June 24, Session on Fluid Dynamics, J. L. BARNEs, UCLA, chairman: 

“On the HELMHOLTZ problem of conformal representation’”’ by A. WEINSTEIN, U. S. 
Naval Ordn. Lab. and Univ. of Maryland. 

“Aspects of conformal mapping in aerodynamics” by I. E. Garrick, NACA. 

“On THEODORSEN’s method of conformal mapping’ by ALEXANDER OsTROWSK!, 
Univ. of Basle and NBS. 

“On conformal mapping of variable regions” by S. E. WARSCHAWSKI, Univ. of Minne- 
sota and NBS. 

“Fluid dynamics, conformal mapping, and numerical methods” by ANDREW VAZSONYI, 
U. S. Naval Ordn. Test Sta. Pasadena, Cal. 
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IV. June 24, Session on Theory of Conformal Maps, W. T. Martin, MIT and UCLA, 

chairman: 

“Some remarks on variational methods applicable to multiple-connected domains” by 
D. C. SPENCER, Stanford Univ. 

“A variational method for simply-connected domains” by A. C. SCHAEFFER, Purdue 
Univ. 

“Kernel functions and conformal mapping’? by STEFAN BERGMAN, Harvard Univ., & 
MENAHEM SCHIFFER, Stanford Univ. 

“The kernel function and canonical conformal maps” by ZEEV NEHARI, Washington 
Univ. 

“A new proof of the Riemann mapping theorem” by P. R. GARABEDIAN, Stanford Univ. 


V. June 25, General Session, J. W. GREEN, UCLA, chairman: 
“Conformal mapping applied to electromagnetic field problems’ 
Polytechnic Inst., Brooklyn. 
“An extremal method in conformal mapping’”’ by L. V. AHLFors, Harvard Univ. 
“Some generalizations of conformal mappings occurring in gas dynamics” by L1ipMAN 
Bers, Syracuse Univ. 
“The use of conformal mapping in the study of flow phenomena at the free surface of an 
infinite sea’’ by EUGENE P. Cooper, U. S. Naval Ordnance Test Sta., Pasadena, Cal. 
“Recent contributions of the Hungarian school to conformal mapping” by G. SzEG6, 
Stanford Univ. 


’ 


by Ernst WEBER, 


VI. June 25, Session on Numerical Methods, H. F. BoHNENBLUsT, Cal. Inst. Techn., 
chairman: 

“Relaxation methods as ancillary techniques” by Sir RICHARD SOUTHWELL, London, 
England. 

“The use of conformal mapping to compute flows with free streamline’ by D. M. 
Youn, Harvard Univ. 

“An approximation method for conformal maps” by LEE H. Swinrorp, Univ. of Cali- 
fornia, Berkeley. 

The second Symposium dealt with Probability Methods in Numerical Analysis, with 
special reference to the techniques now known among physicists by the name, ‘“‘ Monte 
Carlo.” It was held on June 29-July 1, under the joint sponsorship of NBSINA and the 
Ranp Corp. with the cooperation of the Oak Ridge Nat. Lab. 

A word of explanation concerning the ‘‘Monte Carlo’’ technique is in order here. It is 
well known in the theory of stochastic processes that the probability distributions associated 
with certain random walk or random flight problems satisfy certain classical integro- 
differential equations. This fact provides a novel method for approximate numerical integra- 
tion of such equations through the medium of building up large samples of trials of 
corresponding random walks. It also suggests the possibility and desirability of directly con- 
structing appropriate random walks as discrete mathematical models for given physical 
situations, to be used in place of the more classical continuous models. This method of solu- 
tion of problems in mathematical physics by sampling techniques based on random walk 
models constitutes what is known as the “Monte Carlo” method. The method as well as the 
name for it were apparently first suggested by JOHN VON NEUMANN and S. M. ULam. 
“Monte Carlo” techniques are now being used rather extensively in connection with high- 
speed automatic digital computing machinery. 

Many mathematical and physical problems remain to be solved in the theory of ‘‘Monte 
Carlo” techniques. Most of the current ‘Monte Carlo’’ applications are classified; how- 
ever, the proceedings of the Symposium were unclassified. The program was as follows: 


I, June 29, Orientation Session, J. H. Curtiss, chairman: 


“The Monte Carlo method” by S. ULAm, Los Alamos Sci. Lab. 
“Discussion of the Monte Carlo method” by J. von NEUMANN, Inst. Adv. Study 
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II. June 29, Session on Physical Applications of Stochastic Methods, J. von NEUMANN, 
chairman: 

“Introductory remarks” by the chairman 

“Showers produced by low energy electrons and photons” by Rosert R. WILSON, 
Cornell Univ. 

“Nomograms for Monte Carlo solution of the Milne problem” by B. I. Sprnrap, 
Argonne Nat. Lab. and G. H. GortzeL_, New York Univ. 

“Neutron age calculations in water, graphite, and tissue’ by A. S. HouSEHOLDER, Oak 
Ridge Nat. Lab. 

“Methods of probabilities in chains applied to particle transmission through matter’ 
by Lewis NEtson, Oak Ridge Nat. Lab., & WeNpELL DeMarcus, Fairchild Engine and 
Airplane Corp. 

““Multiple-scattered gamma rays” by UGo Fano, NBS 

“Stochastic methods in statistical mechanics” by GiLBert W. Kinc, Arthur D. 
Little Co. 


III. June 30, Session on Physical Applications of Stochastic Methods (continued), FRANK 
C. Hoyt, Argonne Nat. Lab., chairman: 


“Calculations on a water shield for fast neutrons, I’’ by Marta Mayer, Argonne Nat. 
Lab. 

“Calculations on a water shield for fast neutrons, II’’ by Preston HAMMER, Los Alamos 
Sci. Lab. 

“A Monte Carlo technique for estimating particle attenuation in bulk matter” by B. A. 
SHoor, Northrop Aircraft; Lewis NELtson, Oak Ridge Nat. Lab.; WeENDELL DEMArcus, 
& Ropert L. Ecuors, Fairchild Engine and Airplane Corp. 

“Neutron transmission through thick slabs’”’ by WILLIston Suor, U. S. Navy. 

“Estimation of particle transmission by random sampling” by HERMAN KAHN, RAND 
Corp. & T. E. Harris, Ranp Corp. 


IV. JuNE 30, Session on Random Digits, JeERzy NEYMAN, Univ. of California, Berkeley, 
chairman: 
“‘Generation and testing’ by GEoRGE W. Brown, RAnp Corp.; NicHOLAS METROPOLIS, 
Los Alamos Sci. Lab.; & GEorGE E. ForsytHe, NBS. 
“Various techniques used in connection with random digits’’ by ]. von NEUMANN. 
General discussion from the floor. 


V. Jury 1, Mathematical Session, JoHn W. TUKEY, Princeton Univ., chairman: 

“The connection between stochastic processes and partial differential equations” by 
Wit Fetter, Cornell Univ. and NBSINA; Mark Kac, Cornell Univ. and NBS; & J. L. 
Doos, Univ. of Illinois. 

Discussion from the floor. 


VI. July 1, Round Table Discussion: Critique of Applications and Discussion of Possible 
New Directions for Research, JoHN W. TUKEY, discussion leader: 


Panel: J. NEYMAN; M. A. Girscuick, Stanford Univ.; W. A. WALLIs, Univ. of Chicago; 
D. H. BLAcKwELL, Howard Univ.; T. W. ANDERson, Columbia Univ.; & L. J. Savace, 
Univ. of Chicago. 

Also on June 27-28, a condensed, relatively elementary course on automatic computa- 
tion was held at NBSINA. The purpose of this course was to introduce interested persons 
to the logical theory and performance characteristics of automatic digital calculators being 
developed at the present time. Preparation of problems for the calculators was emphasized. 
It was composed of two lecture and two laboratory sessions. 
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BIBLIOGRAPHY Z-IX 


1. R. T. Brown, “Notes on the use of the log-log slide rule in connection 
with trigonometric functions,”’ Jn. Engin. Educ., v. 38, 1948, p. 393-394. 
io Me 2e.0 CM. 


2. A. E. CARTER & D. H. SADLER, ““The application of the National Account- 
ing Machine to the solution of first-order differential equations,”’ Quart. 
Jn. Mech. Appl. Math., v. 1, 1948, p. 433-441. 15.4 X 23.4 cm. 


Summary: MILNe’s formula for approximate quadrature is used as the basis of a method 
for the solution of first-order differential equations on the National machine. The method, 
which is illustrated by a numerical example, enables the machine to form the required 
dependent variable without the necessity for conversion from a sum to an integral. 


3. HaroLtp P. Knauss, “Slide rule calculations of radioactive decay,” 
Science, v. 107, 26 Mar. 1948, p. 324. 


Selected sentences: The decay of a radio-active substance is described by the equation 
N = Noe = Nol0~’!, in which N is the number of particles remaining at the time ¢; No 
is the initial number of the particles; \, the decay constant or fraction of the number 
present disintegrating per unit of time; and \’ = .4343d. On the log log rule a single setting 
of the slide gives the values of N/No for the time. Three significant figures are obtained 
with a 10” slide rule. 


4. H. W. RicHarps & D. H. SADLER, “The application of the National 
Accounting Machine to the calculation of apparent places of stars,” 
RAS, Mo. Not., v. 108, 1948, p. 154-161 + a folding tabular plate. 
17 X 25.1 cm. 


Summary: A method is developed for the systematic calculation of apparent places 
of stars, using the National machine. Checked copy, printed by the machine in precisely 
the form required by the printer, is produced in a time which compares favourably with 
other methods. 


5. V. VAND, “A mechanical calculating machine for X-ray structure factors,”’ 
Nature, v. 163, 29 Jan. 1949, p. 169-170. 17.4 & 25.5 cm. 


First sentences: A large mechanical machine for calculation of X-ray crystallographic 
structure factors has now been completed in our laboratories [Lever Bros. and Unilever 
Ltd., Port Sunlight, Cheshire] and is running satisfactorily. The machine is of a tide- 
predictor type, and it deals with up to 24 harmonic components at a time. The expression 
calculated is F(hkl) = = fi cos 24(hx; + kyi + Is;). 

. 


6. HEINZ WITTKE, Die Rechenmaschine und thre Rechentechnik. Eine 
Einfiihrung und ein Ubungsbuch mit ausgewahlten Anwendungsbeis pielen 
aus der Geoddsie, Geometrie und angewandten Mathematik. 2. Auflage. 
Mit 52 Abbildungen, 12 Formularen und 1 Tabellenanhang. (Sammlung 
Wichmann, v. 12). Berlin, Herbert Wichmann, 1948, viii, 160 p. 17.4 
X 24.4 cm. “Verdéffentlicht unter Lizenz-Nr. 136 der Sowjetischen 
Militaérverwaltung in Deutschland.” 


This is the so called second edition of which we listed the first edition in MTAC, v. 3, 
p. 390. The first edition was neatly bound in full boards. The second edition is a poor paper 
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bound offset print of the first edition, authorized by license of the Soviet military powers. 
The first 70 p. include a description of various desk calculating machines, of problems with 
which they may deal, and also lists of references to literature appearing before 1942. 


» Ga 
NOTES 


106. GEORGE NEVILLE WATSON—TABLE MAKER.—This distinguished 
British Mathematician and Table Maker was born in Westward Ho, 
Devonshire, 31 January 1886. He was a student at Trinity College, Cam- 
bridge, senior wrangler, 1907; class I (div. ii), mathematical tripos, part I], 
1908; Smith’s Prizeman, 1909; fellow of Trinity, 1910-1916. 

He was an assistant professor of pure mathematics at University Col- 
lege, London, 1915-18; professor of mathematics at the University of 
Birmingham since 1918, of pure mathematics since 1937; gold medallist of 
the Danish Royal Academy of Sciences, 1912; a fellow of the Royal Society 
of London, 1919; honorary secretary, 1919-33, and president, 1933-35, of 
the London Mathematical Society ; president of the Mathematical Associa- 
tion, 1932-34; SYLVESTER medallist of the Royal Society, 1946; and Dr- 
MorGAN medallist of the London Mathematical Society, 1947. He is also a 
Sc.D. of Univ. of Cambridge, D.Sc. Univ. of London, M.Sc. Univ. of Bir- 
mingham, Hon. LL.D. Univ. of Edinburgh, Hon. Sc.D. Univ. of Dublin. 

The following is a list of Professor Watson’s tables, and of some papers 
with numerical results of some importance : 


1. (with G. A. Schott) “Asymptotic formulae occurring in electron theory,” Quart. Jn. 
Math., v. 47, 1916, p. 311-333. Various numerical results are given. 
2. “The sum of a series of cosecants,”” Phil. Mag., s. 6, v. 31, 1916, p. 111-118. 


n—l 
ys 


S, = 2 csc(mm/n) is tabulated for mn = [2(1)30(5)100, 360, 1000; 5D]. 


m=1 

3. [Tables connected with gamma functions], BAAS, Report, 1916, p. 123-126. Four 
10D tables: (a) 10 + In (1 + x), x = .005(.005)1; (b) 10 + SJ log r(1 + Adt, x = .01- 
(.01)1; (c) ¥(x) = din I'(x)/dx, x = 1(1)101; ¥(x) for x = 1.5(1)100.5. 

4. “The zeros of Bessel functions,” R. Soc. London, Proc., v. 94A, 1918, p. 190-206. 
Tables of J+4(x), U3(x), — Vi(x), — Vi(x)/U3(x), for x = [0(.05)2(.2)8; 4D]. 

5. “Bessel functions of equal order and argument,” Phil. Mag., s. 6, v. 35, 1918, p. 
364-370. Table of nJ/i'J,(nx)dx, n = [1(2)23; 7D]. 

6. A Treatise on the Theory of Bessel Functions, Cambridge, England, 1922; second ed., 
New York, 1944; reprinted Feb. 1945, Apr. 1948; tables, p. 665-752. For tabular errors in this 
v. see MTAC, v. 2, p. 49-51. 

T. I: Jo(x), Yo(x), Ji(x) and Yi(x), for x = [0(.02)16; 7D]. The values of Jo(x), Ji(x) 
up to 15.5 were taken from MEIsseEL’s 12D table (1889) while the values of Yo(x) and 
Y,(x) were computed partly by interpolation in ALpts’ table of Go(x) and G,(x) (1900). 
T. I gives also, for the same range of argument, 7D values of |H,"(x)|, H,(x), » = 0, 1, 
and of arg H,(x), m = 0, 1 to the nearest 0’’.01. 

T. II consists of tables of e~*Io(x), e~*Ji(x), e*Ko(x), and e*K,(x), e*, for x = [0(.02)16; 
7D]. The 8S or 9S table of e* was constructed with the help of Newman’s 12-18D table of 
e* (1883). 

T. III consists of J3(x), Y3(x), | Hy (x)|, and arg H}(x), of the same scope as T. I; 
a table of e*K}({x) is also included. 

T. IV gives 7D values of J,(x), for m=2(1)5, x=.1(.1)5; 6D values of J,,(x), m=0(1)20, 
for x = 1(1)12; 7S at least, or 7D values of Y,(x) for nm = 0(1)10, and x = 0(.1)5, and 7D for 
n = 0(1)13, x = 6(1)12; 7D values of e*J,,(x) for m = 2(1)5, x = .1(.1)5; 7S at least, or 7D 
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values of K,,(x) for m = 0(1)10, x = 0(.1)5. The 6D values of J,,(x) are taken from LoMMEL’s 
table (1885) with some corrections, and some values of K,(x) were taken from IsHERWOOD’S 
table (1904). 

T. V gives Lommel’s 6D table (1886) of J,(x), nm = — 63(1) + 64, x = 1(1)50, and 
nm = 7$(1)183, x = 1(1)20, and FRESNEL’s integrals, x = .5(.5)50, and x = .02(.02)1 with 
some modifications and corrections, and with the first 16 maxima and minima of the integrals. 
There are about 30 errors here. 

T. VI gives the 7D values of J,,(m), —Yn(m), Jn'(n), Yn'(nm) and niJ,(n), —n'Y,(n), 
niJ,'(n), ni Vn'(n) for n = 1(1)50. 

T. VII gives the first 40 zeros of J,(x) and Y,(x) to 7D, for n = 0(1)5, 4; and also of 
J—4(x) + J3(x), and of J—4(x) — J3(x). 

T. VIII gives the 7D values of 4/¢7Jo(t)dt, 4/7 Yo(t)dt for x = .02(.02)1, with the first 
16 maxima and minima of the integrals, to 7D. 

7. ‘The sum of series of cosecants,”’ Phil. Mag., s. 6, v. 45, 1923, p. 577-581. Compare 

n-1 
no. 2. S, = = csc? (mm/n), table for nm = [2(1)30(5)100, 360, 1000; 5D]. 
m=1 

8. ‘‘Theorems stated by Ramanujan (V): approximations connected with e?,’’ London 
Math. Soc., Proc., s. 2, v. 29, 1929, p. 293-308. Solutions of ue!~“* = Ue!-Y = e~ with values 

135 U u 
of ¢(t) = 3 {—~, - (1 =a \, tables of u, U, o(t), —t“ In (8), for ¢ = [0(.02)- 
1(.1)11; 7D]. 

9. “The constants of Landau and Lebesgue,” Quart. Jn. Math., Oxford s., v. 1, 1930, 
p. 310-318. Contains certain numerical results. 

10. ‘“‘The use of series of Bessel functions in problems connected with cylindrical wind- 
tunnels,’’ R. Soc. London, Proc., v. 130A, 1930, p. 29-37. A few numerical results. 

11. “Du Bois-Reymond’s constants,”’ Quart. Jn. Math., Oxford s., v. 4, 1933, p. 140-146. 
Various computed results are given. 

12. ‘Two tables of partitions,” London Math. Soc., Proc., s. 2, v. 42, 1937, p. 550-556. 
Tables of the number of partitions of m into odd parts and into distinct odd parts for 
n = 0(1)400. 

13. “‘Tabellazione di una particolare funzione definita da un integrale improprio,” 
R. Accad. naz. Lincei, Cl. d. sci. fis., matem. e nal., Rendiconti, s. 6, v. 27, 1938, p. 525-528. 
Table of o(a) = A*Lx + $ — (x? + x)!] sin Axdx, \ = [0(.01).4(.1)32; 7D]. 

14. ‘Three triple integrals,” Quart. Jn. Math., Oxford s., v. 10, 1939, p. 266-276. 
Certain numerical results. 

15. “A table of Ramanujan’s function r(m),’’ London Math. Soc., Proc., s. 2, v. 51, 
1949, p. 1-13. See MTAC, v. 3, p. 468. Gives a table of r(m) for m = 1(1)1000; also r(m)n-"/2 
to 5D. 


We are indeed grateful to Professor BROADBENT, the editor of the 
Mathematical Gazette, for consenting to allow us to reproduce the portrait of 
Professor Watson published as a frontispiece to the Gazette, v. 18, 1934. 


THE EDpitTors 
107. THE LEMNISCATE CONSTANT.—By this constant @ is usually meant 
twice the number 
& = A = 2(2)-{T)}? = BBG, D) = FQ, 21) 
—i\)-1 1 
{( | = 9r02(0|¢) = 1(2)M(—3)/M(—2) = f (1 — x‘)~! dx. 
4 0 


Nie 


It derives its name from the fact that it plays the réle of z in the rectifica- 
tion of the lemniscate. In fact the length of the lemniscate of diameter D is 
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#D. This constant occurs in a good many connections in the theory of elliptic 
functions and in number theory. The values 


(1) A = 1.311031, 
(2) A = 1.31102 87771 46059 87, 
(3) A = 1.31102 87771 46059 90680 3207, 


were given by EULER, STIRLING! and Gauss? respectively. The value (3) 
is given in the FMR Index. WreNcH [MTAC, v. 3, p. 202] has pointed out 
that this value is incorrect beyond the 17th decimal. 

In order to verify certain identities involving this constant with a number 
of numerical functions the writer has made two independent redetermina- 
tions of A using the two lacunary series 


io} wo 
A = dr{ zu (—1)"e-*"*}? =_ 2ime—i{ pit is asian 2 
the first of which was used by ‘Gauss. These calculations agree to 52D and 
give 
A = 1.31102 87771 46059 9% 


523 24 
94945 55970 68413 77475 71 


197 
1581. 
Gauss gave also the constant 
1 
B = 3x/o = 3in/A = (2n)-3{T(3)}? = f x2(1 — x*)-Idx 

correctly to 20D. More accurately we have 

B = .59907 01173 67796 10371 99612 

46140 16193 91136 06331 60783. 


As a byproduct of these calculations we obtain easily the following values of 
r(4) and r(2): 


(2) = 3.62560 99082 21908 31193 06851 
55867 67200 29951 67682 88007 
1(3) = 1.22541 67024 65177 64512 90983 


)3362 89052 68512 39248 10807. 
: oe. &. 
1 J. StiRL1NG, Methodus Differentialis, London, 1730, p. 58; English transl. by F. HoLtt- 


pay, London, 1749, p. 51. 
2K. F. Gauss, Werke, v. 3, Gottingen, 1866, p. 413. 


108. NEWMAN AND GEORGE ELIot.—We have made various references 
to the table-maker, FRANCIS WILLIAM NEWMAN (1805-1897), brother of the 
Cardinal; see for example MTAC, v. 1, p. 454-459, v. 3, p. 201, 257, 451. 
At the first of these references we noted (p. 455) that he was professor of 
Latin at University College, London, 1846-1869. In the Times Literary 
Supplement, London, v. 48, April 30, 1949, p. 281, there is an interesting 
letter by KATHLEEN TILLOTSON! noting, what had escaped her biographers, 
that GEorGE ELroT (1819-1880) was one of the earliest students of the Ladies’ 
College at Bedford Square, London, which opened in October, 1849 and is 
now Bedford College for Women. During 1851 at least, writes Mrs. Tillotson, 
she took courses there in History, German, French, Elocution, and Mathe- 
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matics. Her mathematics, and possibly her history, was with Newman, who 
taught mathematics at the College January—June 1851. In a letter of March 
27, 1874 George Eliot remembered her interest in Francis Newman's books 
“in far-off days . . . and . . . the awe | had of him as a lecturer in Mathe- 
matics at the ‘Ladies’ College.’’’ Mrs. Tillotson continues: ‘‘He had lectured 
in Ancient History since the opening of the college, and it is natural to 
suppose that it was partly his name . . . that drew MAry ANN EvVaANs to 
the Ladies’ College.” 

In the Times Literary Supplement, June 3, 1949, p. 365, Prof. GoRDON S. 
Haicut of Yale University quotes from a letter “‘which has just come to 
light” written by George Eliot Jan. 28, 1851. A passage quoted is as follows: 
“Il am attending Professor Newman’s course of lectures on Geometry at the 
Ladies’ College every Monday and Thursday. You will say that I can’t 
afford this, which is ‘dreadfully true’—but the fact is that I happened to 
say I should like to do so and good-natured Mr. CHAPMAN? went straightway 
and bought me a ticket which he begged me to accept. | refused to accept it— 
and have paid for it—wherefore I must stint myself in some direction— 
clearly in white gloves and probably in clean collars.”’ 





R. C. A, 


1 Mrs. TILLOTSON is a coeditor of the Works of Michael Drayton, and wife of GEOFFREY 
TiLLotson, professor of English Literature in the University of London. 

2 JoHN CHAPMAN, the publisher, at whose home George Eliot resided Jan. 8 to March 
24, 1851. 


109. RHETICUS, WITH SPECIAL REFERENCE TO HIS Opus PALATINUM.— 
Before the twentieth century Rheticus was the greatest calculator of mathe- 
matical tables who ever lived. He was born just a century before Napier 
published his work on logarithms in 1614, and his tables were posthumously 
published in 1596 and 1613. 

Since the ordinary inquirer would in more than one direction find con- 
siderable difficulty in deriving a fairly complete outline of the life and work 
of Rheticus, such an outline is here presented, together with a somewhat 
extended list of sources (everything has been personally inspected) from 
which fuller information may be derived. In this outline some errors in these 
sources have been corrected, a few new facts delineated, and our previous 
sketch of Pitiscus (MTAC, v. 3, p. 390-397) amplified. 

Rheticus was one of many mathematicians named after the places of 
their birth. GEORG JOACHIM was born 16 Feb. 1514 at Feldkirchen in Vorarl- 
berg, in the westernmost part of Austria-Hungary, close to the Swiss border 
and formerly in the Roman province of Rhaetia. Georg’s mother seems to 
have been connected with an Italian noble family of means, named Porrt,?° 
and thus a trip to Italy with his parents as a child was what might be ex- 
pected, and suggested a family in comfortable circumstances. Joachim’s 
zest for advanced education seems to have been stimulated by studies at 
Ziirich, Switzerland. Hence in 1532 the fame of the comparatively new 
University of Wittenberg, where the great reformer PHILip MELANCHTHON 
(1497-1560) was professor of Greek, and MARTIN LUTHER (1483-1546) 
was professor of Theology, drew the eighteen-year old youth to register as 
a student. In April 1536 he received the equivalent of a doctor of philosophy 
degree, the official records stating that he was ‘‘mathematicus excellens.” 
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It was exactly at this time that he first added the surname Rheticus. (It 
may be noted that in an outstanding German reference work" there are two 
biographies of Rheticus, by different authors, one under “Joachim” and 
the other under ‘‘Rheticus.’’) 

In May 1536 Rheticus went to Nuremberg to study with JOHANN 
SCHONER (1477-1547) professor of mathematics at the gymnasium there, 
and a former pupil of REGIOMONTANUS and WERNER.” Later he went to 
Tiibingen where he received the call back to his University.’ The professor 
of mathematics at Wittenberg having died in 1536 Melanchthon, then being 
Dean of the faculty, appointed for that autumn two new professors: one 
ERASMUS REINHOLD (1511-1553) for the higher mathematics (including 
astronomy), and the other, Dr. Rheticus, then 22 years of age, for elemen- 
tary mathematics. After teaching there with some repute for about three 
years Rheticus received leave of absence in May 1539 and went to Frauen- 
burg, Prussia, as an assistant to COPERNICUS (1473-1543), whose doctrines 
he soon advocated with much zeal and personal risk. 

Copernicus was very reluctant to publish his epoch-making astronomical 
work, De revolutionibus orbium caelestium, but publicity achieved for its 
contents by Rheticus led him to change his mind and to permit publication. 
This publicity was a survey of the principal features of the new astronomy, 
cast in the shape of a letter to his former teacher SCHONER at Nuremberg. 
This was published anonymously at Dantzig in 1540, as a Narratio prima, 
and reprinted at Basel in 1541 with the name of Rheticus as author. Of this 
Narratio there were nine later editions, the eleventh being that of ROSEN 
in English.?* Rheticus left Frauenburg at the end of September 1541, re- 
sumed his teaching in Wittenberg, and served as dean of the arts faculty 
there in the early months of 1542.*5 

Chapters 12, 13, 14 of the first of the six books of the De revolutionibus of 
Copernicus were devoted to the trigonometrical material, useful for the 
later discussion of astronomical problems. In the 12th chapter was a 5D 
table, with differences, of the sines of all angles in the quadrant at interval 
10’. In 1542 Rheticus had chapters 12 and 13 printed at Wittenberg with 
the table replaced by his own 7D table of sines and cosines for every minute 
in the quadrant. 

Rheticus left Wittenberg in 1542 and went to Nuremberg where the 
great work of Copernicus was being printed ; the early part was set up under 
his direction. Rheticus placed the first completed copy of the work in the 
hands of Copernicus a few hours before he died in May 1543. The second 
edition of De revolutionibus, edited by Rheticus appeared in 1566, with the 
appended third edition of his Narratio. 

In securing the publication of De revolutionibus Rheticus bears to Coper- 
nicus a relation similar to that of HALLEY to NEwrTOoN, in the following 
century, for assuring the publication of the Principia. 

The following statement of Miiller?* has been frequently quoted: 
“Without Rheticus we would know no Copernicus, without Copernicus no 
Kepler, without Kepler no Newton!” 

At Nuremberg in 1542 Rheticus also published Orationes de astronomia, 
geographia et physica. He taught at the University of Leipzig from 1542 to 
1551. In this latter year was published his Canon doctrinae triangulorum.'®. 
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Leipzig, 24 p., to which we have earlier referred,*® p. 394. In this quarto- 
format publication of 12 leaves, not listed in any of the bibliographies, 
p. 4-17 are occupied with a 7D Canon for all six trigonometric functions, at 
interval 10’, and this is followed by a dialogue between Philomathes, a 
supposed friend of Rheticus, and Hospes, his pupil. This was the first table 
(a) in which all trigonometric functions (or even sin, tan, sec) were brought 
together, or (b) in which the semiquadrantal arrangement was used. Both 
this canon and the undated reprint at Basle (apparently in 1580), are in the 
British Museum. With Rheticus our names for the trigonometric functions 
were not used; sines and cosines were perpendiculars and bases to a hypot- 
enuse 10 000 000; the secants and tangents were hypotenuses and perpen- 
diculars to a base 10 000 000; and the cosecants and cotangents were hypot- 
enuses and bases to a perpendicular 10 000 000. Rheticus was the first to 
define his functions by means of the right-angled triangle without any refer- 
ence to the circle. 

In 1551 Rheticus went to Prague; and in 1554 he lectured at the Uni- 
versity of Vienna,”‘ but in 1557 he settled in Cracow, where he spent most 
of the rest of his life. He died during a visit to Cassau, in Hungary, on 4 Dec. 
1576 in the sixty third year of his age.*! The great scarcity of the works of 
Rheticus (the second edition of De revolutionibus is much rarer than the 
first) is due to the fact that in 1550 all of his works were put on the Index 
Expurgatorius. 

During the last 25 years of his life Rheticus devoted himself principally 
to the calculation of monumental trigonometric tables. Contributions of 
financial assistance from MAXIMILIAN II of Austria (who died in 1576) and 
several Hungarian nobles enabled Rheticus to procure a number of com- 
puters for 12 years (3 and WoLF"). 


The Opus Palatinum. During this period of intense activity Rheticus 
had long known of a young mathematician, a native of Magdeburg, Lucius 
VALENTIN OTHO (1550?-1605?) well known!’ to PRATOoRIUS at the Univer- 
sity of Wittenberg. Towards 1575 Rheticus gladly accepted the services as 
assistant which Otho offered. And thus it turned out that when Rheticus 
died in December 1576 Otho became heir to all of the scientific mss. of 
Rheticus. Rheticus has told us that when Otho came to him he was the same 
age (25) as when he (Rheticus) became an assistant to Copernicus (‘‘Pro- 
fecto in eadem aetate ad me venis, qua ego ad Copernicum veni’’). In this 
way we learn that Otho was born about 1550. He became a professor of 
mathematics (following Pratorius) at the University of Wittenberg, and was 
continually seeking to procure the means for publishing the great tables 
which were practically complete. This was finally arranged with FREDERICK 
IV,*! Elector Palatine, to whom the work is dedicated. Thus 20 years after 
the death of Rheticus, in 1596, his Opus Palatinum de triangvlis a Georgio 
Ioachimo Rhetico coeptum: L. Valentinvs Otho, Principis Palatini Friderici IV. 
Electoris Mathematicvs consummavit. was published at Neustadt, now in 
Bavaria, in about 1446 folio pages, 22.4 X 36.7 cm. It seems to have been 
originally bound in vellum in two volumes. Five of the seven parts were in 
the first volume. 

After the title-page (which we here reproduce) and 9 preliminary leaves 
of dedication and preface by Otho are 7 main sections each (after the 
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second) separately paged, and for each of the first 5 sections there is at 
least one separately dated title-page. The first three of these sections, by 
Rheticus, are as follows: 


1. De Fabrica canonis docirinae triangulorvm. The title page, with verso blank, is followed . 
by p. 3-86. P. 45 is printed and folded so as to occupy the verso of a blank leaf and the recto 
of p. 46, making an extra leaf. Pages 34, 35, 37 and 38 are incorrectiy numbered, while the 
verso of p. 85 also bears the number 85 instead of the correct number 86. Here are three books 
on the construction of the canon. 

II. Trigvetrum rectarvm linearvm in planitie cum angvlo recto, p. 86-104. One book on 
plane triangles. 

III. De triangzlis globi cum angvlo recto. 146 p. Four books on right-angled spherical 
triangles. The verso of p. 21 is 23, the recto of p. 23 is 24, and the verso of p. 24 is blank. 
P. 57 is numbered 60. P. 72-88 are numbered 81, 83, 82, 85-98. Then p. 89-106 are respec- 
tively numbered: 100-103, 95, 105, 106, 98, 108-115, 101, 102. Here a blank leaf is inserted 
which we shall call p. 107-108. Then corresponding to p. 109-146 are pages numbered 103- 
140. In addition there is inserted between p. 115(109)—116(110) a printed half-page contain- 
ing matter left out of the text in the printing. 

Following this material by Rheticus are 5 books by Otho on oblique spherical triangles 
and 3 subsidiary astronomical tables which Otho calls meteoroscopia. 

IV. De triangvlis globi sine angvlo recto libri qvingqve. The title page with verso blank is 
followed by p. 1-264, then commencing on the next recto the pages run from 264 (repeated) 
to 341 on a verso. The recto of the next leaf contains the following: ‘‘Neostadii in Palatinatv. 
Excudebat Matthaeus Harnisius [printer’s device ]. Anno Salutis MD XCVI.” The verso 
of the last mentioned leaf is blank and it is followed by a blank leaf. 

V. Meleoroscopivm numerorvm primum. p. 1-91. Meteoroscotivm numerorvm secundym. 
p. 92-105. Meieoroscopivm numerorvm tertivm. p. 106-121. Each of these three parts has its 
own dated title page. Then follows a closely packed recto page entitled ‘‘Errata sic corrige,”’ 
corrections to be made in this volume I. 

Volume II is devoted to tables, each with an undated title page. The first of these is the 
great Rheticus 10D tables at interval 10”, for all six of the trigonometric functions. 

VI. Georgii Ioachimi Rhaetici Magnus canon doctrinae triangvlorvm ad decades secvn- 
dorvm et ad partes 50000000000. All of this is printed in red, after the name of Rheticus in 
black, 554 p. The table occupies p. 2-541. Then come Errata p. 542-554, incorrectly num- 
bered 142-147, 147 repeated, 549-554. P. 210 reads 201, 290 reads 390, 309 reads 314, 398 
reads 498, 512 reads 125. 

Following these errata is a 7D table at interval 10”, of cotangents and cosecants for 
the first half of the quadrant: 

VII. Tertia series magni canonis doctrinae triangulorym in qvo trigqvetri cum angvlo recto in 
planitie minvs latvs includentivm angulum rectum ponitur partium 10000000. After the first 
two words in black, the rest of the title is red. 181 p. P. 83 is incorrectly numbered 27. 


Concerning this section VII DeMorgan remarks!" that its insertion dis- 
played ‘‘merely the editor’s want of judgment; it is clearly nothing but a 
previous attempt made before the larger plan was resolved on, and is much 
less accurate than the great table to ten places.”’ 

DELAMBRE,’ p. 2, states that in II!I—-1V Rheticus and his disciple are the 
most prolix and obscure authors he had ever met; and that their 500 folio 
pages could be reduced to 10. 


Pitiscus Revision of the Opus Palatinum Canon. Shortly after the 
Opus Palatinum was published it was found that in VI the tangents and 
secants near the end of the quadrant were very inaccurate. Possibly through 
representations of Pitiscus, then chaplain and former teacher of Frederick 
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IV, the Elector Palatine put Pitiscus in charge of the revision of this part 
of the canon. From the account which Rheticus gave of the construction of 
the canon in I, especially book 3, it is clear that the computation of the 
tangents and secants were based on a table of sines and cosines, at interval 
10’’, and extent > 10D. Hence Pitiscus, perhaps about 1602 (when he was 41 
years old), sought the ms. of Otho “then an old man” (say La Lande, 
Delambre,® p. 17, DeMorgan," Glaisher'*—based on the statement of 
Pitiscus: “ob memoriae senilis debilitatem’’—who could not remember 
where the ms. was). He thought that perhaps he had left it at Wittenberg, 
and accordingly Pitiscus sent a messenger there to search for it; but after 
considerable expense had been incurred he returned without it. In 1602 
Otho was about 52 years old, just 11 years older than Pitiscus. While the 
condition of his health may have affected his memory it is certainly im- 
proper*! to refer to Otho at this time as an “old man.” 

After the death of Otho, at Heidelberg, about 1605, when the mss. of 
Rheticus which had been in his possession passed into the hands of the 
orientalist and astronomer at Heidelberg, JAcop CHRISTMANN, he found 
that they included not only the sine canon which had been considered lost, 
but also the original manuscript of the De revolutionibus of Copernicus, now 
preserved in Prague. 

As soon as Pitiscus learned of this canon discovery he examined the mss. 
page for page, although they were in very bad condition, and to his great 
satisfaction he found 


(1) the ten-second canon of sines, 15D, A‘; 
(2) sines for every second of the first and last degrees of the quadrant, 15D, A*; 


(3) the commencement of a canon for every ten seconds of tangents and secants, 
15D, A?; 
(4) a complete minute-canon of sines, tangents, and secants, 15D. 


From this list of tables, taken in conjunction with VI and VII, we are en- 
abled to view as a whole the colossal sixteenth century computations which 
Rheticus achieved, and which must later have been of extraordinary value 
to the astronomer. 

With the canon (1) in hand Pitiscus recomputed to 11D all of the tan- 
gents and secants of VI in the defective region from 83° to the end of the 
quadrant. Then 86 pages were reprinted: p. 1, the new title-page, and 2-86, 
the revised table. For the new title-page to that already given in VI was 
added, 1. 

Recens emendatus a4 Bartholomaeo Pitisco Silesio. Addita est brevis com- 
monefactio de fabrica & vsu huius Canonis. Quae est summa doctrinae, & 
quasi nucleus, totius operis Palatini. Canon hic, und cum brevi commonefac- 
tione de eius fabrica & vsu etiam separatim ab opere Palatino venditur. In 
Bibliopoleio Harnischiano. 

This is a sort of fly-title, undated ; but the accompanying 19-page publi- 
cation has the following title page, 2: Bartholomaei Pitisci Griinbergensis 
Silesii Brevis et Perspicua Commonefactio de Fabrica et vsv magni canonis 
doctrinae Triangulorum Georgii Ioachimi Rhetici. [Printer’s decoration ]. 
Neostadii, Typis Nicolai Schrammii. MDCVII. 

The 86 pages of 1 and the 19 pages of 2 were issued, as stated above, 
as a separate publication, but also with copies of the complete work with 
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the original pages 1-86 in VI eliminated and replaced by 1, and with 2 
following the revised VI, or VII.* It is the latter arrangement which is found 
in the only three existing copies of 1-2, which I have been able to trace, 
and to which I shall presently refer. 


Otho’s Death. Since on page 5 of 2 there seems to be definite reference 
to the Rheticus sine canon being in the hands of Christmann, this suggests 
that Otho had died before 1607 and hence I suggested that he may have 
died about 1605 and that Pitiscus may have questioned him concerning the 
table about 1602. While I arrived at this conclusion as to a possible date of 
Otho’s death independently I was later interested to find that TROPFKE had 
reached the same result in the second ed. of his Geschichte der Elementar- 
Mathematik, v. 4, 1923. The only definite older reference to the death of 
Otho is in the ‘“‘Ad lectorem”’ of the work we shall now consider, namely : 


3. The Rheticus-Pitiscus Thesaurus Mathematicus. After his discovery 
of the new Rheticus tables Pitiscus started to prepare a second work which 
was finally published in 1613 and contained the following four parts. 


A. (Rheticus) Canon of sines for every 10’, 15D, A?, p. 2-271; 

B. (Rheticus) Sines for 0(1’)1°, 89°(1’")90°, 15D, A?, p. 2-61; 

C. (Pitiscus) The fundamental series from which the rest were calculated to 22D, 
p. 1-10; 

D. (Pitiscus) The sines to 22D, A‘ or A‘, for every tenth, thirtieth and fiftieth second in 
the first 35 minutes, p. 11-15. 


Each of the parts C, D has its own correctly dated title-page 1613; but the 
title-page for B as well as that for the whole volume is given incorrectly as 
MDXIII. The volume title-page, in black and red, begins Thesaurus Mathe- 
maticus sive Canon Sinuum ad radium 1.00000.00000.00000. . . . Frank- 
furt. The volume is a large one 24X35 cm. It appeared shortly before the 
death of Pitiscus 2 July 1613 in his fifty second year. From what has been 
presented above the great contribution which Pitiscus made in uncovering, 
and presenting work of Rheticus in meticulously correct form, was indeed 
a notable contribution to science. 

More information concerning this volume is given in MTAC,* p. 395- 
396. 

One had to wait for nearly three and a half centuries before the tables 
of Rheticus were finally superseded by those of ANDOYER (1915-1918). 


Known copies of the 1607 Rheticus-Pitiscus work. 


The three known copies of this work are in Edinburgh, at the Crawford 
Library of the Royal Observatory, which we shall call the ECO copy; in 
Washington, at the Naval Observatory, the WNO copy; and in the private 
library of Mr. Witt1aAm D. Moraan of St. Paul, Minnesota—the StPM 
copy. 

For the preparation of this article I have been deeply indebted to Mr. 
Morgan for allowing me to have his StPM copy for leisurely study. In at 
least one respect it is the finest copy in a single library. It consists of three 
volumes, purchased at different times from three different firms: (i) almost 
certainly the volume of I-V, bound in contemporary white vellum with 
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leather tiers as originally issued in 1596, 36.7 cm. tall, and thus indicating 
that ai least some copies of the original were bound in two v. and divided 
in that way; (ii) containing V—VII, as published in 1596, 1.7 cm. less tall 
than (i); (iii) containing 1 along with the remaining 454 pages of the trigo- 
nometric canon, 2, and 3. This third v. is 3 cm. less tall than (i), and at least 
1.3 cm. less tall than 3 when it was originally published. In StPM the many 
fundamental changes in VI of Rheticus which Pitiscus made in 1 may be 
readily observed. In this respect and in the superb condition of (i) StPM is 
outstanding. 

As noted above, in both ECO and WNO there is replacement of the first 
86 pages of VI by 1. The great kindness of Mrs. Grace O. SAavaGeE, the 
librarian at the Naval Observatory, in responding to my numerous queries 
concerning WNO, has been much appreciated. This copy, as tall as (i) and 
with wider pages, is also bound in vellum but with different ordering of 
material. The Library of Congress Rheticus lacking the Pitiscus revision, 
36.5 cm. tall, is practically as tall as (i). Brown University has film copies 
of the Library of Congress Rheticus, and of Mr. Morgan’s 1607 Pitiscus 
items. 


We have earlier (MTAC, v. 3, p. 396) indicated the reason for the 
scholar’s great interest in the ECO copy. 


R. C. A. 


1 ADRIAAN VAN ROOMEN (or ADRIANUS RoMANUs), Ideae mathematicae pars prima, sive 
methodus polygonorum. | Louvain and mes 4 1593. In preliminary material Roomen has 
(i) “‘lectori philomathi”’ in which he has words of praise for 15 of the principal mathematicians 
of his time (including Rheticus and VALENTIN OrHo), and (ii) a letter written by Rheticus 
in 1568 to Petrus Ramus (1515-1572), telling not only of his great tabular calculations but 
also of such works as the following which he planned to publish: (a) Observation of Phe- 
nomena; (b) German astronomy; (c) Natural Philosophy; (d) Foundations of Chemistry. 
See H. Bosmans, “Romain (Adrien),” Biogr. Nationale . . . de Belgique. Brussels, v. 19, 
1907, col. 855; also A. MOLLER,” p. 40. 

J. C. F. pe LaLanpe, “Lettre sur tables de sinus extrémement rares,” Jn. des 
Scavans, Paris, Sept., 1771, p. 579-584. The Rheticus-Pitiscus Thesavrvs Mathematicvs, 1613. 
2A. =. 3. KASTNER, Geschichte der Mathematik. Géttingen, v. 1, 1796, p. 561f, 590f; em 


1797, 

GE He. Hutton, Mathematical Tables: containing Common, Hyperbolic, and Logistic Loga- 
rithms. . To which is prefixed a large and original history of the discoveries and writings 
relating to those subjects. . . . London, 1785, p. 9-11; also in his Tracts on Mathematical 
and fey on Subjects. V. :. London, 1812, p. 290-293 

5 JEAN BERNOUELI, “Analyse de I’ Opus palatinum de _ & du Thesaurus mathe- 
maticus de Pitiscus,”’ Akad. d. Wissen., Berlin, Nouveaux Mémoires, 1786, p. 10-33. 

B. DELAMBRE, “Rapport sur les grandes tables trigonométriques décimales du 
cadastre," Acad. d. Sci. Paris, Mémoires, Sci. Math. et Phys., v. 5, 1804, p. 56-66. 

"4's RICHE DE PRoNy, “Eclaircissemens sur un point de l'histoire des tables 

trigonométriques,” Acad. d. Sci., Mémoires, Sci. Math. et Phys., v. 5, 1804, p. 67-93. 
DELAMBRE, Histoire de I’ Astronomie Moderne, v. 2, Paris, 1820, p. 1-127: 
“Grandes tables trigonométriques en nombres naturels—Rheticus.” 

® A. DEMorGav, “Rheticus,” Penny Cycl., v. 19, 1841, also The English Cycl., Biography, 
v. 5, London, 1857. 

10 A. DEMorGan, “On the almost total disappearance of the earliest trigonometrical 
canon,” RAS, Mo. Not., v. 6, 1845, p. 221-228; reprinted with an addition in Phil. Mag., s. 3, 
v. 26, 1845, p. 517-526. 

1 A, Deabtoncam, (i) “Table” in Su x to the Penny Cycl., v. 2, London, 1846, p. 
pon pak 599-600. (ii) “‘Table,”” English Cycl., Arts and Sci. Sect., v. 7, 1861, cols. 984, 

12 J. C. PoGGenporrF, Biogr.-Literar. Handwérterbuch, v. 2, Leipzig, 1863. 

ht , L. GLAISHER, Report of the Committee on Mathematical Tables, London, 1873, 
p. 4 b 


“F. Hipcer, “Die chorographie des Joachim Rheticus. Aus dem Autographon des 
Verfassers mit einer Einleitung herausgegeben,” Z. Math. Phys., Hist.-literar. Abth., v. 21, 
1876, p. 125-150. 
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15 R. Wor, Geschichte der Astronomie. Munich, 1877, p. 204f, 236f, 242f, 296, 343f. 

16 Ailgemeine deutsche Biographie, Leipzig, v. 14, 1881, p. 93-94: “Joachim, G.” by K. C. 
BRUANS; V. 28, 1889, p. 388-390; ‘‘Rheticus, G. J.” by S. GUNTHER. 

17 LEOPOLD PROWE, Nicolaus Co pernicus. Berlin, 2 v., 1883-1884; v. 1, Das Leben, v. 2, 
Urkunden. “For the biography of Copernicus and the social history of his times the book 
- oe valuable, but Prowe’s judgment in scientific matters was unreliable’ (E. 

OSEN ). 

18M. Curtze, 1. “Zur Biographie des Rheticus,”’ Altpreussische Monatsschrift, v. 31, 
1894, p. 491-496. This article contains an extract from a Munich library ms., Codex latinus 
Monacensis no. 24101, written by JOHANNES PRATORIUsS (1537-1616), inventor of the plane 
table used in surveying, who was professor of mathematics at the University of Wittenberg 
1571-1576, and professor at the University of Altdorf from 1576 until his death. In this ms. 
we are told that in Aug. 1573 Otho, then mathematician for the landgrave of Hesse, came 
to his patron Pratorius with two notable approximations for 7: 

(a) 3.1415926537 > x > 3.14159265365, published by Viera 20 years later, in 1593; 

(b) 355/113, which is correct to 6D. 


The latter he had derived from the approximations 377/120, found by ProLemy, and 22/7 
given by ARCHIMEDES, by subtractions of numerators and denominators. This is the earliest 
sure date for discovery of this approximation. This statement is made with full knowledge 
of the attribution of this result to: (a) the fifth Century Chinese astronomer Tsu Cn’UNG- 
Cutn by Y. Mikami, The Developmeni of Mathematics in China and Japan (Abh. z. Gesch. 
d. math. Wissen., Heft 30). Leipzig, 1913, p. 50; and (b) ADRIAEN ANTHONISz (1527-1607) 
— 1583 by his son ADRIAEN METIUs in his Arithmeticae libri duo et geometriae VI, 1626, 
p. 51. 
2. “Die abgekiirzte Multiplication,” Z. Math. Phys., hist.-literar. Abt., v. 40, 1895, 
p. 7-13. Abridged multiplication and mss. of RHETICUS are here discussed. The names of 
Orno, Piriscus and WERNER, among others, also come up. 

19K. Hunratu, “Des Rheticus Canon doctrinae triangulorum und Vieta’s canon 
mathematicus,”’ Abh. z. Gesch. d. Math., v. 9, 1899, p. 211-240. 

20M. Cantor, Vorlesungen tiber Geschichte der Mathematik. V. 2, second ed., Leipzig, 
1900, p. 472-475, 600-603, etc. 

21A. VON BRAUNMUHL, Vorlesungen tiber Geschichte der Trigonometrie. Leipzig, v. 1, 
1901, p. 140-149, 158-160, 212-226. (Das Opus Palatinum and Pitiscus), etc.; v. 2, 1903, 7 
index references. 

#2 JOHANN WERNER (1468-1528) of Nuremberg, astronomer and meteorologist, De 
triangulis sphaericis libri quatuor. Edited by A. A. BJORNBO; De Meteoroscopiis libri sex, 
cum prooemio G. J. Rhetici. Cracow 1557. Edited by J. Wirscumipt. Abh. 2. Gesch. d. 
Math. Wissen., v. 24, 1907, 1913. 

23 ApOLF MULLER, ‘Der Astronom und Mathematiker Georg Joachim Rheticus,”’ 
Vierteljahrsschrift fiir Geschichte und Landeskunde Vorarlbergs, n. s., v. 2, 1918, p. 5-46. 
Valuable collection and study of material about the life of Rheticus. 

% “‘Vorarlberger an in- und auslaindischen Hochschulen vom Ausgange des XIII. bis 
zur Mitte des XVII. Jahrhunderts,’’ Forschungen zur Geschichte Vorarlbergs und Liechten- 
steins, v. 1, 1920, p. 58, 128-130. 

26 Franz HAEFELE, “Zur Frage der Herkunft des Astronomen Georg Joachim de Porris,”’ 
Schriften des Vereines fiir Geschichte des Bodensees und seiner Umgebung, Heft 55, Fried- 
richshafen a.B., 1927, p. 122-137. Some new material. 

26 MARTIN BILGERI, Das Vorarlberger Schrifttum, und der Anteil des Landes am deutschen 
Geisiesleben. Vienna and Leipzig, 1936. ‘“‘Georg Joachim, genannt Rheticus, der grosse 
Astronom und Mathematiker,” p. 64-70. Nothing original. 

27 JOHANNES KEPLER, Gesammelte Werke, Rerauagegeben me Max Caspar, v. 1, 
Munich, 1938; Rheticus, Narratio prima, tenth edition, p. 88-13 

28 EDWARD ROSEN, Three Copernican Treatises: The tan of Copernicus, The 
Letter against Werner, The Narratio prima of Rheticus, translated with Introduction and Notes. 
New York, Columbia Univ. Press, 1939. Admirable work. Introduction: ‘Georg Joachim 
Rheticus,” p. 4-6; ‘The Narratio Prima,” p. 9-11. Text of “‘Narratio Prima,” p. 107-196, 
eleventh ed. 

29 ErnsT ZINNER, Entstehung und Ausbreitung der Coppernicanischen Lehre. Zum 200- 
jahrigen Jubilium der Friedrich-Alexander-Universitat su Erlangen. Phys.-medizin. Soz. z. 
— Sitz., v. 74, 1943. Very many references to Rheticus in the index. 

~ C. ARCHIBALD, ‘ ‘Bartholomaus Pitiscus (1561-1613),” MTAC, v. 3, p. 390-397, 
498-499, 1949. 

4 It is rare indeed that DeMorgan anywhere ever deviates in the slightest degree from 
exactness of statement; yet in connection with Rheticus three cases of this kind are to be 
noted. In no. 10 p. 228 it is stated that Rheticus died in the “sixty-first year of his age” 
instead of sixty-third. In nos. 9 (1857) and 11 it is stated that the Opus Palatinum was 
published “‘at the expense of the Emperor Maximilian,” who had died 20 years earlier. About 
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1600 when Otho was about 50 years old DeMorgan refers to him in no. 11 as “then an old 


# As DEMORGAN pointed out, the 86 pages of 1 are easily distinguishable by the in- 
feriority of paper and type; the same is true of 2; also the corrected copies may be distin- 
guished from the uncorrected ones in a moment as follows: look at the bottom of page 7, 
at the running titles of the columns. The uncorrected copy will have as it ought to have 
Basis Differentia Hypothenusa. But the corrected copy will have, as it ought not to have, 
Hypothenusa Differentia Basis. 


QUERY 


32. FRENCH AND RUSSIAN TRANSLATIONS OF A VEGA-BREMIKER TABLE. 
In RMT 635, reference was made to the Carl Bremiker edition of Vega’s 
Logarithmisch-Trigonometrisches Handbuch, first published in 1856. It was 
also noted that in 1857 an English translation of this edition by W. L. F. 
FISCHER, and an Italian translation by Luici CREMONA, were published. 
Who were the authors of the French translation of 1857 and of the Russian 
translation of 1858 also mentioned? The first two translations, as well as the 
1857 (second Bremiker) Gernian edition are in the Library of Brown Uni- 
versity. In what library may the French or Russian translations be found? 

 <. A, 


QUERIES—REPLIES 


41. GIRARD AND SNELL TABLES (Q30, v. 3, p. 451).—ALBERT GIRARD, 
Tables des Sinus, Tangentes, & Secantes, selon le Raid de 100000 Parties. 
Avec un traicté succinct de la Trigonometrie tant des triangles plans, que 
spheriques. The Hague, 1626, 240 pages (unnumbered). 6.8 X 11.8 cm. 

An examination of the copy of this book in the Library of Congress 
revealed that the main table contains natural sines, tangents, and secants 
in units of 10-° for every sexagesimal minute of the first quadrant, so ar- 
ranged that the functions of complementary angles appear on facing pages. 
Each page contains functions for a range of half a degree. 

Following this principal table is a section devoted to the statement and 
illustration of rules for the solution of the four standard cases of oblique 
plane triangles. Although Girard frequently resorts to the device of dissect- 
ing oblique triangles into right triangles, he does state and use the Law of 
Sines. In addition, he gives in the form of rules both the Law of Tangents 
and the Law of Cosines,—the latter in a form involving the versed sine. 

A brief section dealing with some general theorems relating to plane 
polygons is followed by a treatment of both right and oblique spherical 
triangles. The discussion of the solution of oblique spherical triangles is 
limited to three cases: (i) all angles given; (ii) all sides given; and (iii) two 
sides and their included angle given. 

The author then gives to five significant figures the length of a side of 
each of the five regular polyhedra when inscribed in a sphere of diameter 
100 000. His results for the regular tetrahedron and regular octahedron 
contain rounding errors. 

A more extensive table is included showing the lengths to the nearest 
integer of the sides of regular n-gons [m = 3(1)24] inscribed in a circle of 
diameter 200 000. Careful examination showed this table to be entirely free 
from error. 
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After a brief section containing the statement of geometrical problems 
involving the solution of equations such as tan @ = 26, sin @ = 6/v2, 
6 — sin 0 = 2/2, tan 26 — tan @ = 1, in modern notation, Girard reproduces 
Ludolph van Ceulen’s 36-figure approximation to 7, and presents a tabula- 
tion of arc length to radius 10° for the following central angles: 1°(1°)10°- 
(5°)100°, 150°, 180°, 1’(1’)10’(5’)55’, 1’7(1")10(5")55”. This table con- 
tains ten rounding errors. 

Two rules are stated for the calculation of close approximations to an 
acute angle of a right triangle in terms of the three sides. In modern trigo- 
nometric notation these rules are equivalent to the inequalities: 


(1) 3 sin 6/(2 + cos 0) < @ < 3(2 sin @ + tan @). 


The book concludes with miscellaneous information relative to triangles, 
including a schematic array showing the manner of dependence of the area 
of a plane triangle on various geometric magnitudes associated with a 
triangle. 


Joun W. WRENCH, Jr. 
4711 Davenport St., N.W. 
Washington, D. C. 


EpiroriaL Note: The inquiring reader will find it of value to read the remarks concern- 
ing this work in A. VON BRAUNMUHL, Vorlesungen tiber Geschichte der Trigonometrie, v. 1, 
1900, p. 237-238. It may be remarked that Dr. Wrench’s description of Girard’s main table 
is exactly equivalent to that of Pitiscus, 1600 (MTAC, v. 3, p. 391). 

The problem leading to the equation tan @ = 26 is as follows (compare N104, LoTHaR 
Cotitatz, MTAC, v. 3, p. 496): What is the arc which equals half its tangent? 

A relation differing from (1) only in that the right member is 2 sin $@ + tan 30 was 
first given by WILLEBRORD SNELL, in his Cyclometricus, Leyden, 1621, prop. XXIX, p. 43 f. 
The first published proof of this right-hand member relation was given by HuyGENs in his 
De Circula Magnitudine inventa, 1654. If 3y = 6 the relation 6 < 2 sin 3@ + tan $0 becomes 
y < 4(2sin y + tan y), as given by Girard. The single approximation @ ~ 3 sin 0/(2 + cos 8), 
was first surmised by NicHotas De Cusa (1401-1464); see his Opera, 1514 or 1565, p. 1120- 
1154. Compare K. T. VAHLEN, Konstruktionen und Approximationen in systematischer 
Darstellung. Leipzig, 1911, p. 188-190. 


42. Pitiscus TABLEs (Q29, v. 3, p. 398).—In 1927 I purchased for $180, 
from Wheldon & Wesley, London, a copy of RHETIcus-Pitiscus, Thesaurus 
Mathematicus, Frankfurt, 1613, bound with the 1607 volume (86 p.) de- 
scribed in MTAC, v. 3, p. 395. I have also acquired (1924-37) a fine copy 
of the two volumes (6 parts) of the original 1596 RuETICUs-OTHO, Opus 
Palatinum de Triangulis. With reference to the quality of paper and printing 
of the 1607 Pitiscus item the inferiority is notably in evidence. 


WiLtiam D. MorGANn 
1764 St. Anthony Ave. 
St. Paul 4, Minnesota 


EpitoriaL Note: Among the many other valuable volumes of mathematical and 
astronomical tables in Mr. Morgan’s Library are the following: 1. Ptolemy’s Table of Chords, 
first published in the Liechtenstein edition of the Almagest, 1515; also the George of Trebi- 
zond edition of 1528, and the Greek edition of Symon Grynaeus, 1538. 2. Regiomontanus, 
Tabulae Directionum, 1552; also Wittenberg, 1606 edition. These volumes contain a table 
of sines at interval 1’ and calculated for a radius of 60 000. 3. Vieta, Canon Mathematicus. 
Paris, 1579, formerly the property of CHarLEs HutTrTon, purchased in 1930 from Sotheran 
for $200. The great scarcity of this work (of which there is a copy at Brown University) may 
possibly be due to the fact that on account of errors in the volume, Vieta bought up and 
destroyed all copies which it was possible for him to locate. 4. B. Pitiscus, Trigonometry, 
1612 edition (MTAC, v. 3, p. 391). 5. H. Briccs, Arithmetica Logarithmica. London, 1624 
(MTAC, v. 1, 8 97-98, 170). 6. A. VLacg, Arithmetica Logarithmica. 1628. 7. H. BricGs, 
Trigonometria Britannica. London, 1633. 8. H. SHERWIN, Mathematical Tables, first edition. 
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London, 1717. 9. J. Dopson, Anti-Logarithmic Canon. London, 1742. 10. G. von VEGa, 
Thesaurus Logarithmorum Completus. Leipzig, 1794. 11. Astronomical Tables of King 
Alphonse X of Castile (1203-1284), first printed edition, Venice, 1483; also editions of 1492, 
1518, 1524, and Paris, 1545, 1553. 12. E. Hattey, Tabulae Astronomicae. London, 1749. 13. 
J. Cassini, Tables Astronomiques. Paris, 1740. 14. Mr. Morgan is also the fortunate owner 
of what appears to be a fifteenth century vellum manuscript of the extraordinarily popular 
Toledo Planetary Tables, which originated with the astronomer Ibn al-Zarqala (1029?- 
1087?). This manuscript consists of 124 leaves measuring 17 X 25 cm. (See G. Sarton, In- 
troduction to the History of Science, v. 1, Washington, 1927, p. 758-759.) 


CORRIGENDA 
V. 1, p. 466, for Burden, J. A., read Burdon, J. A. 
V. 2, p. 384, for e see Napierian Base, read e 54-55, 454; II 15, 68-69, 300; p. 387, delete 


= 


p. 

2 * 

Napierian Base e 54-55, 454; II 15, 68-69; p. 389, for Weiertrass, read Weierstrass. 

V. 3, p. 385, for M. R. Hestenis, read M. R. Hestenes; p. 466, 1. 2, for (1.570780, read 
1.570780; p. 475, 1. —12, for Re(S), read Re(s); p. 496, 1. 2, for 10-5, read <10-*. 








Index to MTAC III: 21-28 
1948-1949 


The Indices for this volume bring together references to all Articles, RMT, MTE, UMT, 
ACM, OAC, N, Q, QR, Corrigenda, Subjects, and Names. Mrs. D. H. LEHMER has once 
more continued her labor of love by preparing the Subject Index, and for this we are deeply 
grateful. Our octogenarian friend, S.A.J., again has paid special attention to indexing 
Errata but his notable contribution here and in connection with the whole Index Section 
is but a tiny fraction of what he has done towards making the volume as a whole more 
worthy of publication. 


HEADINGS FOR INDICES 


Numbers issued Automatic Computing Machinery 
Articles Other Aids to Computation 
Recent Mathematical Tables Notes 
Corrigenda Portraits 
Mathematical Tables—Errata Queries 
Mathematical Tables—Errata, Author Queries—Replies 
Index Subject Index 
Unpublished Mathematical Tables Name Index 


NuMBERS ISSUED 


21, 1948, Jan. (1-68) ; 22, April (69-148) ; 23, July (149-228) ; 24, Oct. (229-332) ; 25, 1949, 
Jan. (333-398); 26, Apr. (399-452); 27, July (453-500) ; 28, Oct. (501-596). 


ARTICLES 


F. L. Alt, A Bell Telephone Laboratories’ computing machine, 1-13, 69-84 

J. Laderman, The square root method for solving simultaneous linear equations, 13-16 

H. E. Salzer, Coefficients for expressing the first twenty-four powers in terms of the Legendre 
polynominals, 16-18 

D. F. Ferguson & J. W. Wrench, Jr., A new approximation to x (conclusion), 18-19 

W. J. Eckert, The IBM pluggable sequence relay calculator, 149-161 

H. F. Mitchell, Jr., Inversion of a matrix of order 38, 161-166 

H. E. Salzer, Coefficients for expressing the first thirty powers in terms of the Hermite 
polynomials, 167-169 

A. Fletcher, Guide to tables of elliptic functions, 229-281 

A. V. Haeff, The memory tube and its application to electronic computation, 281-286 

R. M. Bloch, R. V. D. Campbell & M. Ellis, The logical design of the Raytheon com- 
puter, 286-295 

J. C. P. Miller, The BAASMTC now RSMTC, 333-340 

F. E. Snyder & H. M. Livingston, Coding of a Laplace boundary value problem for the 
UNIVAC, 341-350 

A. Hillman & I. Sherman, Complex zeros of Yo(z), Y¥i(z), and Y1'(z), 351-352 

G. W. King, A method of plotting on standard IBM equipment, 352-355 

R. A. Fairthorne & J. C. P. Miller, Hilbert’s double series theorem and principal latent 
roots of the resulting matrix, 399-400 

J. O. Harrison, Jr., Piecewise polynomial approximation for large-scale digital calculators, 
400-407 

H. E. Stelson, The accuracy of linear interpolation in tables of the mathematics of finance, 
408-412 
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F. M. Verzuh, The solution of simultaneous linear equations with the aid of the 602 calcu- 
lating punch, 453-462 


G. D. McCann, The California Institute of Technology electric analogue computer, 501-513 


RECENT MATHEMATICAL TABLES 
Reviews by numbers, authors, pages. There is no number 607 


450 (Mitchell) 19; 451 (Schulze) 19-20; 452 (Miiller & Rajna) 20; 453 (Proke3) 20-21; 
454 (Kulikov) 21; 455 (Wijngaarden) 21; 456-457 (Gloden) 21-22; 458 (Gupta) 22-23; 
459 (Lahiri) 23; 460 (Mahler) 23; 461 (Robert) 23; 462 (Walfisz) 24; 463 (Ziaud-Din, 
Kerawala, Kerawala & Hanafi) 24; 464-465 (NBSCL) 25-26; 466 (Reiz) 26; 467 (King) 
27-29; 468 (Freeman) 29; 469 (Hage) 29-30; 470 (Ebsen) 31; 471 (Germany, Deutsche 
Seewarte) 31-32; 472 (Germany, Reichsluftfahrtministerium) 32-34; 473 (U. S. Hydrogr. 
Off.) 34-36; 474 (Emmons) 36-37; 475 (MIT) 37-40; 476 (Barlow) 84; 477 (Cazzola) 
84-85; 478 (Gauss) 85-86; 479 (Comrie) 86-87; 480 (Levallois) 87; 481 (Gloden & Bonneau) 
87; 482 (Cracknell) 87; 483 (Rohrer) 88; 484 (NBSCL) 88-89; 485 (Spenceley & Spenceley) 
89-92; 486 (Anan’ev) 93; 487 (Marchisio) 93-94; 488 (Great Britain, Naut. Alm. Off.) 
94-95; 489 (Duarte) 95; 490 (Ferrier) 95; 491 (Gloden) 96; 492 (Pipping) 96; 493 (Yarden) 
96-97; 494 (Zavrotsky) 97; 495 (Salzer) 97-98; 496 (Kelley) 98; 497 (Palm) 98-99; 498 
(Blewett) 99; 499 (Bouwkamp) 99-101; 500 (Charney) 102; 501 (Harvard Univ., Comput. 
Lab.) 102; 502 (Kober) 103; 503 (Magnus & Oberhettinger) 103-105; 504 (Morgan) 105- 
107; 505 (NBSCL) 107; 506 (U. S. Navy, Off. Res. & Inventions) 107; 507 (Salzer) 107-108; 
508 (Dorodnifsyn) 108; 509 (Placzek) 108; 510 (Skinner-Truscott) 108; 511 (Tsukhanova 
& Salamandra) 109; 512 (Khristianovich) 109; 513 (Deutsche Seewarte) 109-110; 514 
(Japan, Hydrogr. Dept.) 110-112; 515 (Kopal) 112-114; 516 (Rodrigues) 114-115; 517 
(Hickerson) 115-116; 518 (Moretti) 116; 519 (Sachs) 169-170; 520 (Stewart) 170; 521 
(Lambert) 170-172; 522 (Peters) 172; 523 (Bosshardt) 172-173; 524 (NBSCL) 173; 525- 
526 (Rybner) 173-175; 527 (Moon) 175-176; 528 (Beeger) 176; 529 (Kraitchik) 177; 530 
(Lahiri) 177; 531 (Pettit) 178; 532 (Carslaw & Jaeger) 178; 533 (Bell Tel. Labs.) 179; 534 
(Charlier) 180; 535 (Cambi) 180-181; 536 (Campbell & Foster) 181-183; 537 (Doetsch) 
183-184; 538 (Gray, Merwin & Brainerd) 184-185; 539 (Harvard Univ., Comput. Lab.) 
185-186; 540 (Kuznefsev) 186; 541 (Lane & Sweeney) 186; 542-543 (NBSCL) 187-188; 
544 (King) 188; 545 (MacLane) 188; 546 (Pollard & Present) 188; 547 (Irwin) 189-191; 
548 (Tsesevich) 191-195; 549-550 (Japan, Hydrogr. Off.) 195-197; 551 (MIT) 197-198; 
552 (Nat. Adv. Comm. for Aeronautics) 199; 553 (Ellis & Walsh) 296; 554 (Gifford) 296— 
298; 555 (Gloden) 298; 556 (Gupta) 298; 557 (Yarden & Katz) 299; 558 (Akad. N., Matem. 
Institut) 299; 559 (Birge) 300; 560 (Jackman) 301; 561 (Kelley) 301-302; 562 (Nielsen & 
Goldstein) 302-303; 563 (Bennett) 303; 564 (Bouwkamp) 303-304; 565 (Cambi) 304- 
305; 566 (Dwight) 305-306; 567 (Gray & Schelkunoff) 306; 568 (Kaplan) 306-307; 569 
(Kourganoff) 307; 570 (Rice) 307-308; 571 (Tietze) 308; 572 (Rice) 308; 573 (Thomas) 
308; 574 (Peake) 308-309; 575 (France, Service Technique Aéronautique) 309-311; 576 
(Japan, Hydrogr. Off.) 311; 577 (Harvard Univ., Comput. Lab.) 311-314; 578 (Schomann) 
355; 579 (Uhler) 355; 580 (France, Inst. Géogr. National) 355-356; 581 (Beall) 356; 582 
(Italy, Consiglio Naz. d. Ricerche) 356; 583 (Jordan) 356-357; 584 (Kariakin) 357; 585 
(Anon.) 357; 586 (Gloden) 357; 587 (Kraitchik) 357-358; 588 (Salzer) 358; 589 (Lienard) 
358; 590 (Finney) 359; 591 (Finney & Stevens) 359-360; 592 (Fisher & Yates) 360-361; 
593 (Japan, Hydrogr. Dept.) 361-362; 594 (Mainland) 362-363; 595 (Nair) 363-364; 596 
(Emde) 364-366; 597 (Bordoni) 366-367; 598 (Great Britain, ACS) 367; 599 (Harvard 
Univ., Comput. Lab.) 367; 600 (Kitover) 368; 601 (Kober) 368; 602 (Magnus & Oberhet- 
tinger) 368-369; 603 (Yuan) 369; 604 (Lecolazet & Pluvinage) 369; 605 (Ugarov) 369; 606 
(Guldhammer) 369-371; 608 (Goldschmidt & Pitt) 412; 609 (Zimmermann) 412-413; 610 
(Donnay & Hamburger) 413; 611 (Robbins & Smith) 414; 612 (Kerawala & Hanafi) 414; 
613 (Salzer) 414-415; 614 (Woodward) 415; 615 (Bose) 415; 616 (Collie, Hasted & Ritson) 
415; 617 (Ghizzetti) 415-416; 618 (Goldman) 416; 619 (Greenwood & Miller) 416; 620 
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(U. S. Navy, Naval Res. Lab.) 417; 621 (Jaeger) 417; 622 (Zubkov) 417; 623 (Bullard & 
Cooper) 417-418; 624 (Hay & Gamble) 418; 625 (Herget) 418-419; 626 (Italy, Istituto 
Idrogr. d. r. Marina) 419-420; 627 (Lapushkin) 420-422; 628 (Merritt) 462; 629 (Mit- 
rinovich) 463; 630 (Sassenfeld & Tschunko) 463; 631 (Comrie) 464-466; 632 (Larsen) 466; 
633 (Boll) 466-467; 634 (Wittke) 467; 635 (Vega) 467; 636 (Watson) 468; 637 (Milne) 
468-472; 638 (Westberg) 472; 639 (Salzer) 472-473; 640 (Copsey, Frazer & Sawyer) 473; 
641 (Kopal) 473; 642 (Rosser) 474; 643 (Deschodt) 474; 644 (Harvard Univ., Comput. 
Lab.) 474-475; 645 (Macfarlane) 475-477; 646 (Mitchell) 477; 647 (NBSCL) 477-478; 
648 (NBSINA) 478; 649 (Stankiewicz) 479; 650 (Arenberg & Levin) 479; 651 (Harvard 
Univ., Comput. Lab.) 479-482; 652 (Schelkunoff) 482; 653 (Godfrey) 482; 654 (Grébner & 
Hofreiter) 482; 655 (Goodwin & Staton) 483; 656 (Stevenson) 483; 657 (NBSCL) 483-484; 
658 (Italy, Istituto Idrografico) 513; 659 (Dale) 514; 660 (Sainte-Lagiie) 514; 661 (Danjon) 
514-515; 662 (NBSCL) 515-516; 663 (Huskey & Hartree) 516; 664 (Rybner & Sgrensen) 
516-517; 665 (Harvard Univ., Comput. Lab.) 517-518; 666-667 (Lehmer) 518-519; 668 
(Pipping) 519; 669 (Yarden & Katz) 519; 670 (Salzer) 519; 671 (Jeffreys) 519-520; 672 
(Smirnov) 520; 673 (Marcum) 521; 674 (Grover) 521; 675 (Hartree & Johnston) 521; 
676 (Huckel) 522; 677 (Lowan & Horenstein) 522; 678 (Magnus & Oberhettinger) 522- 
523; 679 (NBSCL) 523-524; 680 (Bouwkamp) 524; 681 (Levine & Schwinger) 524; 682 
(Meixner) 524-526; 683 (Hallén) 526; 684 (LeCaine) 526; 685 (Van de Hulst) 526-527; 
686 (Financial Publ. Co.) 527-528; 687 (Italy, Istituto Geogr. Militare) 528-529; 688 
(Hoehne) 529; 689 (Azeredo Rodrigues) 529-531; 690 (Sadler) 531. 


CORRIGENDA 
P. 67, 147, 227, 332, 398, 452, 499, 563. 


MATHEMATICAL TABLES—ERRATA 


118 (Gudermann) 40-41; 119 (Harvard Univ., Comput. Lab.) 41; 120 (Jahnke & Emde) 
41; 121 (Koshliakov) 41-42; 122 (Spence) 42; 123 (Wrench) 42; also Columbia Univ. 
Press 26, Hage 30, Hayashi 42, Jérgensen 62, Kerawala & Hanafi, 24, King 29, Miiller, 
Rajna & Gabba 20, MIT 38-39, NBSCL 25, ProkeS 20, Reiz 26, Schulze 20, U. S. Hydrogr. 
Off. 35, Wrench 18-19, Ziaud-Din 24; 124 (Airey) 116-118; 125 (Gloden) 118; 126 (Walther) 
119; 127 (Yarden) 119-120; also Ageton 145, Anan’ev 93, Bouwkamp 100, Cazzola 84-85, 
Chebyshev 97-98, Editor 84, 1. -5, -6, 106, 1. 19-20, Ferrier 95, Fontoura & Penteado 145, 
Germany, Deutsche Seewarte 138-140, Great Britain, N. A. O. 95, Magnus & Oberhettinger 
104-105, NBSCL 88-89, Skinner (Truscott) 108, U. S. Navy 107, Whittaker & Robinson 
98; 128 (Airey) 199, 129 (Becker & Van Orstrand) 200-201, 130-131 (Dwight) 201; 132 
(Gauss) 201-203, 133 (NBSCL) 203; 134 (Okaya) 203; 135-137 (Peters) 203-205; also 
Cajori 224, Cambi 181, Doetsch 183, Felkel 171, Ferrari 194, Kraitchik 177, Lambert 
171-172, Lane & Sweeney 186, Lehmer 227, Moon 175, Peters 172, Pettit 178, Russell & 
Shapley 192, Rybner 174, Tsesevich 193; 138 (Adams) 314; 139 (Comrie) 314; 140 (Hardy 
& Wright) 314; 141 (Jahnke & Emde) 314-315; 142 (NBSCL) 315; 143 (Peters) 315; 144 
(U. S. Hydrogr. Off.) 315-316; also Airey 261, Akushskii & Ditkin 299, Bertrand 261-262, 
Dale 262, DeMorgan 308, Ditkin 299, Dwight 262, FMR 276-277, France 309, Gauss 
262, Gifford 297, Glaisher 263, Gossot 263, Greenhill 263, Hammer 308, Hancock 263, 
Harvard Univ., Comput. Lab. 313, Hayashi 263-266, Heuman 266, Hippisley 266-267, 
Innes 266-267, Jahnke & Emde 267-268, Legendre 268-274, Lévy 274-275, Liusternik 
299, Loxhay 280, Meissel 275-276, Merfield 269, Moore 276, Nagaoka & Sakurai 276, 
Pearson 306-307, Plana 277-278, Potin 278, Rosenbach, Whitman & Moskovitz 278, 
Runkle 278, Samoilova-fakhontova 278-279, Schlémilch 279, Soldner 308, Spenceley 280, 
Verhulst 280-281, Wayne 281, Weigand 308, Yarden & Katz 299; 145 (Becker & Van 
Orstrand) 371; 146 (Fisher) 371; 147 (Kraitchik) 372; 148 (Oughtred) 372; also France, 
Inst. Géog. Nat. 356, Guldhammer 371, Jahnke & Emde 365-366, Kariakin 357, Pitiscus 
391-392, Rheticus 395-396, Ugarov 369, Uhler 355; 149 (Adams) 423; 150 (Glazenap) 423; 
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ik 151 (Jahnke & Emde) 423; 152 (Kelley) 423; 153 (Richmond) 423-424; 154 (BAASMTC) 
uto 424; 155 (U. S. Coast and Geod. Survey) 424-425; also Bertrand 449, Buckingham 413, 
fit- Callet 426, Davis 424, Davis & Kirkham 450, Goldman 416, Gray & Mathews 449-450, 
66; Hay & Gamble 418, Herget 419, Kerawala & Hanafi 414, Lapushkin 422, Meissel 449- 
ne) 450, U. S. Navy 417, Zimmermann 412; 156 (Holtappel) 485; 157 (Keiley) 485; 158 (Pear- 
73; son) 485; 159 (U. S. Coast and Geod. Survey) 485; also Boll 466-467, Cunningham 486, 
ut. Hacker 493, NBSINA 478, Newman 477, Powell 477, Seewald 472; 160 (D. H. Lehmer) 
78; 531; 161 (D. N. Lehmer) 531-532; also Boll 518, Bouwkamp 525, Briggs 516, Dale 514, 
ard Feller 520, FMR 551, Gauss 551, Hallén 526, Huckel 522; Lehmer 518-519; Magnus & 
r& Oberhettinger 523, Meixner 525, Rheticus 556-557, Watson 549-550. 

84; 
on) MATHEMATICAL TABLES—ERRATA, AUTHOR INDEX 
nee Adams 314, 423 Gloden 118 
672 Ageton 145 Goldman 416 
21: Airey 116-118, 199, 261 Gossot 263 
22- Akushskii & Ditkin 299 Gray & Mathews 449-450 
682 Anan’ev 93 Great Britain, N.A.O. 95 
27; Becker & Van Orstrand 200-201, 371 Greenhill 263 
688 Bertrand 261-262, 449 Gudermann 40-41 
Boll 466-467, 518 Guldhammer 371 
Bouwkamp 100, 525 Hacker 493 
Brandenburg 43 Hage 30 
BAASMTC 424 Hallén 526 
Briggs 516 Hammer 308 
Buckingham 413 Hancock 263 
Cajori 224 Hardy & Wright 314 
Callet 426 Harvard Univ., Comput. Lab. 41, 313 
de) Cambi 181 Hay & Gamble 418 
Liv. Cazzola 84-85 Hayashi 42, 263-266 
ler, Chebyshev 97-98 Herget 419 
gr. Columbia Univ. Press 26 Heuman 266 
er) Comrie 314 Hippisley 266-267 
85, Cunningham 486 Holtappel 485 
45, Dale 262, 514 Huckel 522 
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rdy Editor 84, 106 Kerawala & Hanafi 24, 414 
144 Felkel 171 King 29 
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uss Ferrari 194 Kraitchik 177, 372 
63, Ferrier 95 Lambert 171-172 
67, Fisher 371 Lane & Sweeney 186 
nik FMR 276-277, 551 Lapushkin 422 
76, Fontoura & Penteado 145 Legendre 268-274, 306-307 
78, France 309, 356 Lehmer, D. H. 227, 518-519, 531 
80, Gauss 201-203, 262, 551 Lehmer, D. N. 531-532 
Jan Germany, Deutsche Seewarte 138-140 Lévy 274-275 
ice, Gifford 297 Litisternik 299 
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MIT 38-39 Rybner 174 

Meissel 275-276, 449-450 Samoilova-fakhontova 278-279 
Meixner 525 Schlémilch 279 

Merfield 269 Schulze 20 

Moon 175 Seewald 472 

Moore 276 Skinner (Truscott) 108 

Miiller, Rajna & Gabba 20 Soldner 308 

Nagaoka & Sakurai 276 Spence 42 

NBSCL 25, 88-89, 203, 315 Spenceley 280 

NBSINA 478 Tsesevich 193 

Newman 477 Ugarov 369 

Okaya 203 Uhler 355 

Oughtred 372 U. S. Coast & Geod. Survey 424-425, 
Pearson 306-307, 485 485 

Peters 172, 203-205, 315 U. S. Hydrogr. Off. 35, 315-316 
Pettit 178 U. S. Navy 107, 417 

Pitiscus 391-392 Verhulst 280-281 

Plana 277-278 Walther 119 

Potin 278 Watson 549-550 

Powell 477 Wayne 281 

ProkeS 20 Weigand 308 

Reiz 26 Whittaker & Robinson 98 
Rheticus 395-396, 556-557 Wrench 18-19, 42 

Richmond 423-424 Yarden 119-120 

Rosenbach, Whitman, & Moskovitz 278 Yarden & Katz 299 

Runkle 278 Ziaud-Din 24 

Russell & Shapley 192 Zimmermann 412 


UNPUBLISHED MATHEMATICAL TABLES 
There is no number 66 


63 (Wrench) 42-43; 64 (Porter) 43; 65 (NBSCL) 43-44; 67 (Poulet) 120; 68 (Salzer) 
120-121; 69 (Salzer) 205; 70 (Spenceley) 205; 71 (Salzer) 205-206; 72 (Salzer) 316; 73 
(Dixon) 316-317; 74 (Salzer) 317; 75 (Univ. California, Statist. Lab.) 372-373; '76 (Faucher) 
425; 77 (Bower) 425-426; 78 (Greenwood) 426; 79 (Bartot Res. Foundation) 485; 80-81 
(Gloden) 486; 82 (Poletti) 532; 83 (Salzer) 532. 


AUTOMATIC COMPUTING MACHINERY 


Technical Developments: See under Articles, Alt; Bloch, Campbell & Ellis; Eckert; 
Haeff; Harrison & Malone; McCann; Snyder & Livingston; Verzuh. 

Discussions: S. Lubkin, Decimal point location in computing machines, 44-50; J. L. 
McPherson, Applications of large-scale high-speed computing machines to statistical work, 
121-126; H. F. Mitchell, Jr., Inversion of a matrix of order 38, 161-166 [also as an Article]; 
R. F. Clippinger, Airflow problem planned for the ENIAC, 206-207; B. L. Hicks 
& H. G. Landau, Nonlinear parabolic equations, 207-208; J. V. Holberton, Laminar 
boundary layer flow in a compressible fluid, 208; J. H. Levin, On the approximate solution 
of a partial differential equation on the differential analyzer, 208-209; M. Lotkin, Com- 
putation of the airflow about a cone cylinder, 209-210; R. M. Bloch, R. V. D. Campbell & 
M. Ellis, General design considerations for the Raytheon computer, 317-323; Dorrit 
Hoffleit, A comparison of various computing machines used in the reduction of Doppler 
observations, 373-377; Florence Koons & S. Lubkin, Conversion of numbers from decimal 
to binary form in the EDVAC, 427-431; E. F. Moore, A new general method for finding 
roots of polynomial equations, 486-488; J. L. Meriam, Procedure for the machine or nu- 
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merical solution of ordinary linear differential equations for two-point linear boundary 
values, 532-539. 

Bibliography: Burks, Goldstine & Neumann; Comrie; Goldstine & Neumann; Hartree, 
50-57. Anon.; Burks, Goldstine & Neumann; Cesareo; Hartree; Juley; Moore School of 
Electr. Engineering, 126-132. Aiken & Hopper; Anon.; Berkeley; Burks; Comrie; Condon; 
Curry & Wyatt; Koehler; Russo; Seeber; Lilley; Tjanste; Tumbleson, 210-213. Burks; 
Djure & Welin; Goldstine; Great .Britain, Post Office Dept.; Huskey; IBM, 323-326. 
Aiken & Hopper; Anon.; Bigelow, Pomerene, Slutz & Ware; Couffignal; Haeff; Lion; Moore 
School Electr. Engin.; Inst. Radio Engin., 377-380. Brillouin; Eckert; Fuller; Harvard 
Univ., Comput. Lab.; Inst. Adv. Study; Irwin; Korn; MacLeish, O’Neal & Tyler; Moore 
School Electr. Engin.; O’Neal & Tyler; Williams & Kilburn; U. S. Air Force, 431-445. 
Booth & Britten; Cohen & Keye; Coombs; Hamilton; Hoberg & Ulman; Page; Patterson, 
Snyder, Tabor & Travis; Pfeiffer; Wilkes; Wilkes & Renwick, 488-490. Anon.; Arenberg; 
Boothroyd & Cherry; Brainerd & Sharpless; Bronzo & Cohen; Chu; Davis; Goldstine & 
Neumann; Hartree, Newman, Wilkes, Williams, Wilkinson & Booth; McCool; Sharpless; 
Wiener, 539-544. 

News: Amer. Acad. Arts Sci. 216, Amer. Phys. Soc. 326-327, Amer. Statist. Assoc. 134, 
Assoc. Computing Machin. 57, 132-134, 216, 380, 445-446, 490, 544-545, Automatic relay 
computer 216, BAAS 132, Electronic digital computing in England 213-216, Institute 
Radio Engins. 491-492, Institute Teachers Math. 380, Institution Electr. Engins. 57, 
IBM 381, IBM Selective Sequence Electronic Calculator 216-217, Intern. Sci. Radio 
Union 327, Mass. Inst. Techn. 492, Math. Assoc. Amer. 58, Math. Computing France 58, 
NBSINA 545-547, NBSMDL 327-328, Nat. Electronics Conference 134, Naval Res. Lab. 
446, Office Naval Res. 446, Personal 328, Swedish State Board Computing Machinery 
492; Symposia on modern Calc. Mach. 381-388. 


OTHER Arps TO COMPUTATION 


G. R. Stibitz, A new class of computing aids 217-221. 
G. R. Stibitz, Electric root-finder 328-329. 

Bibliography: Akushskii; Bartlett; Baudouin; Bergman; Berry & Pemberton; Bruk; 
Eckert; Fiirth & Pringle; Globe-Hilsenrath; Herget; Herrick; Kind-Schaad; King; Renner; 
Reynolds; Shaffer, Schomaker & Pauling; Tolstov, 58-61. Anon.; Baracket; Benjamin; 
Black & Olds; Booth; Great Britain, Ministry of Supply; Harbert; Hooks; Johnson; 
Lafferty; Manley; Meyer zur Capellen; Psychometric Soc. & IBM; Vajda; Waddington; 
Weissenstein, 134-138. Adcock; Akushskii; Anon.; Gavrilko; Leibin; Packh; Weems, 
221-222. Babcock; Goldberg & Brown; Murray; Tiller, 329-330. Camp; Carle; Comrie; 
Gilbert; Othmer, Josefowitz & Schmutzler; Purdie; Redheffer; Richardson; Tea; Wilkas; 
Wittke, 388-390. Bothel; Goldberg; Hardy & Dench; Lo-Ho; Nachod; Rubbert; Strohm 
& DeGroot; Svoboda, 446-448. Hacker; Kiely; Reinhardt; Willers, 492-494. Brown; 
Carter & Sadler; Knauss; Richards & Sadler; Vand; Wittke, 548-549. 


NoTEs 


80. R.C.A., Certain gear-ratio tables, 61; 81. R.C.A. & M. S. Corrington, Guide (MTAC 
no. 7),:62; 82. W. F. Willcox, Herman Hollerith, 62-63; 83. H. E. Salzer, Interpolation of 
functions tabulated at fixed points, 63-64; 84. The National Applied Mathematics Lab- 
oratories, 64-65; 85. NAML, 65; 86. Russia and mathematics, 65-66; 87. J. B. Parker & 
D. H. Sadler, German altitude and azimuth tables,” 138-141; 88. R.C.A., Harry Bateman, 
141-142; 89. H. S. Uhler, The magnitude of higher terms of the Lucasian sequence, 142-143; 
90. Mathematical Table Makers, 143-144; 91. R.C.A., Portraits of table makers, 144-145; 
92. C. H. Smiley, Portuguese navigation tables, 145; 93. A film of part of Kulik’s Magnus 
Canon for sale, 222; 94. R.C.A., Lambertian or lambda function, 222-225; 95. R.C.A., 
Rydberg interpolation table, 330; 96. R.C.A., Bartholomaus Pitiscus (1561-1613), 390-397; 
97. Fritz Emde, 397-398; 98. R.C.A., Mersenne numbers, 398; 99. R.C.A., Daniel Friedrich 
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Ernst Meissel, 449-451; 100. A. Ferrier, A new factorization of 2* + 1, 451; 101. R.C.A., 
Newman’s Mathematical Tracts, 451; 102. Tables of x tan x, 451; 103. Editors, Arnold Noah 
Lowan, 494-495; 104. L. G. Pooler, Roots of certain transcendental equations, 495-496; 
105. A. Ferrier, Note on the factors of 2" + i, 496-497; 106. Editors, George Neville Watson 
—table maker, 549-550; 107. D.H.L., The lemniscate constant, 550-551; 108. R.C.A., 
Newman and George Eliot 551-552; 109. R.C.A., Rheticus, with special reference to his 
Opus Palatinum, 552-561. 
PORTRAITS 


F. Emde, no. 25;'A. N. Lowan, no. 27; D. F. E. Meissel, no. 26; G. N. Watson, no. 28. 


QUERIES 


24. C. W. Horton, Integral of a Struve function, 66; 25. R.C.A., Russian Bessel function 
tables, 66; 26. R.C.A., Some clothoid or Euler spiral tables, 146; 27. R.C.A., First table of 
tan- x, 226-227; 28. A. Fletcher, Moseley on definite integrals, 331; 29. R.C.A., Pitiscus 
tables, 398; 30. R.C.A., Girard and Snell tables, 451-452; 31. R.C.A., Japanese World 
War II editions of mathematical tables, 498; 32. French and Russian translations of a 
Vega-Bremiker table, 561. 

QUERIES—REPLIES 


There is no no. 36. 


33. J. Ginsburg, Editor, Portraits and biographies of British mathematical table makers, 
66-67; A. J. Thompson, Henry Briggs, 67; 34. R.C.A., An Engel table, 146; 35. J. Todd, 
Tables of tan! (m/n), 227; 37. N. G. W. H. Beeger, Hendrik Anjema, 331-332; 38. Log 
log tables, 398; 39. R.C.A., Some clothoid or Euler spiral tables, 452; 40. A. Fletcher, Mrs. 
G. O. Savage, Editor, Pitiscus tables, 498-499; 41. J. W. Wrench, Jr., Girard and Snell 
tables, 561-562; 42. W. D. Morgan, Editor, Pitiscus tables, 562-563. 


SUBJECT INDEX 





Accounting 389, 539, 548 

Acoustics 37, 134-135, 416, 540 

Actuarial Tables 466 

Addition and Subtraction Logs 20, 85, 513 

Aerodynamics see Fluid Mechanics 

Air Navigation see Navigation 

Airy’s Integral 337, 339, 418 

Algebraic Equations 13-16, 453-462, 469, 
486-488 

Algebraic Equations, Mechanical Solutions 
of 60-61, 81, 137-138, 221, 299, 329-330, 
436, 494, 545 

Aliquot Series 120 

Alpha Functions 112-114 

Altitude-Azimuth 31-36, 59, 110-112, 114— 
116, 138-141, 309-311, 530 

Analogue Machines 56-57, 61, 127, 208-209, 
212, 217, 323, 327-329, 382, 441, 491- 
494, 501-513, 539-540, 542-544 

Anger Functions 103, 365-366 

Annuity 29-30, 408-412, 465-466 

Antenna Theory 303-304, 524 

Antilogarithms 84-85, 463, 465, 514 

Arc Cosine, etc., see Cosine, Inverse, etc. 


Arc Length 20, 85, 94, 171, 463, 513, 562 
Astronomy 21, 31-32, 85, 109-116, 189-195, 
418-419, 440-441, 513, 526, 531, 556 

Astrophysics 526-527 

Automatic Computing Engine (ACE) 214- 
215 

Automatic Relay Computer (ARC) 216 

Automatic Sequence Controlled Calculator 
(ASCC) 54-56, 128, 130-131, 149, 210, 
213, 314, 377-378, 493, 518 

Azimuth Computer 58-59 


Babbage’s ‘‘Engines” 53-54, 210, 215, 432, 
493 

Ballistics 37-40, 127, 197-199, 252-258, 311- 
313, 373-377, 437 

Bell Telephone Relay Computer 1-13, 54, 
69-84, 207-208, 210, 375-377, 433 

Bernoulli Functions 368 

Bernoulli Numbers 514 

Bessel Functions 
Asymptotic Expansions 335, 365 
Combinations of 305-306, 367, 415, 417 
Integrals of 303, 308, 340, 417-418, 550 
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Logarithms of 62 
Modified 105 
of Complex Argument 25, 333, 335, 352 
of First Kind and Integral Order 41, 65, 
88, 102-103, 107, 116-117, 180-181, 
183, 185-186, 199, 253, 333, 335, 337, 
339, 366-367, 416, 449-450, 474- 
475, 514, 549 
of Fractional Order 187, 339, 466-467, 
484, 514, 550 
of Purely Imaginary Argument 105-107, 
365-366, 415, 478, 514, 549 
of Second Kind and Integral Order 107, 
187-188, 203, 315, 333, 549 
Spherical 26 
Turning Points of 26 
Zeros of 26, 102, 305-306, 340, 351-352, 
366, 496, 514, 549-550 
Bessel Interpolation Coefficients 466 
Beta Function 485 
Binary Computers 48, 128, 287-288, 379, 
427-431, 441-443, 488, 491-492, 540 
Binary Quadratic Forms 22-23, 43 
Binomial Coefficients 87, 171, 339, 463, 472, 
514 
Binomial Congruences 335 
Binomial Distribution 179 
Biography and Bibliography, Portraits: 
Airey 143-144, Anding 143, Anjema 
331-332, Babbage 143, Barlow 66-67, 
Bateman 141-142, Bauschinger 143, 
Becker 143, Bessel 143, Bierens de Haan 
143, Borda 143-144, Briggs 67, Brown 
143, Biirgi 143-144, Burrau 143, Charlier 
180, Chebyshev 144, Cohn 143, Comrie 
143-144, Cunningham 143-144, Dase 
143, Davis 143, Dickson 143, Dodson 67, 
Dwight 143, Eckert 145, Eliot, George 
551-552, Emde 144, oppo. 333, 397-398, 
Farley 66, Filipowski 66-67, Gauss 85- 
86, Gifford 296-297, Glaisher (J.) 143, 
Glaisher (J. W. L.) 143-144, Gray 66-67, 
Hoéne-Wrofski 143, Hollerith 62-63, 
Hoppe 143, Hutton 143, Jacobi 143, 
Kepler 143-144, Kraitchik 143, Kruse 
144, LaLande 143, Lambert 145, 
Legendre 143, Lehmer (D. H.) 143-144, 
Lehmer (D. N.) 143, Lodge 143, Lohse 
143, Lommel 143, Lowan 143, oppo. 453, 
494-495, Markov 143, Martin 143, 
Meissel oppo. 399, 449-451, Miller 
143, Miiller 20, Napier 143, Newman 
551-552, Nielsen 143, Ostrogradskii 145, 
Oughtred 372, Pearson 143, Peirce 143, 
Peters 143-144, Pitiscus 390-397, 498- 


499, 562, Rheticus 552-562, Rivard 143, 
Rydberg 330, Sang 143, Segner 226, 
Sharp 143, Sheppard 143, Stevin 143- 
144, Stieltjes 143, Tallqvist 143, Thomp- 
son 143, Turner 143, Uhler 143, 
Viéte 143, Watson oppo. 501, 549-550, 
Wright 223-225, 372. Also 50-57, 58-61, 
126-132, 134-138, 210-213, 221-222, 
250-260 (Elliptic Functions), 323-326, 
329-330, 377-380, 388-390, 431-445, 
446-448, 488-490, 492-494, 539-544, 
548-549 

Bivariate Interpolation 524 

Bivariate Normal Distribution 372-373, 474 

Bond Values 29-30, 408-412, 527-528 


Calculating Machines see Computing Ma- 
chines 

CIT (Cal. Inst. Techn.) Electronic Ana- 
logue Computer 501-513 

Chebyshev Polynomials 17, 64, 103, 119- 
121, 300, 405-407, 471 

Chebyshev Quadrature 97 

Chi Square Distribution 302, 363, 466, 515, 
519-520 

Christoffel Numbers 26, 41, 416, 473 

Class Number 22 

Clausen Integral 368 

Clothoid 146, 452 

Coding 48-51, 54-56, 78-79, 320, 444, 488- 
490, 541-542, 545 

Complex Numbers 493, 516-517 

Computational Methods 131-132, 216, 299, 
381-382, 400-407, 434-440, 443, 468- 
472, 486-488, 493-494, 532-539, 545- 
547 

Computing Machines 1-13, 44-58, 69-84, 
121-138, 149-161, 206-222, 286-295, 
299, 317-330, 373-382, 388-390, 426- 
448, 486-492, 501-513, 539-548 

Cone Function 103, 105 

Conformal Mapping 103, 250, 368, 493, 545- 
546 

Congruences 96, 177, 357, 486 

Constants 20, 42-43, 94, 254, 302, 369, 463, 
466-467, 478, 513-515 

Continued Fractions 96, 519 

Conversion 20, 94, 340, 463, 465-467, 513- 
515 

Coordinates, Polar 94, 340 

Correlation 60, 371, 389 

Cosecant, Log of 19, 371, 513 

Cosecant, Natural 19, 93-94, 172, 371, 513- 
515 

Cosecants, Sums of 549-550 
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Cosine, Hyperbolic 86-87, 93, 170-171, 205, 
336, 463, 465-466, 514 

Cosine, Hyperbolic Inverse 517-518 

Cosine, Inverse 85, 87, 463, 465 

Cosine, Log of 19-20, 85-87, 355-356, 370, 
463-467, 513-514 

Cosine, Natural 19-20, 85-87, 93-94, 146, 
172, 370, 413, 425-426, 463, 466-467, 
485, 513-516 

Cosine Integral 104, 336, 338, 417, 479- 
482, 514, 526 

Cosine Integral, Iterated 526 

Cotangent, Hyperbolic 88-89, 296 

Cotangent, Inverse 85, 87, 465 

Cotangent, Log of 19-20, 85, 355-356, 371, 
463-467, 513 

Cotangent, Natural 19-20, 85, 88-89, 93- 
94, 172, 296, 371, 463, 466-467, 513 

Crystallography 60, 136, 413, 540 

Cube Root 19, 84-85, 98, 302, 340, 370, 463, 
466, 494, 513-514 

Cubes 19, 84-85, 171, 370, 463, 466, 513-514 

Cubes, Endings of 298 

Cubic Equations 171, 388, 467 

Cubic Forms 24 

Cuboids 177 

Cybernetics 543-544 

Cylinder Functions 103-104, 369 

Cylinder Functions, Elliptic 101 

Cylinder Functions, Parabolic 104, 365 


Decimal Equivalents of Fractions 19, 462, 
514 

Derivatives, Numerical 414-415 

Desk Computers 136-137, 493-494, 548-549 

Determinants 469 

Difference Engine see Babbage’s ‘‘Engines’’ 

Difference Equations 471 

Differential Analyzer 54, 56-57, 127, 137, 
208-209, 212, 491, 493 

Differential Equations 61, 81, 132-133, 160, 
184-185, 206-210, 304-306, 330, 436- 
437, 470, 494, 501-507, 511, 522, 524— 
526, 532-540, 548 

Differentiation, Numerical see Numerical 
Differentiation 

Digamma Function 472 

Digital Machines 1-13, 44-58, 69-84, 121- 
134, 149-161, 206-208, 210-217, 286- 
295, 317-328, 373-381, 427-446, 488- 
492, 540-542, 544-548 

Diophantine Equations 95, 177, 333 

Divisors 334, 337, 339 

Dixon Function 229, 249-250 

Doppler Method 373-377 

Duodecimal Tables 23, 357, 388 
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Eigenvalues 524-526 

Elasticity 29, 93, 357, 368, 483, 511, 545 

Electrodynamics 369 

Electromagnetism 237-239, 250, 253, 256, 
308-309, 483-484 

Flectronic Delay Storage Automatic Cal- 
culator (EDSAC) 490, 542 

Electronic Discrete Variable Computer 
(EDVAC) 54, 57, 129-131, 213, 379-380, 
382, 427-431, 437, 440-441, 444-445, 
489-490, 540-541 

Electronic Numerical Integrator and Com- 
puter (ENIAC) 1-4, 7, 53-54, 57, 83- 
84, 126, 129, 131, 133, 185, 206-208, 210, 
212-213, 323-326, 330, 375-377, 433, 
437, 440, 446, 493, 518, 540-541 

Elliptic Functions 89-92, 104, 229-281, 334- 
335, 365, 369, 449, 514, 521, 551 

Elliptic Functions, Cylinder see Cylinder 
Functions, Elliptic 

Elliptic Integrals 27-28, 89, 99, 104, 229, 
231-281, 304-307, 369 

Emden Functions 336, 339 

Equation Solvers 61, 137-138, 221, 323, 329- 
330, 447, 491 

Equations, Algebraic see Algebraic Equations 

Equations, Transcendental see Transcen- 
dental Equations 

Error Function, see Probability Integral 

Euler Numbers 467, 514 

Euler Product 340 

Everett Interpolation Coefficients 187, 339, 
366, 418-419, 470, 472, 477-478 

Exponential Function 85, 87, 103, 173, 334, 
371, 426, 463, 465-466, 485, 514 

Exponential Integral 104, 336, 338, 514, 526 

Exponents 176, 357 


Factor Tables (general) 85, 87, 95, 171-172, 
222, 334, 336, 339, 462-463, 466-467, 532 

Factor Tables (special) 21, 23, 95, 118-119, 
299, 357, 451, 486, 496, 518-519 

Factorials 85, 205, 340, 355, 357, 466 

Factorials, Logs of 85, 466, 514 

Factorials, Ratios of 26 

Factorials, Reciprocals of 514 

Farey Series 333, 339 

Fermat’s Converse 314, 451, 496-497, 518- 
519 

Fermat's Last Theorem 493 

Fibonacci Sequence 96-97, 119-120, 299 

Figurate Numbers 171 

Fisher’s z-Test 519-520, 524 

Fluid Mechanics 36-40, 109, 116, 136, 197- 
199, 206-207, 209-210, 299, 311-313, 
356, 369, 511, 522, 545 
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Formulae 94, 251, 369, 466, 474, 513-515, 
521-523 

Fourier Coefficients 241, 251, 413, 471 

Fourier Integrals 21, 181-183, 475-476 

Fourier Series 16, 82, 251, 434 

Fourier Transformation 104 

Fourier Transformer 61 

Fresnel Integrals 104, 146, 336, 417, 467, 479 

Functions (special) see under Alpha, Anger, 
Bernoulli, Bessel, Beta, Cone, Cylin- 
der, Digamma, Dixon, Elliptic, Emden, 
Exponential, Gamma, Hankel, Har- 
monic, Hermite Probability, Hyperbolic, 
Hypergeometric, Intensity, Jacobi Zeta, 
Kelvin, Kummer, Laguerre, Lamb- 
da, Legendre, Lommel, Lommel-Weber, 
Mathieu, Planck’s, Polygamma, Psi, 
Radiation, Ramanujan’s, Riemann Zeta, 
S and T, Scattering, Spherical, Sphe- 
roidal Wave, Stokes’, Struve, Symmetric, 
Tetragamma, Theta, Toroidal, Tri- 
gamma, Trigonometric, Weber, Wedge, 
Weierstrass Elliptic, Whittaker, Zeta. 


Galois Fields 97 

Games 389, 490-491 

Gamma Function 98, 103, 307, 317, 336, 
338, 369, 466, 472, 478, 514, 521, 549, 551 

Gamma Function, Complex Argument 415- 
416 

Gamma Function, Incomplete 104, 302 

Gamma Function, Log of 340, 465, 549 

Gaussian Quadrature 97, 470 

Gear Ratios 61, 462 

Geodesic Tables 86-87, 513, 528-529 

Grades 94 

Graeffe Process 487-488, 494 

Graphs 138, 299, 306, 388 

Gravitation 20, 417-418 

Gregory-Newton Interpolation Coefficients 
358, 470 

Gudermannian 466, 514 


Hankel Functions 299 

Harmonic Analysis 61 

Harmonic Analyzer 137-138 

Harmonic Functions 61, 182 

Harmonic Series 159-160 

Harmonic Syntheses 413 

Heat Equation 109, 240, 299, 511 

Hermite Polynomials 103, 167-169, 180, 416, 
471, 473, 521 

Hermite Probability Function 338 

Hydraulics, Hydrodynamics, see Fluid Me- 
chanics 

Hyperbolic Functions 371, 517-518 


Hypergeometric Function 103-104, 248- 
249, 307, 338, 369, 522 

Hypergeometric Function, Confluent 102, 
336, 340, 366, 369, 520, 522-524 


Ideals, Reduced 336, 339 

Indices 372 

Integral Equations 104, 390 

Integrals, Misc. Definite 21, 29, 188, 307- 
308, 331, 369, 418, 472-474, 479-482, 
524, 526-527 

Integrals, Multiple 479, 550 

Integrals, Tables of 87, 188, 251, 466, 482, 
514 

Integraph 138, 306, 390 

Integration, Numerical see Numerical In- 
tegration 

Integration of Differential Equations see 
Differential Equations 

Intensity Function 484 

Interest, Compound 408-412, 465-466, 514, 
528 

International Business Machines (IBM) 
Selective Sequence Electronic Calculator 
216-217, 326-328, 330, 367, 381-382, 540 

international Business Machines (IBM) 
Sequence Relay Calculator 83, 133, 149- 
161, 207-208, 375-377 

Interpolation 26, 63-64, 82, 85, 205-206, 
212, 330, 339, 358, 361-362, 366, 408- 
412, 414-415, 418-419, 466, 468-472, 
481, 494, 515, 519 

Isograph 138, 491 


Jacobi Nome q 229, 234-235, 256 

Jacobi Polynomials 103 

Jacobi Zeta Function 90-91, 241, 244-246, 
256-257 

Julian Year 514 


Kelvin Function 340 
Kummer Function 104 


Lagrangean Coefficients 98, 302, 414, 466, 
469, 472 

Lagrangean Polynomials 469 

Laguerre Functions 104, 365 

Laguerre Polynomials 103, 532 

Lambda Function, see Lambertian 

Lambertian 170, 222-225 

Laplace Coefficients 229, 237, 252, 470 

Laplace Equation 106, 239, 299, 343-346 

Laplace Transforms 104, 178, 183-184, 476, 
523 

Lattice Points 22-23 

Least Squares 189-191, 300, 302-303, 470- 
471, 513 
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Legendre Functions 103, 252, 365, 525 

Legendre Polynomials 16-18, 99, 103, 167, 
186, 334, 337, 339, 369, 470-472, 484, 
514, 521 

Lemniscate Constant 550-551 

Linear Equations 13-16, 60-61, 81, 436, 
453-462, 545 

Linear Forms 21-22 

Linkages 222, 447-448 

Logarithmic Integral 104, 308 

Logarithms, Common 19-20, 84-85, 87, 340, 
355, 365, 370, 425-426, 463, 465, 467, 
513-514 

Logarithms, Natural 85-87, 142-143, 171, 
360, 372, 463, 465, 513-514 

Logarithms, Addition and Subtraction see 
Addition and Subtraction Logarithms 

Logs of Primes 463 

Log Log 398, 426, 548 

Lommel Functions 103, 186 

Lommel-Weber Functions 103, 336 


M, the Modulus of Common Logarithms 20, 
85, 465 

Magic Squares 531-532 

Magnetism 99, 237-239 

Mark I 58, 69, 79, 83, 432-446, 493, 545 

Mark II 2, 58, 83-84, 210, 433-434, 437- 
438, 446, 493 

Mark III 446, 492, 499, 540 

Mathieu Functions 103-104, 366 

Matrices 159-166, 399-400, 436-437, 453- 
456, 469, 473, 544 

Mellin Transform 475-476 

Memory Tube 281-286 

Mersenne Numbers 142, 398 

Meteorology 85, 102 

Monte Carlo Problem 544, 546-547 

Mortality Tables see Actuarial Tables 

Motor-boating 508 

Multiplication Tables 299, 355, 465 


Navigation 32-36, 43-44, 58-59, 138-141, 
145, 195-197, 221, 309-311, 361-362, 
369-371, 419-422, 529-531 

Network Analyzer 222, 417, 539 

Newton Interpolation Coefficients 470, 472 

Nomographs 135-137, 173-175, 390, 447, 
514, 529-531, 547 

Normal Distribution 515 

Normal Ordinate 466 

Number Theory 21-24, 95-97, 142-143, 
176-177, 298-299, 333-340, 357-358, 
451, 468, 486, 496-497, 518-519, 531- 
532, 550-551 


Numerical Differentiation 466, 468-470, 494 
Numerical Integration 60, 81-82, 97, 107- 
108, 132, 466, 468-470, 494, 516 
Numerical Integration of Differential Equa- 
tions see Differential Equations 
Numeroscope 432 


Optics 99-101, 412, 482, 524 
Opus Palatinum 554-558 
Orthogonal Polynomials 300, 360, 369, 472 


Partitions 340, 519, 550 

Pearson Integral 335 

Pell Equation 335 

Physical Tables 85, 253-254, 513 

Physics 20, 107-109, 175-176, 237-239, 250, 
253, 256, 308-309, 369, 483-484, 522- 
523, 547 

Pi 18-19, 20, 172 

Pi, Multiples of 85, 463, 513 

Pi, Powers of 42-43, 85 

Planck’s Function see Radiation Function 

Planimeter 138 

Poisson Distribution 301 

Poisson-Laplace Equation 61 

Polar Coordinates 172-173 

Polygamma Function 334, 336 

Polygons 561 

Polyhedra, Regular 467, 561 

Polynomial Approximations 400-407, 426 

Powers, Fifth, Decomposition into 336, 339 

Powers, Fourth 87, 466, 514 

Powers, Fractional 466, 485 

Powers, Higher than Fourth 85, 167-169, 
171, 334, 337, 339, 463, 466, 514 

Powers, Negative 466, 485 

Powers, Sums of 358, 514 

Prime Pairs 532 

Primes, Distribution of 532 

Primes, Lists of 85, 87, 171-172, 357, 467, 532 

Primes, Logs of 463 

Primes, Number of 95 

Primes of Special Forms 119, 357-358 

Probability Integral 87, 98, 178, 180, 301- 
302, 316-317, 336-337, 339, 359-360, 
365, 466, 474, 514 

Probability Integral, Integrals and Deriva- 
tives of 336, 338, 521 

Probability Integral, Inverse of 180, 301 

Psi Function, see Digamma Function 

Punch Card Machines 2-4, 7, 60, 135, 150- 
151, 221, 325, 352-355, 376-377, 381, 433 

Punch Card Methods 27-29, 60, 135, 255, 
352-355, 456-462 

Punch Card Tables 60 
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Quadratic Forms 22-23, 43 
Quadrature see Numerical Integration 
Quantum Mechanics 27-29, 188 
Quarter Squares 514 

Quartic Congruences 21, 96, 486 


Radiation 107, 484 

Radiation Function 175-176 

Radioactivity 548 

Railway Transition Curves see Clothoid 

Ramanujan’s Function r(m) 23, 177, 298, 
468, 550 

Random Numbers 547 

Raytheon Computer 286-295, 317-323, 327, 
540 

Reciprocals 84-85, 87, 94, 169-170, 466, 513 

Reciprocals of Factorials see Factorials, 
Reciprocals of 

Recurring Series 142, 519 

Reeves Electronic Analogue Computer 
(REAC) 501 

Reeves Instrument Corporation Analogue 
Computer (REEVAC) 327, 382 

Relay Computers 1-13, 54, 69-84, 127, 149- 
161, 164-166, 207-208, 216-217, 376- 
378, 545 

Riemann Zeta Function 21, 366, 475-476 

Root-finder 328-329 

Roots, Higher than the Third 340, 466 

Rounding-off Errors 437, 544 


S and T Functions 20, 85-86, 464-465, 467, 
515 

Scattering Function 107, 484 

Secant, Inverse 465 

Secant, Log of 19, 513 

Secant, Natural 19, 93-94, 172, 513, 515, 
556-561 

Selectron 378-379, 444, 544 

Series, Double 399-400 

Series, Summation of 87, 170, 466, 475-476, 
520 

Sine, Hyperbolic 86-87, 93, 170-171, 205, 
336, 463, 465-466, 514 

Sine, Hyperbolic Inverse 517-518 

Sine, Inverse 85, 87, 463, 465 

Sine, Log of 19-20, 85, 87, 355-356, 370, 
463-407, 513-514 

Sine, Natural 19-20, 85-87, 93-94, 172, 335, 
370, 413, 463-467, 485, 513-516, 556- 
561 

Sine Integral 104, 334, 336, 338, 417, 479- 
482, 514, 526 

Sine Integral, Iterated 526 


Slide Rule 19, 58, 134-137, 221, 329, 388- 
389, 446-447, 493, 548 

Snedecor’s Ratio F 363-364, 466 

Spectroscopy 27-29, 255 

Spherical Area and Volume 463 

Spherical Functions 369 

Spherical Harmonics 365 

Spherical Triangles 32, 556 

Spheroidal Wave Functions 99-101, 524-526 

Square Endings 21-22 

Square Root 19, 60, 84-85, 87, 98, 171, 302, 
340, 370, 463, 466, 494, 513-514 

Square Root of Sum of Squares 447 

Squares 19-20, 84-85, 87, 171, 370, 412-413, 
463, 466, 513-514 

Statistics 98-99, 125-126, 133-134, 179- 
180, 300-303, 359-364, 372-373, 415, 
426, 466, 485, 515, 519-521, 524, 546-547 

Stirling Numbers 463 

Stochastic Methods 546-547 

Stokes’ Functions 366-367 

Struve Functions 66, 103, 340, 365-366 

Student's t-Distribution 466, 520, 524 

Sum of Powers of Divisors 23, 177 

Surveying 493 

Symmetric Functions 24, 414 


Tangent, Hyperbolic 86-89, 296, 463, 465- 
467, 514 

Tangent, Hyperbolic Inverse 517-518 

Tangent, Log of 19-20, 85-87, 355-356, 370, 
463-465, 513 

Tangent, Natural 19-20, 85-89, 93-94, 146, 
172, 296-298, 370, 463, 466-467, 513- 
515, 556-561 

Telephony 98-99 

Tetragamma Function 483 

Tetrahedral Numbers 316, 423-424 

Thermodynamics 511 

Theta Functions 89, 104, 229, 235-236, 239- 
241, 244-245, 257, 334, 369 

Toroidal Functions 103, 250 

Transcendental Equations 173-175, 178, 
414, 469, 495-496, 562 

Triangles, Solutions of 556, 561-562 

Trigamma Function 483 

Trigonometric Functions (special) 19-20, 
21, 88, 93-94, 296, 356, 451, 467 

Trigonometric Series 368 

Trigonometric Sums 472 


Universal Automatic Computer (UNIVAC) 
133, 341-350, 378, 382, 540 
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Van der Poi’s Equation 108 
Versed Cosine 94 
Versed Sine 20-21, 94 


Waring’s Problem 24, 423-424 

Wave Equation 106 

Weber Functions 103, 365-366 

Wedge Functions 105 

Weierstrass Function, Elliptic 104, 229, 247- 
250, 254 


Weights and Measures 85, 369, 513 
Whirlwind Computer 444-445, 540 
Whittaker Functions 104 


Zeta Function, Jacobi see Jacobi Zeta Func- 
tion 

Zeta Function, Riemann see Riemann Zeta 
Function 

Zonal Harmonics 514, 521 


NAME INDEX 


There are about 15 names where the determination of appropriate initials was beyond 


our powers. 


Abel, N. H. 104 

Abraham, M. 481 

Abramowitz, M. 477 

Abf’l Hasan Ali 226 

Adams, E. P. 88, 314, 423, 463 

Adams, O. S. 249-250, 263 

Adcock, W. A. 221 

Adriaen Metius 560 

Ageton, A. A. 36, 114-116, 145 

Ahlfors, L. V. 546 

Aiken, H. H. 1, 2, 4, 54, 58, 102, 130, 161, 
164, 185, 188, 210, 311, 327, 367, 377, 
382, 439, 446, 491-492, 499, 517, 545 

Airey, J. R. 25, 102, 107, 116-118, 143-144, 
199, 233, 239, 250, 261, 268, 314, 335- 
336, 338 

Airy, G. B. 337, 339 

Aitken, A. C. 300, 399-400, 469 

Akhmatov, V. V. 420 

Akushskii, I. fa. 60, 221, 299, 314 

Aldis, W. S. 549 

Alexander, S. N. 128, 380, 438, 545 

Alger, P. L. 25 

Al-Khowarizmi 227 

Allan, F. E. 300 

Allen, E. S. 232, 250 

Allen, H. F. 382 

Alphonse X 563 

Alt, F. L. 1, 13, 44, 69, 84, 121, 133, 159, 161, 
179, 380, 382, 440, 446 

Ambrosio, A. 382 

Ambrosio, B. F. 382 

Ames Aeronautical Laboratory 199 

Anan’ev, I. V. 93, 116 

Anderson, E. W. 34 

Anderson, R. L. 360 

Anderson, Ruth K. 382 

Anderson, T. W. 547 

Anding, E. 25, 143 


Andoyer, M. H. 43, 94, 237, 250, 296-297 
426, 515-516, 558 

Andrews, E. G. 2, 133, 179, 445, 490 

Andrews, E. J. 382 

Anger, C. T. 103, 365-366 

Anjema, H. 172, 331 

Anno, S. 382 

Anthonisz, A. 560 

Aquino, F. R. de 35-36 

Archibald, R. C. 20-21, 25-26, 61-62, 66, 
85-87, 89, 92, 94, 102, 105, 117, 142-143, 
145-146, 170, 172-173, 180, 186-187, 
203, 225, 227, 229-230, 255, 270, 272, 
276, 278, 298, 307, 330, 356-358, 366- 
367, 369, 372, 396, 398, 415-418, 425, 
439, 451-452, 463, 467, 475, 482, 493, 
495, 498, 514-517, 520-521, 523, 526, 
549, 559-561 

Archimedes 560 

Arenberg, D. L. 479, 540 

Arenz, R. F. 382 

Arndt, L. 248, 250 

Arnold, H. A. 481 

Arnold, Phyllis 161 

Arnquist, W. N. 382 

Aroian, L. 373 

Arrow, K. J. 382 

Asplund, S. E. 382 

Atkinson, E. 256 

Avery, H. T. 382 

Ayachit, G. R. 360 

Azeredo Rodrigues, O. A. de 114, 529 


Babbage, Charles 53, 143, 210, 215, 396, 
432, 493 

Babbage, R. H. 432 

Babcock, A. B. 329 

Babcock, H. A. 382 

Badellino, M. 25, 415 
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